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ABSTRACT 

This  thesis  presents  an  analysis  of  transonic  potential 
flow  through  an  oscillating  unstaggered  thin  plate  cascade. 
A  collocation  technique  is  used  involving  the  superposition 
of  adjacent  blade  "isolated  foil"  potentials  with  inter- 
ference potentials  of  unknown  strength.   Imposition  of 
flow  tangency  requirements  leads  to  integral  equations  for 
the  unknown  source  distributions  of  the  interference 
potentials.   Results  presented  include  the  interference 
coefficients  for  the  steady  and  oscillating  cases. 

The  steady  case  is  extended  to  the  determination  of  the 
potentials  along  the  blade  and  compares  favorably  with  a 
parallel  solution  using  a  Laplace  transform  approach  to 
the  sonic  wind  tunnel  wall  interference  problem. 
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I.   INTRODUCTION 

Considerable  advances  are  still  being  made  in  the 
development  of  aircraft  turbopropulsion  systems.   In  the 
quest  for  higher  thrust  to  weight  ratios,  one  must  look  to 
higher  stage  loading  as  well  as  significantly  increased 
turbine  inlet  temperatures.   From  the  compressor  and  fan 
design  point  of  view,  this  translates  primarily  into  increased 
tip  speeds.   Moving  toward  better  blade  design,  the  inves- 
tigation of  transonic  and  supersonic  relative  flows  is 
becoming  more  important  not  only  for  increased  performance 
but  also  for  flutter  prediction.   As  shown  by  Pratt  and 
Whitney,  various  flutter  problems  have  already  been  encoun- 
tered [Pig.  1  as  extracted  from  Ref.  1]  and  more  such 
problems  can  be  expected  in  the  future. 

Recognizing  that  the  actual  three  dimensional  flow 
fields  encountered  in  a  transonic  stage  are  complicated  by 
the  effects  of  rotation,  finite  blade  thicknesses  and  cambers, 
boundary  layers,  and  spanwise  Mach  distributions,  it  is 
excessively  difficult  to  model  these  effects  simultaneously. 
One  can  hope  that  eventually  a  complete  analysis  will 
become  practical,  but  at  this  time  the  problems  must  be 
faced  piecemeal.   One  simplified  technique  historically 
used  in  modeling  such  flows  is  the  two  dimensional  cascade. 

The  following  investigation  was  concerned  with  transonic 
cascade  flow.   This  problem  had  previously  been  considered 
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by  Hamamoto  [Ref.  2]  and  Gorelov  [Ref.  3].   Also  of  impor- 
tance is  the  related  problem  of  transonic  wind  tunnel 
interference.   Here  the  studies  of  Sandeman  [Ref.  4]  and 
Murasaki  [Ref.  5]  must  be  mentioned.   For  a  more  general 
exposition  of  this  problem,  reference  is  made  to  Berndt 
[Ref.  6]  and  Goethert  [Ref.  7]. 

Hamamoto  based  his  analysis  on  the  linearized  oscilla- 
tory transonic  flow  equations.   Gorelov  attempted  to  include 
the  effect  of  blade  thickness  by  approximating  the  nonlinear 
transonic  perturbation  equation.   This  approximation, 
however,  is  valid  only  for  low  supersonic  flows.  In  the  present 
investigation,  a  transonic  flow  approximation  was  intro- 
duced based  on  previous  work  by  Oswatitch  [Ref.  8]  and  Teipel 
[Ref.  9].   The  technique  permitted  a  general  formulation 
for  both  the  steady  choke  flow  problem  (in  conjunction  with 
Spreiter's  local  linearization  technique  [Ref.  10])  and 
the  oscillatory  problem.   The  resulting  boundary  value 
problem  was  solved  in  a  manner  similar  to  that  used  by 
Gorelov.   The  problem  formulation,  solution  technique, 
numerical  evaluation  of  the  recurring  integrals,  and 
representative  results  are  described  in  the  following 
sections. 
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II.'  UNSTEADY  TRANSONIC  FLOW  THEORY 


A.   DEVELOPMENT  OP  THE  GOVERNING  TRANSONIC  PLOW  EQUATIONS 


U 


QO 


2C 


Figure  2.   COORDINATE  SYSTEM 


X 


Consider  planar  flow  of  an  inviscid  adiabatlc  perfect 
gas  along  the  positive  X  axis  which  is  perturbed  by  a  thin 
body  (or  thin  bodies)  executing  small  transverse  unsteady 
motions.   The  perfect  gas  equation  is  applicable 


p  =  pRT 


(II-D 


where : 


p  =  pressure 

p  =  density 

T  =  absolute  temperature 

R  =  gas  constant 


Ik 


Furthermore  the  conservation  conditions  for  mass,  momentum, 
and  energy  must  hold.   Denoting  the  velocity  components  in 
the  X  and  Y  direction  respectively  as  U  and  V,  the  conserva- 
tion of  mass  may  be  expressed  as 


i£  +  ^UP)  +  9(VP)  =  0  (II-2) 

at    ax    a  y    u  KX±  C) 


Conservation  of  momentum  is  expressed 


£u  =  au  +  TT3U   vau.    i  ap 
Dt  "  at   uax   ay    p  ax 


Dt   at    ax    ay    p  aY 


Conservation  of  energy  is 


or 


where 


(II-3) 


DT^e  +  ^  =  -  F  tfx(up>  +  Iy(vp^       (ii"4) 


DT  =  °  (JI-5) 


srr-  =  substantial  derivative 

e   =  energy  for  the  individual  particle 

s   =  entropy  for  the  individual  particle 

2     2     2 
W   =  IT  +  V^ 
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These  equations,  of^ course,  exclude  finite  discontinuities 
(shocks)  in  the  flow  field. 

As  is  well  known,  an  important  role  is  played  by  certain 
functions  which  reduce  the  number  of  dependent  variables. 
Most  often  the  velocity  potential  is  introduced.   Its 
existence  depends  on  the  condition  of  irrotationality  of 
the  flow.   For  transonic  flows  over  sufficiently  thin 
bodies  (small  perturbation  assumption!)  it  is  usually 
possible  to  work  with  the  velocity  potential  since  any 
vorticity  introduced  by  occurring  shocks  is  of  at  least 
third  order  in  the  perturbation  velocities.   It  follows 
that  the  velocity  can  be  expressed  as 


U  =  1$  -  * 

dX        vx 


BY    *y 


(H-6) 


The  continuity  equation  (II-2)  may  be  written 


^  8£  -  «     +  *    .  (n-7) 

p  Dt    xx    yy 


Recall  that  the  velocity  of  sound  is 


a2  -  («)  (II-8) 

K  s=constant 


16 


and  that 


Dp_  =  d£  .  DP  =  1_  DP  ,     s 

Dt    dP    Dt     2   Dt*  K         yj 


Therefore  the  continuity  equation  (II-7)  becomes 


$   +  $   +  J— ($  p  +$  p  +p  )  =  0  .   (11-10) 
xx    yy    2'xxyyt' 
a  p 


Integration  of  the  Euler  equations  along  a  streamline  gives 
the  unsteady  Bernoulli  equation 


2 
•  +  W  +  /  S£  ■  p(t)  (ii-n) 

t    ^      p 


1  2 
where  F(t)  =  ^"U^  for  uniform  parallel  flow  at  infinity. 

Taking  the  partial  of  (11-11)  with  respect  to  t,  one  has 


7  pt  ■  -!t[*t  +  £  •  <JI~12) 


Inserting  (II-3)  and  (11-12)  for  P"  ,  P  ,  and  P   in  (11-10) 

x   y       u 

and  utilizing  (II-6)  one  obtains  the  following  partial 
differential  equation  in  *. 

2  2 

A  ft  2$  ft 

^  a2  ;qxx    v±  a2  ;qyy    &2    wxy 
2$        2<(>         , 

"  -T  *xt  "  -f  *yt  -  "7  *tt  =  °         (II"13) 

a        a        a 
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The  usual  way  to  simplify  this  equation  is  to  neglect 
all  product  terms  obtaining 

?     '  "  M  i 

(1-M   )*        +   <£>        -   2-r-4    ,    -  -4*.  .    =    0 

»     xx  yy  am  xt        a2   tt  (n_iin 

oo 

which  is  the  well  known  unsteady  linearized  potential 
equation  valid  for  the  supersonic  and  subsonic  regimes. 
However,  for  the  case  near  M   =  1.  such  a  "naive"  lineari- 
zation  cannot,  in  general,  be  carried  out.   Consider  the 
following  analysis  for  the  steady  case. 

Define  the  perturbation  potential,  4>,  such  that 

$  =  U  X  +  6.  (11-15) 

oo        * 

For  the  steady  case,  the  Bernoulli  equation  (11-11)  is 


2 

-tt-  +   f —  =  constant  .  (II-16) 

2     p 


Now  since 


;^L   =  c   T  =  -^r  R  T 
P     P     Y-l 


and 


a2  =  Cff)   -  YR  T 

H  s 


Equation  (II-16)  becomes 
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W2    a2 

-—  +  — -  =  constant  (11-17) 

2    y-1 


The  constant  is  evaluated  by  noting  that  at  the  critical 
Mach  number 


W*  =  a* 


and  (11-17)  becomes 


2   Y+l   *2    y-1  TT2  ,-,._  -,0\ 

a   =  J-p—  a   -  J-p—  W  (11-18) 


With  (11-18)  and  (11-15)  in  (11-13),  one  has 


<{)      (a2-U2)    +   <f>      a2   =    d>      [(y+1)U   d>   +*~  4    2 
Yxx     «     °°  Yyy  °°        Yxx      '  «>Yx     2      Tx 


+  ^77^  d>    2]   +    2*      [U   A      +    6    d>    ] 
2      Yy  TxyL   «°yy        yxYy 


+   <f>      [(y-1)U  $  +^^4)   2   +  ^^  d>   2]  (11-19) 

yy      '  °°  x      2    Yx  2      Yy 


As  was  mentioned  previously,  under  the  Prandtl  assump- 
tion, one  would  neglect  the  product  terms  on  the  right-hand 
side  (RHS).   However,  for  M  s  1,  a  more  careful  analysis 
is  in  order.   To  this  end  one  introduces  the  Prandtl-Glauert 
transformation: 
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A  =   e<J>(X3Y)  (11-20) 

where   X  =    X 

Y  =   Y  /l-M2  (11-21) 

CO  v  ' 

u 

M  =    —  (11-22) 

oo     a 

CO 

Here,  one  requires  <j>  to  be  of  order  unity,  therefore   e 
assumes  the  order  of  magnitude  of  A. 
Equation  (II-18)  becomes 

T        ,,    .,2,    2  -        (l-M2) a2    =    g2(y+1)U  ^— $- 

ed> —   (l-M    )a      +    z6 —  °°     °°  '  °°yxxYx 

Yxx  °°     °°  Yyy 


+  e3^-  A — A-  +  e3Xo^  (l-M2)  A — A- 
2   YxxYx       2      oo '  YxxYy 


+  2e2d) —  (1-M2)U  A-  +  2e3A — A-A-  (l-M2) 
Yxy     °°  <»Yy      YxyYxYy     oo' 

2 
+  £2(1-m2)(y-1)U„?-?-  H-e^U-M2)2*^ 

+  e3  i^pli-  (l-M2)^-  (11-23) 


Now  the  terms  on  the  RHS  may  be  neglected  for  suitably 
small  e.   To  investigate  the  validity  of  such  a  measure, 

one  may  compare  the  largest  term  on  the  RHS  with  a  term 

2 
from  the  LHS.   This  term  is  the  first  term  of  order  e  ,  the 

3      2     2 

remaining  terms  being  of  order  e   or  e  (l-M  ).   Consider 

CO 
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the   ratio 

U   e(y+D    *- 

CO  '  T  Y 

E   =   -« 5 (11-24) 

•    a     (l-ir) 

CO        v  CO  ' 

It  is  clear  that  the  nonlinear  term  becomes  of  higher 
order  for  the  subsonic  and  supersonic  cases  but  must  be 
retained  for  transonic  flows.   This  consideration  can  be 
generalized  to  unsteady  transonic  flows  thus  leading  to 
the  transonic  small  perturbation  equation  as  the  governing 
equation : 

M^  M 

[1-M  -(Y+l)    tt-  4>    ]4>        +   <$>      -2   —  $    .— %  -$..    =    0 
L        °°      '  U      YxJYxx        ^yy        a        xt      2    Ytt 

co  "  "  oo  n 

00 

(II-25) 

A  rigorous  justification  for  this  equation  was  first  given 

by  Lin,  Reissner,  and  Tsien  [Ref.  11]  and  by  Landahl  [Ref.  12]. 

A  further  simplification  of  this  equation  is  possible  if 
the  transonic  flow  is  sufficiently  unsteady.   Although  this 
can  be  shown  by  mathematical  order  of  magnitude  considera- 
tions in  these  same  references,  it  is  instructive  to  consider 
the  physical  aspects  of  the  transonic  wave  propagation 
problem.   Figure  3  shows  an  acoustic  source  in  transonic 
flow.   It  may  be  seen  that  the  parts  of  the  wave  fronts 
moving  upstream,  referred  to  as  receding  waves,  move  away 
very  slowly  from  the  source.   The  downstream  or  advancing 
waves  move  with  a  velocity  approximately  equal  to  twice  the 
speed  of  sound. 
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Uo» 


Figure  3.   TRANSONIC  WAVE  PROPAGATION 


Assuming  now  a  continuous  spatial  distribution  of  distur- 
bance sources,  it  is  readily  understood  from  Fig.  4  that  the 
receding  waves  will  tend  to  cancel  provided  the  wavelength 
(given  by   2tt  j  a— U  |  /oj  ,  where  w  is  the  circular  frequency  of 
the  disturbance)  is  very  short  compared  to  body  reference 
length,  C.   However,  if  the  sources  oscillate  very  slowly, 


Source  A  and  B  with  Respective  Waves 


The  Resulting  Wave 
Figure  4.   WAVE  CANCELLATION 
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pressure  waves  of  the  same  sign  will  have  sufficient  time  to 
interact  and  possibly  to  build  up  to  finite  amplitudes. 
Obviously,  these  receding  waves  may  create  large  phase  lags 
between  motion  and  pressure  on  the  airfoil  causing  severe 
aeroelastic  instability  problems. 

We  can  draw  the  conclusion,  therefore,  that  sufficiently 
unsteady  transonic  flows  can  be  analyzed  by  linearized 
theory,  Eqn.  (11-14),  whereas  only  mildly  unsteady  transonic 
flows  must  be  investigated  using  Eqn.  (11-25).   The  condition 
for  linearization  is  generally  stated  as  [Ref.  12] 

K  »   I1'"  mlocaJ 

where  |l  -  ^t or AT  I  is  the  larSest  deviation  of  the  Mach 
number  from  unity  in  the  flow  and  K  is  the  reduced  frequency 
(K  =  wC/Uj. 

B.   DEVELOPMENT  OF  THE  BOUNDARY  CONDITIONS 

To  completely  describe  the  problem,  proper  boundary  and 
initial  conditions  must  be  specified.   The  most  important 
ones  are: 

1.  The  flow  on  the  surface  of  the  body  must  be  tangent 
at  any  instant  of  time. 

2.  Sommerf eld's  radiation  condition  must  be  satisfied 
at  any  instant  of  time,  i.e.,  waves  propagate  away  from 
sources  of  disturbance. 

Further,  certain  edge  and  wake  conditions  must  be 
satisfied,  e.g.,  smooth  flow  at  trailing  edges  must  be 
established  (Kutta  condition). 
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The  general  flow  tangency  condition  can  be  stated  as 
follows.   If  the  equation  of  the  surface  of  a  body  v/hich  is 
moving  in  a  time  dependent  manner  is  given  by 

p(x,y,t)  =  o,  (11-26) 

then  the  flow  tangency  condition  requires  that  the  total 
derivative  of  the  function  P  is  zero,  i.e.. 


~=F.  +UF   +VP  =0.  (II-27) 

Dt    t      x      y 


Consider  an  airfoil  that  is  placed  in  a  Cartesian  coordinate 
system  in  such  a  way  that  if  the  airfoil  is  restricted  to 
small  amplitude  motions  out  of  the  Y  =  0  plane,  then  its 
upper  and  lower  surface  can  be  described  by 


P   =  Y  -  Y  (X,t)  -  0 
u        u   ' 


PT  =  Y  -  YT (X,t)  -  0 


(11-28) 


Therefore,  from  the  boundary  condition  one  has 


V  =  Y    +  U  Y  at   Y  =  Y„(X,t) 

Ut  x                 U 

V  =  YL   +  U  YL  at   Y  =  YL(X,t) 

t  x 


(11-29) 


One  now  introduces  the  small  perturbation  assumption, 

$  =  U  X  +  <J>  (11-15) 


00 


where  the  disturbance  velocity  components  are 

u  =  4) 


x  (11-30) 

v  =  d> 
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such  that 

U  =  UM  +  u  (II-3D 

V  =  v  (II-32) 

and  u.v  <<  U 

'  CO 

In  (11-29)  the  term,   u  Y    ,  can  then  be  neglected  compared 

x 

to  the  much  larger  term,  UWY    .   The  boundary  conditions 

x 
then  simplify  to 


v  =  Yu   +  U00Yu     for  Y  =  Yu(X,t) 
u        x 


and  v  =   YL      +   UTOYL  for  Y   =   YL(X,t). 


(II-33) 


For  small  amplitude  oscillations,  (11-33)  can  be  further 
simplified  by  expressing  the  flow  tangency  condition  in  the 
plane  Y  =  0.   Taylor  series  expansion  gives 

v(X,Yu,t)  =  v(X,0+,t)  +  Yu  3^X>0  ?t)  +.... 

and  v(X,YL,t)  =  v(X,0",t)  +  YL  3%yX?°  ill    + 

Using  the  same  argument  as  above,  the  product  terms  can 
be  neglected  so  that  the  flow  tangency  condition  can  be 
stated  as 

v  =  Y    +  U   Y  for  Y  =  0+ 

t        x  (11-34) 

v  =  YT   +  U  YT  for  Y  =  0" 

Lt    °°  Lx 

For  a  thin  foil  performing  a  time  dependent  motion  about 
the  Y  =  0  plane,  the  airfoil  thickness  terms  can  be  dis- 
regarded by  similar  arguments.   The  above  conditions  become 
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v  "  H  +  U»  H         for   Y  =  0     (11-35) 

where  h(X,t)  describes  the  time  dependent  chord  line  of  the 
blade. 

C.   THE  DERIVATION  OF  THE  UNSTEADY  EQUATIONS 
Define  the  following  nondimensional  terms. 


X 

= 

X 

c 

Y 

= 

Y 
C 

t 

= 

t    U 

00 

c 

? 

= 

* 

u  c 

00 

V 

= 

V 

h 

= 

h 
C 

(11-36) 


where  C  is  the  chord  of  the  foil.   Equations  (11-25) 
and  (11-35)  may  now  be  written  in  nondimensional  form: 

[1    _    M2    _    (y+1)M2^"-]^ — ■  +? 2M2?-r-  -    M2for  =    0 

L  °°        VI      '   a>Yx     xx       yy  °°  xt  °°   tt 


(11-37) 


and  v  =   hr-  +   h-  (11-38) 

X>  X 


Since  the  main  interest  of  this  thesis  centers  on  the 
determination  of  the  aeroelastic  stability  of  a  foil,  only 
small  perturbations  from  the  steady  state  need  be  studied. 
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If  one  restricts  oneself  to  such  small  amplitude  perturba- 
tions, it  is  possible  to  write  the  perturbation  potential  as 
the  sum  of  a  steady  state  and  a  time  dependent  potential. 
Also,  in  the  case  of  a  harmonic  oscillation,  one  may  describe 
the  time  dependent  quantities  as  the  product  of  a  complex 
amplitude  function  and  a  harmonic  time  factor.   Therefore 
one  defines  the  potentials  <j>,  and  \p   such  that 

4(X,Y,t)  =  4)1(X,Y)  +  iKX,Y,K)e1KF  (II-39) 

and  the  foil  chord  function  and  vertical  velocity  are 

h(X,t)  =  h1(X)  +  h(X,K)elKt  (11-40) 

v(X,t)  =  Vl(X)  +  v(X,K)eiKF  (11-41) 

where        K  =  U~  (11-42) 

CO 

It  may  be  shown  that  (11-37)  becomes 

[1-M2-M2(y+1  )*-,-]*•-, —  +  <J>n —  =  0  (11-43) 

L   oo  «»^  i   yyixJVlxx    lyy 

which  is  the  steady  state  equation,  and  the  unsteady  equation 

is 

[1-M2]^ M2(y+D—  (<f>-,-*J>-)  +  M2K2ip 

°°   XX    °°  '   \—        lxTx      °° 

-2iKM2ip-  +  <h»  =   0  .  (11-44) 

w  x     yy 

The  boundary  equation  (11-38)  becomes 

v-,(X)  =  h,-(X)  (11-45) 

1       lx 

and     v(X,K)  =  h-(X,K)  +  iKh(X,K)  (11-46) 
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On  the  foil  surface  equation  (11-43)  must  satisfy  the 
boundary  equation  (11-45);  (11-44)  must  satisfy  (11-46). 

For  convenience,  the  overbar  and  hat  will  now  be  dropped 
and  the  steady  and  unsteady  equations  with  boundary  equa- 
tions become 


[1-M2  -  M2(Y+l)<k  ]$_    +  <f>,    =0     (II-47) 
oo    c»v   /Ylx  Ylxx    lyy 


subject  to 


v1  =  hlx   9i      at  Y=0  (11-48) 


and 


[1-M2]i|)   -  m2(y+1)I— (4>n  *  )  +  M2K2^ 
00  Txx    <*>' '   '8x  Ylxyx     °° 


-  2iKM2ij^  +  *   -  0  (II-49) 


subject  to 


v  =  h   +  iKh   ,   at  Y=0  .  (11-50) 


Note  that  (11-47)  is  a  nonlinear  equation  in  <J>,.   Equation 
(11-49)  is  linear  in  ty  but  the  existence  of  the  variable 
coefficient  <t>,   considerably  complicates  its  solution. 

JL  X 

D.   A  BRIEF  REVIEW  OF  POSSIBLE  APPROACHES  TO 
UNSTEADY  TRANSONIC  FLOW  ANALYSIS 

The  following  is  a  presentation  of  various  approaches  to 
the  solution  of  equation  (11-49). 

1.   Method  of  Teipel  and  Hosokawa 

An  approximation  of  equation  (11-49)  was  proposed 
by  Teipel  [Ref.  9]  and  Hosokawa  [Ref.  13]  for  oscillating 


28 


airfoils  in  sonic  flow  (M^  =  1).   For  steady  sonic  flow, 
Oswatitch  [Ref.  8]  introduced  a  linearized  approximation  by 
setting 

(Y+lHlxx  =  X  =  constant,   A  >  0  ,     (II-5D 
allowing  (11-43)  to  be  written  in  the  easily  solved  form 

*lyy  "  X   *lx  "  °  •  ^-^ 

Teipel  and  Hosokawa  generalized  this  procedure  to  unsteady 
sonic  flow  by  approximating  equation  (11-49)  with  (neglecting 
the  term  in  ^Yv) 

-X  *x  +  ^yy  -  2iK^x  +  K2^  -  0  .        (II-53) 

Defining  Q 

Q  =  X  +  2iK  ,  (11-54) 

one  has  from  (11-46) 

^yy  -  Q  ^x  +  K2  i>   =  0  .  (11-55) 

This  equation  allows  the  study  of  sonic  flow  past  a 
single  oscillating  airfoil  in  unbounded  flow  over  the  entire 
frequency  range.   It  contains  the  steady  solution  (Oswatitch 
approximation)  as  a  special  case.   For  large  K  (Refer  to 
page  23)  the  term  containing  <j>,    becomes  negligibly  small; 

J.  AA 

therefore,  retention  of  X  at  high  reduced  frequencies 
becomes  unnecessary.  Thus,  when  one  takes  X  =  0,  the 
Teipel-Hosokawa  theory  also  contains  the  linearized  theory. 
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2.  Local  Linearization  Technique  for 
Steady  Transonic  Flow 

It  was  shown  that  the  Oswatitch  approximation  of  the 

steady  transonic  small  perturbation  equation  is 

*lyy  ~  X  .*lx  =  °  •      X  =  constant  >    °  •       (H-52) 

The  solution  to  this  equation  for  sonic  flow  past  a  single 
airfoil  is 

X 

♦nCX.Y)  -  -  -+=,       fS|.  ^Z§T  ds  .       (11-56) 

\fllX       -1  /x-s 

Of  particular  interest  is  the  velocity  on  the  profile  which 
is  given  by 

4lx(X}0)  =  -  JL  ^  /X  Xt|i   ds  .  (11-57) 

1X         /irV  ax  -1  /X-S 

Recalling  the  definition  of  A,  (11-51),  Spreiter  argued  that 
the  original  variable,  $1       ,  could  be  reintroduced  into 
equation  (11-57)  thus  linearizing  the  equation  locally. 
This  leads  to  an  ordinary  nonlinear  differential  equation  in 
(|>1   which  can  be  solved  analytically.   The  agreement  with 
experiment  is  remarkably  good  [Ref.  10].   Its  application 
to  the  choked  wind  tunnel  wall  interference  problem  was 
given  by  Sandeman  [Ref.  4],   It  is  expected  that  an  extension 
of  this  technique  to  the  unsteady  case  may  prove  fruitful. 

3.  WKBJ  -  Method 

Rubbert  [Ref.  1.4]  applied  the  .WKBJ  method  to  this 
problem.   Differentiating  the  axisymmetric  counterpart  of 
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equation  (H-I19)  for  Mtt  =  1,  multiplying  by  <J>   and  applying 
the  Hankel  transform,  the  equation  can  be  reduced  to  the 
form 

^xx  +  N2q(x,N)Y  =  0  (11-58) 

where  N  is  a  large  parameter.   Unfortunately  no  inversion 
of  the  resulting  transform  integrals  has  been  found. 

k .   Method  of  Geometrical  Acoustics 

A  third  approach  has  recently  been  proposed  to  obtain 
the  thickness  effects  for  oscillating  wings  in  transonic 
flow  [Ref.  15].   This  procedure  begins  with  a  given  non- 
uniform steady  transonic  flow  field  which  has  been  deter- 
mined either  analytically  or  experimentally  (In  the  latter 
case  an  analytical  expression  must  be  obtained  which  matches 
the  measured  pressure  coefficient.).   Superimposing  upon 
this  flow  field  small  unsteady  perturbations,  the  path  and 
transmission  time  of  these  acoustic  perturbations  are 
determined  by  the  methods  of  geometrical  acoustics.   Again, 
only  preliminary  calculations  have  been  carried  out  to  date. 

5 .   Unsteady  Transonic  Flow  with  Shock  Waves 

Typically,  transonic  flows  around  an  airfoil  exhibit 
locally  supersonic  regions  which  are  terminated  by  shock 
waves.   Often,  the  interactions  of  small  unsteady  perturba- 
tions with  these  shock  waves  may  be  responsible  for  aero- 
dynamically  induced  instabilities.   Eckhaus  [Ref.  16]  made 
the  problem  mathematically  tractable  by  assuming  a  normal 
shock  extending  from  the  surface  to  infinity  with  uniform 


31 


steady  supersonic  flov;  in  front  and  uniform  steady  subsonic 
flow  behind  the  shock.   Performing  a  small  perturbation 
analysis,  i.e.,  superimposing  small  time  dependent  perturba- 
tions upon  this  idealized  steady  transonic  flov;,  he  was  able 
to  obtain  equations  for  the  interaction  of  the  sound  waves 
with  the  shock.   However,  only  a  few  numerical  results  were 
obtained. 

For  the  following  study,  the  Teipel-Hosokawa 
approach  was  used. 
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III.   TRANSONIC  UNSTEADY  CASCADE  FLOW  PROBLEM 


A.   PROBLEM  FORMULATION  FOR  SONIC  FLOW  PAST 
OSCILLATING  UNSTAGGERED  CASCADES 


■*- 


2b 


■>"  £ 


Figure  5.   UNSTAGGERED  CASCADE 

Consider  an  infinite  cascade  of  thin  blades  arranged  as 
shown  in  Figure  5.   Blade  chord  is  2b  and  blade  spacing  is 
P.   Subscript  zero  referes  to  the  reference  blade  and 
subscript  one  refers  to  the  adjacent  blade  in  the  positive 
ordinate  (ri)  direction.   The  reference  coordinate  system  is 
fixed  at  mid-chord  with  the  abscissa  (£)  in  the  direction 
of  flow. 

Two  n  coordinates  are  specified  with  n   originating  on 
the  reference  blade  and  ru  on  the  adjacent  blade  such  that 


nx  =  nQ  +  P 


(Iii-D 


Now  the  following  nondimensional  variables  are  defined: 


x  =  £ 

*     b 


(IH-2) 
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Ym  -  •£■     >   m-0,1  (III-3) 

R  -  |  (111-4) 


As  the  governing  equation,  the  Teipel  approximation  is 
used , 

*  -  Q  *x  +  K2  ip  =  0  .  (11-55) 

This  equation  must  be  solved  subject  to  the  following 
boundary  equations: 

1.  Zero  disturbance  upstream  of  the  cascade. 

2.  Flow  tangency  on  the  blades. 

As  shown  in  Section  II-B,  this  latter  requirement  can  be 
expressed  mathematically  as  follows.   Along  the  reference 
blade  equation  (11-50)  may  be  written 

*y(X,Yo=0)  =  vo  =  ^+  iKh0   '        flIII-5) 

and  along  the  adjacent  blade 

9h, 

*  (X,Y0=R)  -  v1  =  ^p  +  iKh1  .        (III-6) 

To  solve  this  boundary  value  problem,  the  potential 
function  ty   is  written  as  the  summation 

♦<XfY0)  =  *°(X,YQ)  +  ^0(X,YQ)  +  *°(X,Y1) 

+  $1(X,Y1)  (III-7) 

where  $   and  $?  are  the  isolated  blade  solutions  for  the 
o      1 

reference  and  adjacent  blades  and  $   and  $>   are  interference 

0  o      1 

potentials. 
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Such  an  approach  has  previously  been  proposed  for  the 
subsonic  and  supersonic  cases  by  D.  N.  Gorelov  [Refs.  3  and 
17].   Also  refer  to  Samoylovich  [Ref.  18]. 

Inserting  (III-7)  into  (III-5)  and  (III-6)  and  realizing 

*my(X>0)  =  Vm  •   m=0>1  ■ 
one  obtains  for  the  reference  blade, 

»   (X,Y  =  0)  +  $?  (X,Yn=  -R) 

oyv  »  o        lyv  *  1         ly 


*ov(X>V  0)  +  *ly(X>Yl=  "R)  +  *lv(X'Yl=  "R)  =  ° 


(III-8) 
and  for  the  adjacent  blade, 


oy   '  o        ly 

(III-9) 


$°  (X,Y  =  R)  +  $   (X,Y  =  R)  +  <K  (X,Y  =  0)  =  0  . 
oy   '  o        oy   '  o        ly   »  1 


The  isolated  blade  solution  is  already  known  from 

Teipel  [Ref.  9]  as  2  y2 

EXP[^(X-S)    -   «   jfr 

^'^  -  "/fe  -I  V°(S)J 7^ dS 

(111-10) 

2 
and  x  EXP[^(X-S)    -   §  x^] 

*?(X,Yt)    =   —        /      v,  (S) dS 

11  vW     -1         1  /x^ 

(III-ll) 


The    interference   potentials  must   have   the    same   form  but 

contain  unknown   source   distributions   u„   and   uu  :    thus, 

o  Is       p       ' 

2  Y 

.,  x  EXP[^(X-S)    -  §■  ^4] 

$    (X,YJ  =   -  —        /      u    fS)    S 5    A"b      dS 


(111-12) 
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2 

and  y2  n     Y, 

-         .      x'  EXP[^r(X-S)    -   jf  T+rl 

•nCX.Y,)   =  -—       /     u-.CS)   Q     . , =L-^-  dS.    (111-13) 

1  -1  /ttQ      -1        x  v/X-S 


These  source  distributions  are  approximated  by  the  power 
series  as  follows: 


L-l      M 

ii  =   I  PH3T  (111-14) 

m        N£Q  mN 


and 

vm  =  I      6mN  ^  '  (HI-15) 


Inserting  these  approximations,  (III-l^)  and  (111-15), 
into  equations  (III-8)  and  (III-9)  results  in 

ll0{    PoNXN  +  I  /?  (°1N+P1N)_/X7^5  EXP[4(X-S)-|||]d^ 

=  0     (111-16) 


for   the   reference  blade   and 

ll0{Vl^+  ?     l/f(60N+P0N'_{X^iEXPt4(X-S)-"  X^dS} 

=   0  (111-17) 

for  the  adjacent  blade. 
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S—  X 


Figure  6.   BOUNDARY  CONDITIONS  FOR  TT  PHASE  LAG 


The  coefficients  in  (111-15)  may  be  evaluated  for  the 
case  of  simple  oscillation,  no  translation,  and  tt  phase  lag 
between  blades  as  follows.   One  arbitrarily  sets  the  angle 
of  maximum  amplitude  for  the  reference  blade  at  a   =  -1 
[Refer  to  Fig.. 6].   Then  the  angle  of  the  adjacent  blade  is 
a-,  =  1.   For  small  a,  equation  (III-5)  may  be  written 


*y(X,Y0=  0)  =  -1  -  iKX 


(111-18) 


and  (III-6), 


IPyU.Y^    0)     =    1    + 


iKX 


(111-19) 
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Hence,  the  coefficients  are 


G    =  -1 
oo 

0  ,  =  -iK 
ol 


.  (111-20) 

lo  =   1 


and  0.,-.  =   iK 


Similarly,  for  in  phase  oscillation,  the  coefficients  are 

6    =  0-   =  -1 

00  lo 

(111-21) 

and  6  n  =  e..n  =  -iK  . 

01  11 

Satisfying  boundary  equations  (III-16)  and  (111-17)  at 
L  points  along  the  blades  generates  a  set  of  2L  equations 
in  2L  unknown  interference  coefficients.   The  coefficients 
in  hand,  ^(X,Y  )  in  (III-7)  becomes  a  known  function. 

A  computer  program  for  the  solution  of  this  problem  in 
the  tt  phase  shift  case  is  developed  in  Appendix  A.   Some 
sample  results  are  presented  in  Table  II  for  the  steady  case 
and  Table  VII  for  the  unsteady  case. 


B.   EVALUATION  OP  THE  INTEGRALS  OCCURRING 
IN  EQUATIONS  (III-16)  and  (111-17) 

The  integrals  occurring  in  (III-16)  and  (111-17)  are  of 

the   form 

In  =  X  o&f*  EXpr4(x-s)  -  §  DPs7]  dS    •        (III-22) 

By  making  the  substitution,  U  =  X-S,  one  obtains 

jn  =  £x+li^r  EXPC4  u  "  f  ir]  dU  (III"23) 
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Recalling  Q  =  A  +  2iK,  I„  is  complex  and  must  be  resolved 
into  real  and  imaginary  parts  for  evaluation. 

The  exponential  term  in  (111-23)  may  be  written 


rrprK2  TT    Q  R2n   _  FYpr  AK2U      AR2 
EXP[-q  U  -  f  -y-]   -  EXPL-2— ~~2    ~  Tip 


Q      "  U  V  +  W 


A^  +  4lC    2U 


One  now  defines  the  following  terms 
A-    *2 


X2  +  'IK2 


B  =  i£- 


E  =  *£ 


Equation  (111-24)  may  be  rewritten, 

EXP[~  U  -  |  ~_]  =  EXP[XAU  -  |]EXP{i[-2KAU-|]}   (111-25) 


Separating  (111-25)  into  its  real  and  imaginary  parts,  one 
obtains  the  following  two  fundamental  integrals  from 
(111-23): 

XNR  =  f      (X~^   EXP[AAU  -  |]  C0S[2KAU  +  |]  dU    (111-26) 

and  IWT  =   (*   ^Mrr   EXP[AAU  -  £]  SIN[2KAU  +  £]  dU    (111-27) 

N  1      Q  ^        rj  3/^  U                       U 

Where    JU    =  TNR  "  1TNI  (IH-28) 
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Expanding    (X-U)    in  the   finite  binomial   series, 
M        N+l    f    rVJ.-l   ...    VN-J+1    TTJ-1 

equations  (111-26)  and  (111-27)  become 

N+l    f    ,  \J-1   KTI    VN-J+1      X+lr    T    «r  u 

hn  •  J3  THrTr-rroLrr  „/   {"   *  ™AU  - |i  ■ 

•  C0S[2KAU  +  |]|  dU 

(111-30) 


and 


N+l  ,  ,  vJ-1 N-J+l   X+l 


N+l  ,  ,  nJ-±  m.  yn-j+±   a+u    j 

•  SIN[2KAU  +  |]|  dU 

(III-3D 

The  fundamental  integrals  in  (111-30)  and  (111-31)  do 
not  lend  themselves  to  analytical  integration.   Numerical 
integration  techniques  may  be  used  but  care  must  be  exer- 
cised in  the  vicinity  of  the  lower  limit.   Here  the  integrand 
is  singular  when  J<5/2.   The  problem  is  further  aggravated 
by  the  rapidly  oscillatory  character  of  the  functions 
SIN[XAU  -  jj]  and  COS[XAU  -  |] . 

Three  specific  cases  are  posed  for  the  treatment  of 
these  integrals. 

1.   The  Steady  Case 

Setting  K  =  0  reduces  (111-30)  and  (III-3D  to  the 
following: 
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N+l  ,  ,  NJ-1  „,  ,rN-J+l   X+l 


hn   "   Z   W|W  J "V-f  EXP[-  |]  dU   (111-32) 


and 


%  "  ° 


The  numerical  integration  of  this  case  is  straight- 
forward.  The  exponential  term  forces  the  integrand  to  zero 
near  the  origin  permitting  one  to  neglect  that  portion  of 
the  integral  very  near  U  =  0.   Equation  (111-32)  may  also  be 
evaluated  by  expressing  it  in  terms  of  the  error  function. 
(Note:   complete  agreement  between  the  two  integration 
techniques  has  been  obtained  thus  verifying  the  numerical 
integration  procedure  for  the  steady  case.) 
2.   The  Oscillatory  Case  with  X  =  0 

In  the  unsteady  case  for  X  =  0,  equations  (111-30) 
and  (111-31)  contain  the  basic  integrals 


X+l   T_r  p 

/    U   *  C0S[2KAU  +  £]  dU  (111-33) 

0  U 


and  |+1  UJ  *  SIN[2KAU  +  |]  dU.  (111-34) 

0 


To  avoid  the  previously  mentioned  difficulties  near 
the  lower  limit,  the  integration  must  be  accomplished  in 
two  steps,  i.e., 


4l 


/         UJ~^  s?nC2KAU   +    l]dU   =    ^   U<J~^  SINC2KAU   +    l]dU    + 


+      /      UJ~"2   g^[2KAU   +   |]dU. 
a 

(111-35) 

For  judicious  selection  of  (a),  the  second  integral 
in  the  RHS  may  be  numerically  integrated  in  a  straightforward 
manner.   The  integral  on  the  limits  0  to  (a),  which  will 
henceforth  be  referred  to  as  the  correction  integral,  may  be 
approximated  for  sufficiently  small  (a)  as  follows. 

Recall  the  following  relations: 


C0S[2KAU  +  |]  = C0S[2KAU]C0S[~]  -  SIN[2KAU]SIN[|]      (111-36) 

and 

SIN[2KAU  +  |]  = SIN[2KAU]C0S[jj]  +  C0S[2KAU]SIN[|]      (111-37) 


and  the  truncated  series  expansions 


2 

C0S[2KAU]=1  -  (2KAU)  (111-38) 

and 

o 

SIN[2KAU]  = 2KAU  -  (2KAU)   g  (111-39) 


Equations  (III-38)  and  (111-39)  are  valid  for  2KAU  suffi- 
ciently close  to  zero.   One  may  now  write 


H2 


/   UJ  *   C0S[2KAU  +  fhdU  - 
0  U 


!      UJ~*  COSCfhdU  -  (2KA)/   UJ~^'  SIN[ff]dU 
0  u  0  u 

•  (2K^   /a  UJ'^  COS[|]dU 
d        0  u 


+  (2K^   /   UJ+iE  SIN[|]dU  (111-40) 


and 


/   IT  *  SIN[2KAU  +  £]dU  = 
0  u 

/   UJ"*  SIN[~]dU  +  (2KA)  /   UJ"SCOS[|]dU 
0  u  0 

-  (2K^   ^  UJ~*  SIN[~]dU 
d         0 

(2KA)3  /   UJ+*  COS[£]dU  .  (111-41) 
S-  0         U 


In  theory  integrals  of  the  form 


ra  TTm/2  COS  rE-,,TT 
/   U    SIN  [U]dU 


may  be  expressed  in  terms  of  the  Fresnel  integrals.   Hence, 
integration  of  (III-JJO)  and  (III-^ll)  can  be  accomplished  in 
"closed"  form.   (In  actuality  there  is  a  computer  accuracy 
limit  on  this  technique  as  described  in  Appendix  A-IV. ) 
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Verification  of  this  technique  for  J  =  1  was  accomplished 
by  checking  it  against  a  known  analytical  function  as 
follows : 


l-k   S?N[2KAU  +  §^U  +  /"-L  °°*[2KAU  +  §]dU 

0    U   *  a   U  * 


00 


1  COSr^nTT       ,       E 


/      -i-  ^2KAU   +   ^dU  (111-42) 

0      \j  %    b   N  U 

where  the  correction  integral  was  evaluated  using  the  Fresnel 

integral  technique,  the  second  integral  was  evaluated  using 

numerical  integration  techniques,  and  the  integral  on  the 
RHS  was  evaluated  using  the  following  relation: 

00  __^__  

For  the  case  K  =  0.01  and  a  =  0.01,  Table  IA  presents 
the  comparison  of  the  value  of  (111-43)  with  that  obtained 
from  (111-42)  for  various  R's.   (The  numerical  integration 
was  truncated  at  an  upper  limit  of  500  and  the  convergence 
criterion  EPS  for  subroutine  ZSUM  was  set  at  1.0E-4.) 

Table  IB  displays  the  computed  integral  using  (111-35) 
with  an  upper  limit  of  X+l  =  2.0  for  various  values  of  (a). 
(The  convergence  limit  placed  on  ZSUM  was  set  at  1.0E-5.) 
Table  variables  are  listed  below. 
DK  =  K 
DL  =  X 
RHO  =  R 
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AL  =  a  (Correction  integral  upper  limit.) 

SCOR  =  Value  of  correction  integral  (III-lJl)  for  J  =  1. 

CCOR  =  Value  of  correction  integral  (111-40)  for  J  =  1. 

SUMI  =  Value  of  numerical  integration  term  from  (a)  to 
specified  upper  limit  for  the  SIN  integral. 

SUMR  =  Value  of  numerical  Integration  term  from  (a)  to 
specified  upper  limit  for  the  COS  integral. 

SINT  =  Value  of  (111-35)  SIN  integral  for  specified  upper 
limit  and  J  =  1. 

CINT  =  Value  of  (111-35)  COS  integral  for  specified  upper 
limit  and  J  =  1. 

SSIN  =  Value  of  (111-43)  SIN  integral. 

CCOS  =  Value  of  (111-43)  COS  integral. 

3.   The  General  Case,  K/  0  and  X  ^  0 

In  the  apparent  absence  of  a  suitable  correction 
integral  for  this  case,  one  must  evaluate  the  full  equations 
(111-26)  and  (111-27)  by  numerical  methods.   The  difficulty 
close  to  the  lower  limit  is  alleviated  to  some  extent  by  the 
exponential  term.   The  exponential  function  drives  the 
amplitude  of  the  integrand  to  zero  more  rapidly  than  a 
negative  power  of  U  can  amplify  the  highly  oscillatory 
behavior  of  the  transcendental  functions.   For  this  reason, 
as  in  the  case  of  K  =  0,  one  may  justify  neglecting  that 
portion  of  the  integrals  very  near  zero. 

Although  no  results  are  presented,  the  computer 
routine  developed  in  Appendix  A  is  designed  to  carry  out 
such  a  case. 
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TABLE  IA  VERIFICATION  OF  CORRECTION  INTEGRAL 


R 

SCOR 

SUMI 

SINT 

SSIN 

0.5 

34.62 

1.016 

35.63 

35.57 

1.0 

14.45 

3.45  8 

17.91 

17.85 

1.5 

4.968 

7.032 

12.  00 

11.94 

2.0 

-1.  114 

10.16 

9.045 

8.987 

R 

CCOR 

SUMR 

CINT 

CCOS 

0.5 

15.48 

19.80 

35.28 

35.32 

1.  D 

-1.782 

19.34 

17.56 

17.60 

1.5 

-5.797 

17.45 

11.65 

11.69 

2.0 

-4.490 

13.19 

8.696 

8.736 

COMPARISON  OF  VALUES  COMPUTED  USING 
EQUATION  ( I  I  1-43)  WITH  THOSE  OBTAINED 
USING  (111-42)  FOR  AN  UPPER  LIMIT  OF  500 
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TABLE    IB    VERIFICATION    OF    CORRECTION    INTEGRAL 


DK=    0.500                  DL^  0.0 

SCOR=       1.57005D    00  CCOR; 

SUMI=       3.96851D    00  SUMR; 

SINT=       5.53856D    00  CINT^ 


RHO  =    0.500 
1.67114D-01 

3.24013D    00 
3.40724D    00 


AL=    0.010 


DK=    0.500  DL=    0.0 

SCOR=       1.66484D    00  CCOR  = 

SUMI=       3.873720    00  SUMR  = 

SINT=       5. 53856D    00  CINT^ 


RHO=  0.500 
2.65483D  00 
6.06207D  00 
3.40724D  00 


AL=    0.050 


DK  =    0.500  DL=    0.0 

SCOR=       3.74310D    00  CCOR 

SUMI=       1.79545D    00  SUMR= 

SINT  =       5.53856D    00  C I  NT 


RHO=  0.500 
■1. 12801D  00 
4.53526D  00 
3. 4072 4 D  00 


AL=    0.100 


DK=    0.500                  DL=  0.0 

SCOR=       4.62026D    00  CCOR; 

SUMI=       9.18297D-01  SUMR= 

SINT=       5.53856D    00  CINT' 


RHO-    0,500 
4.96270D-01 
2.91097D    00 
3.40724D    00 


AL=    0.200 


DK=    0.500                  DL=  0.0 

SCOR=       4.878260    00  CCOR; 

SUMI=       6.602940-01  SUMR; 

SINT=       5.53856D    00  CINT^ 


RHO=    0.500 
1 .27508D    00 
2.13216D    00 
3.40724D    00 


AL=    0.300 
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C.   CALCULATION  OF  THE  POTENTIAL  FUNCTION  AND 
THE  PRESSURE  COEFFICIENT 

The  potential  function  is  obtained  from  equation  (III-7), 

*(X,T0)  =  *°(X,YQ)  +  ^0(X,Yo)  +  ♦J(XfY1)  *  $1(X,Y1) 

(III-7) 
The  interference  potentials  $   and  4>   may  now  be 
computed  using  the  previously  determined  coefficients  (V   N, 
m  =  0,1;  N  =  0,...,L-1)  in  equations  (111-12)  and  (111-13). 
Equation  (III-7)  may  be  rewritten  in  terms  of  the  coeffi- 
cients as  follows: 


+  7==  (e  +p  )  /  -A—  exp[^(x-s)4  r5rqr]ds 

(111-44) 

The  resulting  integrals  may  be  evaluated  using  the  same 
techniques  utilized  in  the  coefficient  derivation. 

The  pressure  coefficient  is  obtained  from  the.  relation 

Cp  =  -2(*  +  IK*)  .  (111-45) 

By  placing  the  partial  with  respect  to  X  of  equation  (111-4*1) 
into  (111-45),  one  may  obtain  the  pressure  distribution  at 
any  point  in  the  flow  field  as  well  as  on  the  blade  surfaces. 
Further,  the  force  and  moment  coefficients  are  easily 
obtained  by  integration  over  the  blade  chord. 
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D.   EVALUATION  OP  THE  STEADY  STATE 
POTENTIAL  ALONG  THE-  BLADE 

The  computation  of  the  steady  state  potential  along  the 

upper  surface  of  the  reference  blade  follows  from  equation 

(111-4*0  evaluated  at  K  =  0,  YQ  =  0,  and  Y-,  =  -R.   The 

development  which  follows  is  for  the  tt  phase  shift  case 

where  [From  Eqns.  (111-20) .  "J 


6    =  -1 
oo 


0lo  =   1 


(111-46) 


Equation  (111-44)  becomes 


WX'0)    " F^ 

/ttX 


1       L~1 

+  "7=i     I      ■ 
/ttV  N=0 


X 


X      SN 


X     K 


x    exp[-  if  t%1 


I     -1-  ds  +  J 

-1  /x-s'  -1  /x-s1 


IN 


-1      \/x-S 


-  EXP[-   \  |^]dS 


X        „N 


oN 


-i     /x~^ 


dS 


(111-47) 


The  computer  program  is  documented  in  Appendix  B-I.   The 
resulting  potential  values  are  compared  with  those  obtained 
from  the  Sandeman  wall  interference  solution  of  Section  III-E 
in  Tables  III  and  IV. 
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E.   LAPLACE  TRANSFORM  SOLUTION  FOR  THE  SONIC  WIND 
TUNNEL  WALL  INTERFERENCE  PROBLEM 

A  parallel  development  to  the  collocation  technique  is 

found  in  the  Laplace  transform  solution  of  the  sonic  wind 

tunnel  wall  interference  problem  [Ref.  4], 


Virtual  Blade 


R  =  2d 


/////////////////////////// 


77  ////////// //////////////// 


Figure  7.   SONIC  WIND  TUNNEL 


Figure  7  shows  the  arrangement  of  a  blade  facing  a  solid 
wall.   Due  to  symmetry,  the  solution  for  the  flow  between 
the  blade  and  the  wall  will  be  the  same  as  that  which  would 
exist  between  the  blade  and  the  virtual  blade  beyond  the 
wall.   Hence  it  is  the  same  solution  expected  of  an  unstag- 
gered  cascade  with  it  phase  shift  between  the  blades. 

The  basic  equation  for  oscillatory  sonic  flow  from 
Section  II  is 


i|>   -  Q^x  +  K"^  =  0  . 


(11-^9) 
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The  boundary-value  problem  can  be  formulated  as  follows: 

1.  There  Is  no  disturbance  upstream  of  the  blade 
leading  edge. 

2.  Flow  tangency  conditions  must  be  satisfied  along 
the  blade  and  along  the  solid  wall. 

This  latter  condition  may  be  stated  mathematically, 

\\>    (X,0)  =  -v(X)     at  Y  =  0  (111-48) 

and  ip  (X,d)  =  0        at  Y  =  ±d         (111-49) 

Denoting  the  Laplace  transformation  of  a  function  with 
respect  to  X  by  an  overbar,  i.e., 

f(p)  =  c£fU)  -  /   e"pX  f(X)  dX  , 

0 

one  obtains  for  the  Laplace  transform  of  Eqn.  (11-^9) 

^yy  "  X?'^  =  °  (IH-50) 

where  X2  =  Qp  -  K2  (111-51) 

The  boundary  conditions  transform  into 

iL  =  -v(p)  Y  =  0  (111-52) 

«y 

and  if  =  0  Y  =  ±d         (111-53) 

A  solution  to  (111-50)  satisfying  (111-52)  and  (111-53) 
has  been  given  by  Miles  [Ref.  19], 


<KP,Y)  =  y   v(p)CSCH(Xd)COSH[A(d-|Y|)]sgn(Y)       (III-^) 
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Restricting  oneself  _to  the  flow  between  the  upper  tunnel 
wall  Y  =  d  and  the  blade  one  has  Y  _>  0 .   Hence  the  absolute 
value  sign  may  be  dropped. 

To  invert  (111-5*0,  a  suitable  expansion  of  the  equation 
was  given  by  Miles  as 


?(p,Y)  =  y  v(p)  I      (EXP[-X(2Nd+Y)]  +  EXP{-X[ (2N+2 )d-Y] } ) 

N=0 


(IH-55) 


Recalling  the  Laplace  transform  of 


f(X)  =  X~1/2  e~a/HX  Rea  >  0 


being 


f(p)  =  |  e  pAf(X)  dX  -  /tP ,  Re  p  >  0  ,        (111-56) 

0  /F 


one  has  for 


f(p,Y)  =  i  EXP[-X(2Nd+Y)] 


where    X  =  ifQp-K 


=  f- 


the  inversion    K2     Q  (2Md+Y)2 

EAf  L-pr  a  -  tt  ^ J 

f(X,Y)  = = - .  (111-57) 

/itQX" 


Therefore,  with  the  aid  of  the  convolution  theorem,  one 
obtains  the  following  solution: 


*(X,Y)  =---=,   /  -~^M-   EXP[^(X-S)] 
/ttQ  -1  /x-s       !" 


.  |   {EXP[_Q  C2NW)2-,  +  EXp[_Q  [(2N.2)d+Y]2]|  dS 


(111-58) 
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On  the  blade  surface  where  Y  =  0,  this  result  simplifies  to 
iKx,o)  =  -  7=,  /' 


/    ~=,  EXP[~  (X-S)] 

m  -i    /xir       Q 


N=0 


EXP[-  §5^]  +  EXP [4  <*»??>  V] 


X-S 


X-S 


I 


dS. 
(111-59) 


For  steady  flov;  where  K  =  0,  one  obtains  the  further 
simplif i cat ion, 

f 

X 


iKx,o)  =  -  -~=z,  i 


/ 


v(S) 


1   /x-s' 


dS 


+  2  f  /  viji  EXP[_X  443  dS 

N=l  -1  V^-S 


X-S 


(111-60) 


This  expression  was  previously  given  by  Sandeman  [Ref.  4] 
in  his  analysis  of  choked  wind  tunnel  flov/.   Equation 
(111-58)  represents  the  generalization  to  oscillatory  flow. 

Appendix  B-II  describes  the  translation  of  (111-60)  into 
a  computer  program.   The  solution  .of  this  function  provides 
the  steady  state  check  of  the  collocation  technique  as  is 
shown  in  Tables  III  and  IV. 
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IV.   PRESENTATION  OF  RESULTS 

Table  II  presents  the  interference  coefficients  for 
the  steady  state  case  obtained  for  the  controlling  param- 
eters indicated  in  the  tables  with  the  computer  routine 
developed  in  Appendix  A. 

Table  III  contains  the  value  of  the  perturbation 
potential  along  the  reference  blade  obtained  from  the 
routine  developed  in  Appendix  B-I  and  utilizing  the  inter- 
ference coefficients  of  Table  II.   These  results  may  be 
compared  with  those  in  Table  IV  obtained  from  the  steady 
state  Laplace  transform  solution  for  the  sonic  wind  tunnel 
wall  interference  problem  of  Section  III-E  and  programmed 
in  Appendix  B-II. 

Table  V  demonstrates  the  erroneous  correction  integrals 
obtained  from  Subroutine  COR  of  Appendix  A  as  the  result 
of  computer  numerical  inaccuracies.   A  comparison  of  the 
interference  coefficients  for  two  cases  before  and  after 
the  anomaly  was  detected  is  presented  in  Table  VI. 

A  sample  of  the  unsteady  interference  coefficients  for 
the  case  X  =  0  is  presented  in  Table  VII. 

Table  VIII  compares  a  sample  steady  state  output  with 
the  output  of  the  general  case  X  i-   0,  K  ?   0,  for  very 
small  K. 
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POTENTIAL     SOLUTION    FOR    THE    FOLLOWING    PARAMETERS 

DL  =       0.200  DK  =       CO  RHO  =       1.000 

THE     INTERFERENCE    COEFFICIENTS 


DEL( 00) = 
DEL(01)= 
DEL (02)= 
DEL ( 03)= 
DEL(04)= 
DEL(05)= 
DELI C6) - 
DEL(07)= 
DEL(08)= 
DEL(C9)= 


-2.0453110E  00 
-1.2504177E    00 

3.0618846E-01 
-1. 5100145E-01 

1.0638958E-01 
-6.7890167E-02 
-2.0572830E-02 

3.3797889E-03 

1.3067639E-01 
-9.7604275E-02 


DEL(10)= 
DEL(ll) = 
DEL( 12)= 
DEL (13) = 
DEL( 14) = 
DEL (15)= 
DEL(16) = 
DEL  (  17)  = 
DEL(18)= 
DEL(19)  = 


2.0453110E  00 
1.2504177E    00 

-3.0618846E-01 
1.5100145E-01 

-1.0638958E-01 
6. 7890167E-02 
2.0572830E-02 

-3. 3797889E-03 

-1.3067639E-01 
9.7604275E-02 


VALUE    OF    THE    POTENTIAL    ALONG   THE    BLADE 


STATION 
X=-0.90 
X=-0.  80 
X=-0.70 
X=-0.60 
X=-0.50 
X=-0.40 
X=-0.30 
X=-0.20 
X=-0.10 
X=-0.  CO 
X=  0.10 
X=  0.20 
X=  0.30 
X=  0.40 
X=  C.50 
X=  0.60 
X=  0.70 
X=  0.80 
X=  0.90 
X=    1.00 


POTENTIAL 

POT=    1.1780910E  00 

POT=    2.1757832E  00 

POT=    3.1754J.12E  00 

POT=    4.1755638E  00 

POT=    5. 1756020E  00 

POT=    6.1755772E  00 

POT=    7.1755638E  00 

POT=    8.1755724E  00 

POT=    9.1755791E  00 

POT=    1.0175587E  01 

POT=    1.1175572E  01 

POT=    1.2175570E  01 

POT=     1.3175570E  01 

POT=     1.4175561E  01 

POT=    1 .5175563  E  01 

POT=     1.6175568E  01 

POT=    1.7175552E  01 

POT=    1.8175476E  01 

POT=     1.9175308E  01 

POT=    2.0176529E  01 
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POTENTIAL    SOLUTION    FOR    THE    FOLLOWING    PARAMETERS 

DL=       0.200  DK  =       0.0  RHO=       2.000 

THE  INTERFERENCE  COEFFICIENTS 


DEL(CO) = 
DEL( 01)= 
DEL (02)= 
DEL(03)= 
DEL( 04)= 
D  E  L  (  0  5  )  = 
DEL(06)  = 
DEL(07)= 
DEL(08)= 
DEL( 09) = 


-8.0229187E-01 
-6. 1007094E-01 

1.4253497E-01 
-6.6098273E-02 

2.7520720E-02 
-1.2785759E-02 

5.6872237E-02 
-5. 7169359E-02 

4.0132940E-02 

5. 0579518E-02 


DEL(10)= 
DEL(ll) = 
DEL(  12)  = 
DEL.  (  13)  = 
DEL(14)=-2 
DEL(15)=    1 
DEL (16)=-5 
DEL(  17)  =    5 
DEL (18)=    4 


8.0229187E-01 

6.1007094E-01 

1. 4253497E-01 

6.6098273E-02 

7520720E-02 

2785759F-02 

6872237E-02 

7169359E-02 

0132940E-02 


DEL (19) =-5. 0579 5 18 E-02 


VALUE    OF    THE    POTENTIAL    ALONG    THE    BLADE 


STATION 
X=-0.90 
X=-0.80 
X=-0.70 
X=-0.60 
X=-0.50 
X=-0.40 
X=-0.30 
X=-0.20 
X=-0.10 
X=-0.00 
X=  0.10 
X=  0.20 
X=  0.30 
X=  0.40 
X=  0.50 
X=  C.60 
X=  0.70 
X=  0.80 
X=  0.90 
X=    1.00 


POTENTIAL 
POT=  8.1581414E-01 
POT=  1.3204842E  00 
POT=  1.8202248E  00 
POT=  2.3200750E  00 
POT=  2.8200855E  00 
POT=  3.3201075E  00 
POT=  3.8201075E  00 
POT=  4. 3200874 E  00 
POT=  4.8200750E  00 
POT=  5.3200960E  00 
POT=  5.8200798E  00 
POT=  6.3200665E  00 
POT=  6.8200712  E  00 
POT=  7.3200731E  00 
POT=  7.8200769E  00 
POT=  8.3200636E  00 
POT=  8.8200598E  00 
POT=  9.3201199E  00 
POT=  9.8202019E  00 
POT=     1.0319364E    01 
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POTENTIAL    SOLUTION    FOR    THE    FOLLOWING    PARAMETERS 
DL=       0.200  DK=       0.0  RHO  =       5.000 

THE    INTERFERENCE    COEFFICIENTS 


DEL(00)= 
DEL(01)= 
DEL{ 02)  = 
DEL(03)= 
DEL (04) = 
DEL(05)= 
D  E  L  ( 0  6  }  = 
DEL(  07)  - 
DEL (08)= 
DEL(09)= 


-1. 1538458E-01 

-1.892513  6E-01 
7. 1354136E-03 
2.2944096E-02 

-2.5751479E-02 
2. 8709877E-02 

-6.2098347E-02 
4. 1337617E-02 
5.2630220E-02 

-5.1 547259 E-02 


DEL(  10)  = 
DEL ( 11)= 
DEL (12} = 
DEL( 13) = 
DEL(14)= 
DEL(15) = 
DEL( 16)= 
DEL(17)= 
DEL(18)  = 
DEL( 19)= 


1.  1538458E-01 

1.8925136E-01 
-7.1354136E-03 
-2.2944096E-02 

2.5751479E-02 
-2.8709877E-02 

6.2098347E-02 
-4.1337617E-02 
-5. 2630220E-02 

5. 1547259E-02 


VALUE    OF    THE    POTENTIAL    ALONG    THE    BLADE 


STATION 
X=-0.  90 
X=-0.80 
X=-0.70 
X=-0.60 
X=-0.50 
X=-0.40 
X=-0.30 
X=-0.20 
X=-0.  10 
X=-0.00 
X=  0.10 
X=  0.20 
X=  0.30 
X=  0.40 
X=  0.50 
X=  0.60 
X=  0. 70 
X=  0.80 
X=  0.90 
X=    1.00 


POTENTIAL 
POT=    7.9726553E-01 

POT=    1.1282368E  00 

POT=    1.38 5 5448 E  00 

POT=    1.6115170E  00 

POT=     1.8232851E  00 

POT=    2.0286140E  00 

POT=    2.2310429E  00 

POT=    2.4321480E  00 

POT=    2.6326466E  00 

POT=    2.8328724E  00 

POT=    3.0329742E  00 

POT=    3.2330236E  00 

POT=     3.4330444E  00 

POT=    3.6330500E  00 

POT=    3.8330574E  00 

POT=    4.0330648E  00 

POT=    4.2330627E  00 

POT=    4.4330502E  00 

POT=    4.6330080E  00 

POT=    4.8337231E  00 
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POTENTIAL     SOLUTION    FOP     THE    FOLLOWING    PARAMETERS 
DL=       0.200  DK=       0.0  RHO=    10.000 

THE    INTERFERENCE    COEFFICIENTS 


DEL(  00)  —  1.  5652659E-03 
DEL(01)=-8.5002407E-03 
DEL(02)=~1.4  88  8  73  0E-02 
DEL( 03)=-5.3168200E-03 
DEL(04)=  7.0907287E-03 
DEL( 05) =-6. 0283695E-04 
DEL( 06)=-4.9  88614  5E-03 
DEL(07)=  3. 8378909E-03 
DEL(08)=  1.6985249E-03 
DEL (09) =-2. 18  47  319E-03 


DEL(10J=  1.5652659E-03 
DEL(11)=  8. 5002407E-03 
DEL(12)=  1.4888730L— 02 
DEL(13)=  5.3168200E-03 
OEL( 14)  =-7. 0907 2 8 7E- 03 
DEL(15)=  6.0283695E-04 
DEL(16)=  4.9886145E-03 
DEL( 17)=-3. 8378  909E-03 
DEL (18) =-1.69 85249 E-03 
DEL(19)=    2. 1847319E-03 


VALUE    OF    THE    POTENTIAL    ALONG    THE    BLADE 


STATION 
X=-0.90 
X  =  -0.80 
X=-0.70 
X=-0.60 
X=-0.50 
X=-0.40 
X=-0.  30 
X=-0.20 
X=-0.10 
X=-0.00 
X  =  0.10 
X=  0.20 
X  =  0.30 
X=  0.40 
X=  0.50 
X=  0.60 
X=  0.70 
X=  C.80 
X  =  0.90 
X=    1.00 


POTENTIAL 
POT=    7  .9786545E-01 

POT  =    1.1283703E  00 

POT=     1.3819685E  00 

POT  =    1.5957613E  00 

POT=     1.7841234E  00 

POT=    1.9544506E  00 

POT  =    2.11U937E  00 

POT=    2.2573242E  00 

POT=    2.3949909E  00 

POT  =    2.5258138E  00 

POT=    2.6510839E  00 

POT=    2.7718019E  00 

POT=    2.8888006E  00 

POT=    3.0027514E  00 

POT=    3.1142006E  00 

POT=    3.2235994E  00 

POT=    3  .3313122E  00 

POT  =    3.4376440E  00 

POT=    3.5428400E  00 

POT=    3.6471376E  00 


TABLE    III 


60 


UJ 

z> 

a 

t— 4 

■z. 

X 

o 

UJ 

\- 

LU 

o 

z 

LU 

a: 

UJ 

Li- 

o 

ar 

o 

UJ 

o 

t- 

t 

_J 

o 

< 

X 

3 
1 

a: 

1 
UJ 

o 

in 

< 

a: 

O 

o 

o 

o 

o 

o 

o 

o 

o 

rH 

I— t 

rH 

■— i 

i—i 

r-l 

r-l 

1— 1 

r-l 

1— 1 

rH 

_l 

UJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

CO 

h- 

UJ 

U-l 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

LU 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

s: 

Is- 

r- 

o 

in 

in 

>*• 

m 

<f- 

m 

o 

r- 

n0 

o 

o 

vO 

vO 

O 

vO 

-o 

vO 

X 

< 

vO 

o 

vO 

^0 

vO 

o 

*0 

vO 

-o 

o 

O 

O 

o 

o 

o 

o 

vO 

o 

o 

vO 

r- 

a: 

o 

o 

o 

s0 

o 

o 

o 

o 

o 

o 

o 

vO 

o 

o 

NO 

o 

o 

o 

O 

vO 

< 

_J 

vO 

O 

o 

vO 

^> 

vO 

n0 

o 

vO 

r— 4 

rH 

rH 

rH 

rH 

rH 

rH 

r-l 

I— 1 

rH 

rH 

o 

O- 

o 

< 

rH 

rH 

T-* 

.—i 

rH 

i—4 

f— ( 

1— 1 

r~l 

o 

rH 

CM 

m 

<f 

in 

n0 

h- 

CO 

0> 

o 

z 

• 

1—4 

o 

o 

o 

1- 

rH 

CNJ 

ro 

<f 

in 

o 

r- 

co 

O 

.— 1 

rH 

rH 

f-$ 

r-l 

r-l 

rH 

rH 

r-l 

r-H 

CM 

_J 

2 

z 

< 

►H 

3 

ii 

UJ 

_l 

O 

^ 

o 

< 

_! 

o 

CL 

II 

II 

n 

ii 

n 

II 

ii 

11 

II 

II 

II 

II 

n 

II 

II 

II 

II 

II 

II 

II 

►—1 

_J 

h- 

1- 

h- 

h- 

t- 

h- 

h- 

h- 

t— 

H 

H 

r- 

H- 

h- 

r— 

1— 

h- 

H 

h- 

I- 

h- 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

a 

o 

a 

O 

O 

O 

2 

u. 

a. 

a. 

Q. 

a. 

CL 

O- 

a. 

a. 

o_ 

o_ 

a. 

a. 

a. 

CL 

a. 

CL 

Q. 

CL 

O. 

G 

UJ 

H 

UJ 

O 

X 

a, 

H 

UJ 

cc 

X 

o 

\- 

u_ 

o 
o 

LL 

z 

C\J 

O 

Q 

• 

z 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

O 

1—4 

o 

a 

o 

CO 

r- 

o 

in 

>3- 

ro 

C\J 

rH 

o 

^H 

c\j 

m 

-4- 

in 

o 

Is- 

CO 

c* 

o 

UJ 

t— 

1—4 

3 

r> 

H 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

O 

o 

rH 

-J 

_i 

II 

< 

i 

1 

l 

1 

1 

1 

I 

I 

1 

1 

< 

o 

_J 

1- 

n 

II 

ii 

II 

ii 

II 

ii 

ii 

II 

II 

II 

II 

ii 

II 

ii 

II 

II 

II 

II 

II 

> 

OO 

o 

oo 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

> 

H 
W 

rH 

< 


61 


Ill 

Z> 

o 

t— 1 

2: 

X 

0 

LU 

H 

LU 

O 

"Z. 

LU 

c£ 

LU 

LL 

O 

CC 

O 

LU 

O 

H 
•J 

< 
Eh 


cm 


_J 

11 

_J 

0 

< 

X 

2 
1 

C£ 

1 

LU 

O 

00 

< 

c£ 

H 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.-< 

_J 

LU 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CO 

f~ 

1 

LU 

IU 

UJ 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

UJ 

LU 

UJ 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

5: 

O 

CM 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CM 

CO 

CM 

1— 1 

■—1 

CM 

r~ 1 

r— 1 

1— 1 

CM 

CO 

X 

< 

r- 

CO 

CO 

CO 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

K- 

cc 

co 

CO 

CO 

CO 

co 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

< 

_j 

t— 1 

CO 

01 

CO 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

0") 

CO 

O 

Q_ 

0 

< 

rO 

CO 

CO 

CO 

co 

CO 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

0 

z: 

• 

n 

0 

O 

0 

k 

CO 

r-l 

r-l 

C\! 

C\J 

CO 

CO 

4- 

>t 

LT\ 

in 

0 

0 

r-~ 

r- 

CO 

co 

c* 

C* 

1— 1 

_j 

z 

z 

< 

3: 

II 

LU 

_i 

O 

^ 

O 

< 

_J 

0 

Q. 

it 

II 

II 

II 

II 

II 

II 

II 

II 

II 

It 

II 

II 

II 

II 

II 

II 

II 

II 

II 

M 

_J 

h- 

(~ 

h- 

t~ 

h- 

h- 

h- 

h- 

\~ 

t- 

1— 

h- 

h- 

h- 

h- 

t— 

h- 

K 

h" 

h- 

1— 

O 

O 

O 

0 

0 

0 

CJ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 

0 

0 

0 

Z 

LU 

O- 

a. 

D_ 

o_ 

o_ 

a. 

CL 

fX 

o_ 

o_ 

a. 

Q. 

o_ 

Ou 

0, 

CL 

CL 

Cl 

a 

a. 

LU 

(~ 

LU 

O 

X 

Q. 

\- 

LU 

CC 

X 

0 

h- 

LL 

0 
0 

UL 

z 

00 

O 

0 

• 

Z 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I — 1 

0 

O 

c> 

CO 

h" 

0 

in 

•4- 

CO 

CM 

1—1 

0 

r-l 

CM 

CO 

vt- 

LT\ 

nO 

r- 

CO 

c> 

0 

LU 

1- 

1 — 1 

Z) 

Z> 

\- 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I— 1 

_) 

_J 

II 

< 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

< 

a 

_l 

t- 

II 

11 

II 

II 

It 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

It 

II 

II 

> 

(/1 

a 

00 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

62 


m 

<3 

X 

o 

LU 


UJ 

o 


LU 
CC 
LU 

Li- 
ar 

LU 


o 
o 
o 


in 


_J 

II 

_J 

O 

< 

X 

3 
1 

a: 

1 
LU 

a 

(/) 

< 

a: 

r-l 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

_j 

LU 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

co 

H 

1 

LU 

LLJ 

LU 

IU 

LU 

LU 

ILI 

LU 

LU 

LU 

LU 

UJ 

UJ 

LU 

UJ 

LU 

LU 

UJ 

UJ 

UJ 

LU 

LU 

s: 

in 

h- 

o 

o 

o 

o 

CM 

CM 

rH 

st" 

tn 

O 

rH 

CxJ 

CO 

CM 

ro 

O) 

CO 

ro 

X 

< 

00 

O 

o 

CO 

in 

CO 

m 

sj- 

o 

H 

CM 

CM 

CO 

CO 

CO 

CO 

ro 

ro 

CO 

ro 

t- 

a: 

oo 

co 

in 

rH 

ro 

CO 

rH 

CM 

Cvl 

Ol 

CO 

CO 

CO 

ro 

CO 

ro 

ro 

01 

ro 

ro 

< 

_J 

r- 

CM 

CO 

r-| 

(M 

CM 

ro 

CO 

CO 

ro 

01 

CO 

CO 

CO 

CO 

ro 

CO 

ro 

0J 

ro 

o 

O- 

o 

< 

o 

rH 

01 

vO 

CO 

o 

CM 

<t" 

O 

co 

o 

CM 

vt 

o 

00 

o 

CM 

<r 

s0 

oo 

■z. 

• 

►— < 

o 

o 

o 

1- 

r» 

I— 1 

i— i 

rH 

rH 

CM 

CM 

CM 

CM 

CM 

CO 

CO 

ro 

01 

CO 

>t 

si- 

NT' 

vt- 

si- 

_J 

^L. 

Z 

_ 

< 

2 

II 

LU 

_l 

O 

*L 

O 

< 

_J 

d 

a. 

ii 

II 

ii 

11 

II 

II 

II 

11 

II 

II 

ii 

II 

n 

II 

II 

II 

ll 

II 

II 

ll 

r-H 

_) 

h- 

Y~ 

l- 

h- 

h- 

1- 

h- 

H 

h- 

H- 

H 

h- 

i- 

h- 

f- 

h- 

J— 

I— 

t- 

h- 

\- 

o 

O 

O 

o 

O 

o 

O 

O 

O 

O 

O 

O 

O 

o 

o 

o 

o 

a 

o 

o 

o 

z: 

LL 

a. 

o_ 

a. 

o_ 

o_ 

Q. 

a. 

0l 

Q. 

o_ 

a. 

O- 

O- 

Q. 

o_ 

O- 

a. 

a. 

CL 

Q. 

LU 

H 

LU 

O 

X 

Q. 

H 

LU 

a: 

X 

o 

h- 

LL 

o 
o 

LL 

z 

CNJ 

O 

o 

• 

z 

o 

o 

o 

o 

O 

o 

Q 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

t— < 

o 

a 

a> 

CO 

r- 

vO 

in 

<f- 

CO 

CM 

.H 

o 

rH 

CM 

co 

«fr 

LT» 

o 

r- 

CO 

c* 

o 

LU 

h- 

l—t 

r> 

Z> 

h- 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

rH 

-J 

_J 

II 

< 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

< 

O 

_1 

f- 

ii 

II 

ii 

ti 

II 

II 

ii 

II 

II 

II 

11 

II 

II 

II 

II 

II 

II 

II 

II 

II 

> 

I/) 

Q 

00 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

> 

M 

W 

-»• 

J 

CQ 

< 

B 

63 


3 
C? 
»— < 

z 

X 

o 


LLJ 

o 

z 

UJ 

a: 

UJ 

LL 

or 

UJ 

\- 


o 
o 
o 

c 

o 


_J 

II 

_J 

O 

< 

X 

3: 
1 

a: 

1 

UJ 

0 

00 

< 

QL 

rH 

0 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

0 

_J 

UJ 

0 

0 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

0 

co 

►- 

1 

UJ 

UJ 

LU 

Ul 

Ul 

LU 

Ul 

UJ 

Ul 

UJ 

UJ 

LU 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ul 

UJ 

21 

in 

CO 

CO 

Is- 

ro 

0 

0 

m 

O 

CM 

O 

rH 

O 

in 

O 

0 

CM 

in 

in 

r-l 

X 

< 

00 

CO 

a> 

Is- 

rH 

<r 

<\l 

m 

O 

CO 

O 

CO 

co 

r- 

CM 

0 

ro 

0 

CO 

rH 

H 

OL 

CO 

CO 

rH 

m 

■+ 

<r 

rH 

Is- 

in 

in 

t— 1 

rH 

CO 

CM 

<*- 

m 

r-t 

Is- 

CM 

r- 

< 

_J 

Is- 

cm 

CO 

O 

CO 

m 

rH 

in 

cr> 

CM 

m 

Is- 

CO 

0 

rH 

CM 

ro 

m 

-^ 

<t- 

O 

Q. 

0 

< 

O 

1—1 

m 

lt\ 

Is- 

0 

i-l 

CM 

m 

in 

0 

Is- 

CO 

0 

rH 

CM 

ro 

-t 

in 

0 

Z 

• 

►—4 

a 

O 

0 

h- 

I'- 

r-l 

<~\ 

rH 

I— { 

rH 

CM 

CM 

cm 

CM 

CM 

CM 

CM 

ro 

ro 

ro 

CO 

CO 

CO 

CO 

_j 

Z 

Z 

< 

►— 1 

2 

11 

UJ 

_j 

O 

*: 

O 

< 

_J 

Q 

O- 

ll 

II 

11 

II 

II 

II 

II 

11 

11 

11 

11 

II 

II 

n 

11 

11 

11 

11 

11 

11 

►— 1 

_J 

h- 

1- 

h- 

K- 

H 

H 

1— 

»- 

h~ 

1- 

i— 

H 

h- 

K- 

h- 

i~ 

1- 

H- 

l~ 

h- 

h- 

O 

O 

O 

O 

O 

O 

O 

O 

0 

O 

0 

O 

O 

O 

0 

O 

0 

0 

O 

O 

0 

z 

LL 

a. 

CL 

CL 

Cl 

CL 

CL 

o_ 

a. 

CL 

0. 

CL 

O. 

CL 

CL 

CL 

CL 

CL 

CL 

CL 

CL 

UJ 

h- 

UJ 

O 

X 

Q. 

h- 

UJ 

ct 

X 

O 

1- 

u_ 

O 
O 

U. 

z 

cm 

0 

0 

« 

z 

a 

O 

O 

O 

O 

O 

0 

0 

O 

0 

0 

O. 

O 

O 

O 

O 

O 

O 

O 

O 

t— 1 

O 

0 

0s- 

co 

Is- 

vO 

in 

sf 

CO 

C\J 

rH 

0 

1—1 

CM 

CO 

<" 

in 

O 

Is- 

CO 

O^ 

O 

UJ 

h- 

1—4 

3 

3 

H- 

0 

0 

O 

O 

0 

O 

0 

0 

O 

0 

0 

O 

0 

O 

0 

O 

O 

0 

O 

rH 

_J 

_l 

11 

< 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

< 

O 

_J 

t— 

11 

II 

II 

II 

11 

II 

II 

11 

II 

11 

II 

II 

II 

II 

11 

II 

II 

II 

II 

II 

> 

l/) 

Q 

00 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

> 

_ 

M 

H 

J 

CO 

< 

Eh 

64 


CM 

O 

o 
o 
o 
in 

CM 


< 
> 

< 


cc 
o 


oo 
O 


o 
o 

in 


ii 

o 

X 


OO 

cc 

UJ 


< 

a: 

< 

CL 

o 

o    o 
2    in 


o 

_i  ii 

_i  *: 

o  o 


UJ 

x 


ii 
or 

< 


o    _i 

OO      D 


< 

t-l 

< 

_J 

CtL 

h- 

O 

»— 1 

z: 

UJ 

s: 

CD 

h- 

>—i 

7Z. 

_j 

X 

t—t 

a: 

o 

•>. 

00 

UJ 

o 

UJ 

O 

CL 

w(\j 

o 

O- 

o 

00O 

ID 

in 

O  I 

a: 

CO 

o 

o 

CD 

rr^in 

U- 

• 

OCM 

t 

UJ 

o 

ZO 

ID 

t—i 

_J 

00  || 

< 

II 

—  UJ 

> 

o 

in 

_l 

o 

CMO 

_J 

< 

1  u. 

< 

a: 

-> 

o£ 

O 

wO 

o 

UJ 

*z 

UJ 

1— 

*•< 

h- 

Z 

X 

z 

o 

^-^ 

v-ro 

K— 1 

UJ 

JO 

oo 

_J 

<  1 

^ 

3 

UJ 

aro 

►— « 

Z 

o 

oo 

00 

2 

00 

o 

UJO 

O 

UJ 

t-O 

•— < 

cd 

o 

2    • 

1- 

U- 

CM 

t-ntn 

o 

CM 

UJ 

00 

CO 

a  ii 

&. 

o 

CO 

t-_j 

cc 

o 

co 

< 

O 

• 

z 

o 

o 

o 

OO 

•> 

< 

h- 

o 

II 

DO 

n 

z 

00 

_J  a: 

_i 

»— t 

oo 

OUJ 

~3 

o 

00 

00 

OOM 

■H       ro       <-0       CO 


o 


o 


cMmNO^om<f"in<f- 

i— I        ,_|        ,-1        r- »^Hi— »        ,— )        i~» 


1 

o 

1 

o 

1 

o 

1 

O 

1 

o 

1 

o 

1 

o 

o 

> 

o 

1 

o 

1 
O 

1 

o 

I 
o 

in 

(M 

<r 

r- 

0> 

in 

in 

CO 

vf 

CO 

a> 

r-l 

a< 

CO 

rH 

CO 

co 

o 

<~i 

in 

o 

rH 

0> 

CM 

r- 

>t 

00 

r-< 

in 

o 

CO 

m 

■+ 

CO 

vO 

in 

CO 

<f 

•*■ 

•-I 

O 

o 

co 

in 

<* 

m 

r-( 

CM 

r- 

o 

r- 

CO 

1— t 

>-* 

o 

•4- 

CM 

t— i 

r» 

r-4 

CO 

1 

CM 

CM 

co 
1 

CM 
1 

comf^-ocM<j-Oom*tinvf 

OOOrHi-lrH.— (r- It-lrHf— I.— I 


o 

1 

o 

1 

a 

1 

1 

o 

1 

o 

1 

o 

1 

O 

1 

o 

1 

o 

1 
o 

1 

o 

1 

o 

CO 

in 

o 

st- 

sf 

h- 

o 

«-H 

CO 

CO 

co 

«*■ 

c* 

o 

r- 

o 

O1 

o 

CO 

CO 

in 

o 

f— 1 

CO 

co 

r- 

CM 

r- 

in 

vO 

1—1 

o 

CO 

Is- 

CO 

rO 

rO 

in 

CM 

i—i 

■4- 

o 

CM 

o 

CO 

CO 

a> 

0^ 

in 

vO 

1~ 1 

^' 

~j 

in 

CM 

•—I 

vO 

CM 

r-i 

r-4 

1— « 

»*■ 

1 

1 

|-H 

cm     co     <t     m     >o 


co     c> 


cm     co 


65 


o 

Q 
O 

o 

in 

CM 


g 


w 

►J 

CQ 

O 

< 

O 

Eh 

O 

• 

in 

ii 

O 

X 

OL 

oo 

DC 

lil 

r- 

LU 

51 

< 

cC 

< 

Q. 

O 

o 

O 

z 

in 

*— 1 

t 

3: 

o 

O 

_J 

(1 

_J 

^ 

o 

o 

LL 

UJ 

X 

H 

a: 

o 

u_ 

00 

■z. 

o 

o 

r-< 

• 

H 

o 

3 

_J 

II 

O 

_J 

oo 

o 

II 

•       < 

t 

a: 


00 


< 

►—I 

< 

t- 

0 

t—4 

z 

UJ 

s: 

o 

K 

t-H 

■z. 

_J 

X 

►—i 

f£ 

o 

^ 

00 

UJ 

o 

UJ 

O 

a. 

<— o 

O 

0. 

o 

000 

Z> 

o 

or 

r- 

o 

o 

r- 

aim 

u_ 

• 

Ocm 

• 

LU 

o 

zo 

Z) 

*—l 

_J 

00  || 

< 

ii 

—,UJ 

> 

in 
•  a: 

_J 

o 

CMO 

~) 

< 

1  u. 

< 

a' 

~3 

or 

O 

~o 

0 

UJ 

#2: 

UJ 

(~ 

£•  < 

h- 

2: 

X 

z 

Q 

t-H 

»— CM 

1— I 

LU 

JO 

00 

-J 

<t  1 

^ 

3 

UJ 

0:0 

►—I 

z 

o 

00 

00 

2 

CO 

o 

LUO 

O 

UJ 

>~o 

r- < 

aL 

o 

z:  • 

H 

U- 

CO 

t-tin 

O 

O 

UJ 

00 

r~ 1 

O  II 

a: 

O 

r» 

1 1 

a: 

o 

4- 

< 

a 

• 

2: 

o 

o 

o 

OO 

r 

< 

1— 

o 

II 

z>o 

ii 

2: 

00 

_ia: 

_i 

»— < 

00 

Ouj 

-5 

a 

oo 

00 

oOM 

cMromvomir\inm(M 


in     co 


0 

1 

O 
1 

O 

1 

0 
1 

0 

1 

0 

O 

O 

O 
1 

O 

O 

O 

O 

1 

0 

1 

O 

1 

O 

1 
0 

1 

Q 

1 

0 

1 

O 

1 
O 

1 

0 

O 

O 

O 

a 

CM 

CO 

CM 

r- 

r- 

r-l 

in 

r-l 

r- 

CO 

^ 

CM 

r-l 

m 

0 

0 

CM 

r» 

in 

CO 

CO 

CO 

en 

O 

in 

CO 

r-i 

r-l 

<}- 

CO 

r- 

vt 

cm 

sf 

r- 

ro 

r-l 

in 

<■€> 

CM 

r-l 

0 

r-l 

r~ 

r- 

co 

r-l 

st- 

vf 

CM 

O 

O 

CM 

r-l 

vO 

1 

1— 1 
r 

r-l 

CM 

^0 
1 

r-l 

1 

r-l 
1 

r-l 

r-l 

rH 

1 

CM 

ro 

in 

m 

in 

m 

in 

in 

CM 

r-l 

t^i 

in 

r- 

O 
1 

0 

t 

O 

0 

0 
1 

0 
1 

0 

0 
1 

0 
1 

O 

0 

0 

0 

1 

a 

1 

0 

1 

0 

I 

1 

1 

0 

1 

O 

1 

0 

1 

0 

O 

0 

a 

0 

m 

0 

CO 

CO 

0 

CM 

en 

in 

CO 

<f 

CO 

OD 

r- 

CJ> 

<f- 

CM 

r- 

CO 

CO 

CM 

in 

vO 

•4" 

<!- 

0 

0 

0 

a 

r- 1 

-4- 

p) 

CO 

CO 

O 

CO 

in 

CM 

■co 

in 

CO 

co 

CM 

m 

CM 

^ 

0^ 

r- 

O 

0 

in 

h- 

•4- 

CM 

r-l 

r- 

r-l 

1 

CM 

1 

1—1 

00 

l-l 

1 

C\J 

CO 

CNJ 
1 

» 

r-i     cm    co     -4-     in     v0 


CO     c^ 


r-i    cm     ro 


66 


00<— If-fr-lr-HCNJCMCM^Hi— 1 

o 

M 
Eh 

O 

w 

Eh 

0<— lrHi-»r- » CN(  CXI  rsj ,— | , — | 

oooooooooo 

IAOOOOOOOoOO 

W 

1  1  1  1  1  1  1  1  1 

O        1     1     1     1     1     1     1     1     1 

Q 

O  U  J  UJ  UJ  UJ  UJ  UJ  UJ  UU  U J  UJ 

OUJUJUJUJUJUJUJUJUJUI 

•  CNJ  O  <~<  i>  h-  m  LT>  f-  O  r-l 

•  cm  o  i-t  o^  r-  in  Is-  <j-  o  cm 

« 

or-  i-t  r*-  irw  m  --t  .-h  m  o 

or-i-tr—  LOsom—trHLno 

O 

o-tmr-tMcsjoc^rH^- 

0<(-inN  c\jrvjO0M-H-4- 

PC 

II    O  Is-  <t"  00  C\|f-.r-<vO  |— r-t 

ii  or-Ntooc-jr^-i-tsor*-'-* 

_J«««ef»*«*e 

<tr-ir—  Ni"»-ti-lv0rOr-<rHi-< 
II              1              III              1 

w 

W 
LH 

<-— !r— ^t-4r- tvOC^-— tf— tf-t 

II      1      III      1 

•-)>-<c\ic\|fMf\jc\jrM(\jcNJ 

f-»f-tCNJ(\ICNJC\|CNJc\jr\JCNJ 

oooooooooo 
i   i   i    i   i   i   i   i   i   i 

OOOOOOOOOO 

1      1      1      1      1      1      1      1      1      1 

< 

1/) 

i   i   i    i   i   i   i   i   i   i 

UJ  UJ  UJ  UJ  Ul  UJ  U  J 1 U  UJ  UJ 

1     1     1     1     1     1     1     1     1     1 

UJUJUJUJUJUJUJUJUJUJ 

H 

r-i  r-  r-f\ji"-r-^-ococ\! 

r-4r-r--rvji>.r~-vtococ\i 

Q 

Z 

oc\JocjNO,-<c\ir0Lr\in 

O  C\J  O  CJN  O  H  t\J  rO  LTV  LH 

S 

UJ 

Ol^vOcOlAcO^OOC^ 

oiAocoinroo^o^oo 

< 

i— i 

vt^ooooomr^-^-sot^- 

^rococ^ooror-vtvOr- 

<_> 

o********** 

O     »<■«•»<:••«» 

W 

»— i 

OsOCr«sO(NJrOroLnoo^O^ 

OvOONvOCNJrnmLnoa>o 

cr; 

u_ 
u_ 

UJ 

m                  til 

in                  III 

o 

w 

O  II    II    II    II    II    11    II    II    II    II 

O  ||    II    II    ||    II    II    ||    II    ||    II 

o 

*-^  •-**"•%  -~^  "^^  .*-"•  *-*.  «-"*  *-*  *—» 

«-S--*^-»--»»— *^«»  *-«»-*■*  <*^*-* 

CQ 

o 

II  o>-ic\iro^"ir>Oi>-coo^ 

ii  Of~tc\jro^Lnvor—coo 

2 

Of-<t-<r-<l— t^Hr-(r-Ht-(r-4'-t 

Oi — ••— tr—Ji— «i— 1«— tr— It— If— Ir— t 

a:^_J_i_j_J_i_i_i_j_J 

00 
Eh 

o 

a:  _J  _j  _j  — i  _j  _i  _i  _i  _i  _j 

f— 1 

UJ  UJ  UJ  U  J  UJ  UJ  UJ  UJ  UJ  UJ 

UJUJUJUJUJUJUJUJUJUJ 

w 

H 

oooooooooo 

oooooooooo 

M 

u 

O 

■z. 

H 

Z> 

Ph 

u_ 

Of— Ir—lt—lr-^CNJCMCNjr-f'H 

Of-tl— lf-tr-lC\JC\JC\Jf-tf-l 

LU 

OOOOOOOOOO 

OOOOOOOOOO 

o 

o 

1  1  1  1  1  1  1  1  1 

1  1  1  1  1  1  1  1  1 

o 

z 

UJ  UJ  UJ  UJ  UJ  UJ  U  J  iJJ  LU  UJ 

UJUJUJUJUJUJUJUJUJUJ 

UJ 

ocxjOr-iCM^-roiAr— o— «• 

ocNjoM-tor—Lnr-.^-ocNj 

w 

CC 

o  r—  ■-<  f~-  lo  no  m  i-t  .-i  m  o 

o  f-  i~i  Is-  in  o  ro  f-i  »-i  m  o 

O 

UJ 

mo<j-Lnr-r\JC\iOc^t-,»t 

LnOvj-inr-cMr\joc7>i-tvj- 

s 

U_ 

or 

UJ 

•  oNvh~«^'Co<\irv-'~'NO^,-H 

•  CM fre0C\jr*-r4<Oh-i-l 

w 

i-t  h-  «^  rH  r-(  s0  CO  —1  l-H  lH 

HN>t'-tHvOr<1H'-tr-l 

f— 1 

H           1        1               1 

11        1      1           1 

Of-<r-IC\|<\!CNjC\|r\Jc\|CNJC\) 

Or-<-~icNJ<Nr\jc\j<\ic\jc\jc\j 

oooooooooo 

1      1      1      1      1      1      1      1      1      1 

OOOOOOOOOO 

1      1      1      1      1      1      1      1      1      1 

Eh 
M 

1     1     1     1     1     1     1     1     1     1 

UJ  UJ  UJ  UJ  UJ  LU  UJ  UJ  UJ  UJ 

1     1     1     1     1     1     1     1     1     1 

UJUJUJUJUJUJUJUJUJUJ 

•Hr—r--  ojr-r-^ooocNj 

rHr-r-cMr-r-Ntc^cocN 

O  cnj  o  C*  o  rH  Cvi  ro  tn  in 

O  cm  o  o>  o  .-I  cnj  ro  m  in 

oirivocoLncoo^OvOcr« 

OLnvOcoLnoooOvoo 

*tr»icooNcoror-«tOf'- 

vj-roooo>oor0f->tor— 

< 

•  ••t*«««ti 

H 

vOOOC\JrorOu,\oO>0> 
1    1    1    1        II           1 

vOfjNvOcNjmroinvOCT'O^ 

I  1    1    1        II            1 

II  n  ii  ii  ii  ii  il  ii  ii  ii 

> 

W 
•J 

< 

M  n  il  H  II  ii  li  n  n  ii 

•o^-tcNirnvt^vor^coa* 

•0»-trNJro>tLn>Ois-cofJ> 

ooooooooooo 

ooooooooooo 

Eh 

— lUJUJUJUJUJUJUJUJujUJ 

ooooooooooo 


_JUJUjUJUJUJUJUJUJUJUJ 

ooooooooooo 


67 


CO 

Z 

UJ 
t—i 

o 


UJ 

o 
o 

■z. 
o 


3 
LL 

O 

UJ 


uj 


ro  r-i  xi-  co  ro  c\j  ro  ro  co  m 

oooooooooo 

I    I    I    I   I    I    I    I    I    I 

OUJUJujLUUJUjUJlULUUJ 

•  O  ro  cOco  o  -4-  O  •-(  cm  co 

OcomcoorM*4-OrHincM 

>t  vj-  <t-  <j-  i*»  t*-  <J-  r*~  f*-  h- 

ii  roovO.~<oc\jc^om^ 

_j...«...... 

<  CM  in  <?"  0'*  C\J  fl  CO  sO  r-lr-4 

1  L      » '  L 

•-HrorOcocNKOcMcOcOvO 

oooooooooo 

I  I  I  I  I  I  I  I  I  I 

UJ  LULU  LUIJJ  UJ  UJ  IJJ  U  1  LU 

-4-^moc\iOrMr~-c\ir«- 
f>-corocomcoor*-mvO 
coco^Qr\|incocomcMvo 

coro^rococxj^-NtfMvO 

0»t«»»e««»« 

OOCT'CM^H^vtr-lCNJsOrO 
O  I  ||  I 

in  ii  II  ii  n  II  ii  ii  ii  ii  il 

ii  o  r-t  cm  ro  ^"  m  vo  r*-  co  t> 

0>— «. — 1< — « I— «  •— I  T— «  !--(>— lr-4,-4 

cjc:  .j  — i  _j  — i  _i  „  j  _i  _i  _j  — i 

LU  LU  UJ  LU  LU  LU  IU  IJJ  IU  LU 

DODQDDOQQQ 


ro  rH  -4"  ro  ro  <m  o'  j  ro  ro  m 

oooooooooo 

I  I  I  I  I  I  I  I  I  I 

U  J  LU  UJ  LU  IU  LU  U  J  U  J  LU  U J 

Oo  roro  coo  »*- o  r-icMoo 
OcomcoOnjvj-vO'HmrM 

•  COOvO'--iC7>rMOOm<J- 
O    t .    «    • 

r\iLO<f-cr'CMr-irOv0r-i.-< 
M  III  III 

it  ►•.•,•.»►►•.»*  •> 
O^rOrocoOsJrocMrOrovO 

OOOOOOOOOO 

I  I  I  I  I  I  I  I  I  I 

LU  LU  LU  LU  LU  tU  UJ  LU  Ul  LU 

**r-<rmCr>iMO'-ir*-c\ir- 
h-  ro  ro  co  in  co  or-  in  vO 
co  co  vo  cMin  co  co  m  cm  sO 
co  ro  in  ro  ro  cm  ■4*  «J*  cm  ^O 

0>>O(Mr-«»-i<f-'-<njv0r0 
III  II        I 

ii  ii  ii  n  ii  ii  ii  ii  ii  ii 
O-*-"- — ■ — >~ 

«o«-ir\jcOvtmvor*-coo*> 
ooooooooooo 


_J  LU  Ui  UJ  IU  LU  LU  UJ  UJ  LU  UJ 

ooooooooooo 


OhhoOOHOhoh 

ooooooooooo 
m  |  | 

O  LUUJLU  LU  LU  UJ  UJl  U  UJUJ 

•  o^"coinr*-m^trocop*~ 

OcgOrv,  v^-r-cMOOcooo 

cosOr-CMOr-mOcoLn 

ii  mrooNcovOvoromo<t- 

_j   ••••••••«• 

<t  H  in  t— •  in  CO  <— J  i— <  CM  r-H  i—l 

II        II  I 

C\|HoOOhOhoh 
OOOOOOOOOO 

I  I 

UJ  UJ  LU  LU  UJ  LU  LU  LU  LU  LU 

rjNO^j-or—coChvDcor— 
cm  vj-in  cm  O  <f-  ro  co  Is-  ■-« 
roc\ico-4T*-o^comoin 
co  co  cm  ro  ro  o -o  O^  ro  m 
o********** 

Oi-^OCM^Of*-— (OCMr-lr-l 

O   I  ||  II 

in  ll  ii  ii  il  ii  ii  ii  ii  ll  ti 
II  Or-4cMCOvtmvOr-coo 

Or-*-— lr-Jr-l,-<l— 1<— 1'-4>— 1-— I 

X  —  —  «—  w  <—  —  w«-  w  w 

Qd  — » I  — »_J I  — t I  — t  — » I 

UJUJIUUJUJLULUIUUJUJ 

oooooooooo 


H^OOOHOrlOrl 

oooooooooo 
I  I 

LU  UJLU  LU  UJ  UJ  UJ  LU  U  J  LU 

o  O  ^t  co  m  r—  in  ^  ro  co  r- 
o  cMOr— ^  r-O  roo  coco 
mcOvor—cMor—  mocoin 

•  mcoOCOsOOcOinOv}* 
o«*«******* 

i— ( in  ~-i  m  ro  i— i  f-*  cm  r~t  r-i 
ii  II  III 

ONt-iOOO^O^OM 
OOOOOOOOOO 

I  I 

UJ  UJUJ  UJ  LU  IU  UJ  UJ  UJ  UJ 

C' o  <4  o  r- ro  o  vo  co  r- 

CM^inC\!CT>*J"COCOr*-r-l 

rocMcOvtr*-ocomOLn 
cocor\jcorocj>vOorom 

rH  O  CM  O  h-  r-i  vO  cm  <-t  — « 

II  II  I 

II  H  II   il  II  H  ll  II  il  II 
o -.^.—-^-»— .~— — — » 

•  o  i-(c\j  co  ^  in  <j  r— co  o 
ooooooooooo 


—ILULULUUJUJUJUJUJUJUJ 

ooooooooooo 


£3 

o 

M 
Eh 
O 
W 

W 
Q 

o 

K 
f£ 
W 

w 

Eh 
< 
Q 
<S 
W 
O 

w 

CO 
Eh 

W 
M 
O 
H 

H 
O 

o 

w 
o 

w 

K 
W 

Ph 

K 
W 
Eh 
12 
M 


CO 
M 
> 

W 

CQ 
< 
Eh 


68 


o 

o 

O 

rH 

r-H 

rH 

rH 

CM 

CM 

CM 

rH 

r-i 

o 

rH 

rH 

rH 

CM 

rH 

>-^ 

rH 

rH 

rH 

rH 

in 

a 

o 

O 

o 

O 

O 

o 

O 

a 

o 

o 

O 

o 

O 

O 

O 

o 

o 

O 

O 

o 

o 

1 

1 

1 

1 

1 

i 

1 

1 

1 

rH 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

o 

UJ 

LU 

UJ 

UJ 

UJ 

Ul 

UJ 

UJ 

UJ 

UJ 

o 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

UJ 

Ul 

UJ 

• 

<\l 

o 

rH 

0* 

r- 

m 

r- 

«4" 

o 

CM 

• 

o 

en 

in 

CM 

en 

Is- 

r- 

CO 

vO 

*± 

o 

r- 

rH 

Is- 

in 

vO 

en 

H 

rH 

in 

o 

o 

fj> 

o 

<f 

O 

CM 

o 

CM 

o 

in 

o 

o 

-4" 

m 

r- 

CM 

CM 

o 

O 

r-i 

"4- 

en 

in 

in 

f- 

o 

CM 

O 

rH 

r- 

0> 

II 

O 

r- 

«t 

co 

CM 

Is- 

r-< 

v0 

Is- 

rH 

II 

o 

rA 

CM 

O 

<i- 

CM 

O 

CM 

<}■ 

m 

_J 

_l 

< 

l 

r- 
I 

<*■ 

r-l 

1 

T-i 

x0 

i 

en 
I 

r~* 
1 

rH 

r-4 
1 

< 

CO 

1 

in 
I 

CM 

i 

r~f 

rH 

1 

en 

1 

en 

»t 

l 

rH 

rH 

CNJ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

r-i 

CM 

CM 

CM 

t-i 

f~\ 

rH 

CM 

rH 

o 

1 

o 

1 

O 
I 

o 

1 

o 

1 

O 
i 

O 
i 

O 
i 

o 

1 

O 
i 

O 
I 

o 

1 

o 

1 

O 
■ 

O 
i 

o 

1 

o 

1 

o 

1 

O 

o 

l/> 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

l 

UJ 

F 

UJ 

1 

UJ 

i 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

LU 

1 

UJ 

1 

LU 

1 

LU 

1 

UJ 

»- 

rH 

1- 

r- 

r\j 

»- 

r- 

vt 

o 

CO 

CM 

m 

00 

r-< 

o 

o 

-0 

^0 

O 

o 

in 

z 

o 

CM 

o 

O 

o 

rH 

<\J 

en 

in 

in 

o 

o 

a> 

«fr 

in 

in 

r- 

rH 

en 

vO 

L1J 

o 

in 

o 

CO 

in 

00 

O 

o 

O 

O 

in 

00 

* 

CM 

CO 

en 

>J- 

tf 

o 

en 

r-< 

<»- 

en 

00 

O 

CO 

en 

r- 

■4- 

o 

r- 

rA 

CM 

co 

r- 

r- 

CM 

^ 

sO 

CM 

in 

o 

o 

o 

►—I 

o 

vO 

o 

o 

CM 

en 

en 

m 

o 

o> 

0s 

o 

o 

vj" 

in 

H 

en 

rH 

rH 

en 

Cr> 

CM 

LL 

in 

1 

1 

1 

o 

1 

1 

1 

1 

UL 

• 

LU 

O 

o 

II 

II 

ii 

ii 

ii 

ii 

ii 

ii 

ii 

ti 

l-< 

n 

II 

ii 

II 

ii 

n 

II 

ti 

II 

ii 

ii 

o 

r-! 

C\J 

co 

vt- 

in 

«Q 

Is- 

00 

o 

II 

o 

rH 

CM 

en 

<t- 

in 

o 

r- 

CO 

o 

o 

rH 

►H 

r-l 

rH 

rH 

r-i 

i-h 

r-f 

rH 

rH 

o 

rH 

r-A 

rH 

H 

r-l 

r~i 

rH 

rH 

rH 

rH 

z 

X 

Vj. 

«^ 

%_* 

w 

w 

*-* 

■wr 

w. 

^-- 

w 

X 

w « 

»-« 

*^ 

%-* 

^» 

V-* 

•^ 

«-^ 

*m* 

<-» 

o 

CL 

_l 

_l 

_l 

_J 

_l 

-J 

_J 

_l 

-J 

_J 

cc 

„J 

_J 

_J 

_J 

_J 

_l 

_J 

_J 

_J 

_J 

•— < 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

LU 

LU 

1- 
u 

o 

Q 

O 

o 

O 

Q 

Q 

a 

D 

a 

o 

o 

a 

o 

o 

Q 

O 

Q 

o 

Q 

u_ 

o 

rH 

r-l 

rH 

r~4 

CM 

CM 

CM 

rH 

rH 

r-i 

r-i 

r-t 

CM 

t-i 

l-{ 

l-l 

rH 

rH 

r-t 

UJ 

o 

o 

o 

1 

O 
1 

o 

1 

O 
I 

O 
i 

O 
1 

o 

o 

1 

O 
I 

O 
| 

O 
1 

o 
i 

o 

1 

o 

1 

o 

o 

1 

o 

1 

O 
1 

o 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 
UJ 

1 

UJ 

i 

UJ 

i 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

i 

UJ 

1 

UJ 

1 

UJ 

Ul 

1 

UJ 

1 

UJ 

1 

LU 

1 

UJ 

LU 

o 

CM 

o 

rH 

o 

h- 

in 

Is- 

•4- 

o 

CM 

o 

0s- 

en 

in 

CM 

en 

r- 

r- 

CO 

s0 

•4" 

on 

o 

h- 

r-l 

r- 

in 

v0 

en 

rH 

r~< 

in 

o 

o 

o 

o 

vj- 

0> 

CM 

o 

CM 

o 

in 

o 

UJ 

m 

o 

vt- 

m 

r- 

CM 

CM 

o 

o> 

rH 

«* 

in 

en 

in 

m 

h- 

o 

CM 

>o 

r-i 

r- 

0^ 

u. 

• 

0s- 

h- 

*fr 

CO 

CM 

r- 

r-4 

vO 

f- 

T-\ 

• 

0> 

r-4 

CM 

-o 

<f 

CM 

O 

CM 

<r 

LO 

a: 

o 

o 

UJ 

r-< 

h- 

4: 

1— 1 

rH 

vO 

en 

rH 

rH 

T-* 

CO 

in 

CM 

o 

rH 

i-\ 

en 

en 

-4- 

-t 

h- 

ii 

1 

1 

1 

II 

1 

1 

1 

1 

z 

:*: 

*- < 

D 

r-l 

rH 

CM 

CNJ 

CM 

CM 

CM 

CM 

CM 

CM 

D 

CM 

r-< 

CM 

CM 

CM 

rH 

rH 

r-t 

CM 

r~* 

O 
| 

o 

1 

o 

1 

o 

1 

O 
■ 

o 

1 

O 
1 

o 

1 

O 
i 

o 

o 

1 

o 

1 

o 

1 

O 
1 

O 
1 

O 
1 

o 

1 

O 
i 

o 

1 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

i 

UJ 

1 

LU 

1 

UJ 

1 

UJ 

i 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

LU 

1 

UJ 

1 

LU 

1 

LU 

l 

U_i 

1 

UJ 

r-4 

h- 

r- 

<M 

h- 

f- 

>t 

O 

CO 

CM 

in 

co 

<-A 

o 

O 

s0 

s0 

<G 

o 

in 

o 

cm 

o 

0> 

o 

r-t 

CM 

en 

in 

in 

o 

o 

O 

>t 

m 

in 

Is- 

rH 

en 

o 

o 

in 

o 

CO 

in 

CO 

0s- 

o 

«o 

0^ 

in 

00 

sf 

CM 

CO 

en 

^f 

<* 

o 

en 

h 

■4- 

en 

co 

O 

co 

en 

r- 

«4- 

O 

f- 

rH 

CM 

00 

Is- 

r- 

CM 

* 

v0 

CM 

m 

> 

vD 

0s- 

O 

CNJ 

en 

en 

in 

n0 

0s- 

o> 

o 

«* 

in 

rH 

cn 

t— i 

r-4 

en 

0> 

CM 
■ 

H 

1 

1 

1 

1 

1 

I 

I 

1 

1 

1 

1 

1 

1 

tJ 

II 

ii 

11 

II 

ii 

n 

ii 

ii 

II 

ii 

II 

II 

ll 

II 

ii 

ii 

II 

n 

ii 

II 

CQ 

o 

< 

t 

o 

r-l 

CM 

fO 

<j- 

in 

o 

Is- 

00 

0s- 

t 

o 

i-t 

CM 

en 

«4- 

in 

o 

r- 

CO 

0> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

ii 

_l 

_l 

_l 

_J 

_j 

_j 

_J 

_j 

_J 

_j 

II 

-J 

_J 

_J 

_j 

_i 

_j 

_l 

_i 

_j 

_J 

_j 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

_J 

UJ 

UJ 

UJ 

UJ 

LU 

LU 

LU 

LU 

LU 

UJ 

Q 

o 

O 

o 

o 

o 

a 

a 

a 

o 

O 

o 

o 

a 

O 

o 

a 

o 

o 

o 

O 

o 

69 


O 

^M* 

.— 

<-. 

#""* 

O 

r-< 

-H 

CM 

r-< 

0 

O 

O 

rH 

O 

0 

0 

rH 

rH 

O 

0 

0 

0 

rH 

rH 

O 

O 

in 

O 

O 

O 

O 

0 

O 

O 

O 

O 

0 

0 

O 

O 

O 

0 

0 

0 

O 

O 

O 

O 

rH 

1 

1 

1 

1 

CM 

1 

1 

O 

UJ 

LU 

lU 

LLI 

UJ 

LU 

LU 

UJ 

UJ 

UJ 

O 

UJ 

LU 

LU 

UJ 

UJ 

UJ 

LU 

UJ 

LU 

UJ 

• 

0^ 

co 

ro 

0") 

r~l 

Sf 

CO 

O 

r~ 

<t 

» 

r-\ 

rH 

O 

CM 

0 

in 

CM 

CO 

rH 

CO 

O 

CM 

sf 

C\J 

CO 

in 

CO 

CM 

in 

^0 

0 

O 

rH 

0> 

O 

CO 

CO 

r- 

CO 

CO 

CM 

co 

LO 

CT' 

O 

CO 

CO 

O 

O 

r-l 

CM 

CO 

rH 

CO 

sf- 

CO 

0 

co 

sT 

rH 

in 

CM 

II 

vO 

CM 

sf- 

r~ 

rH 

CM 

CM 

r-l 

0 

CO 

II 

rH 

c> 

O 

CO 

CO 

CO 

CO 

CO 

r- 

O 

-J 

_J 

< 

in 
1 

1 

m 

sj- 

r-l 

sj- 

1 

1 

r-t 

>* 

CO 

1 

< 

CO 

1 

st- 

1 

1 

f-i 

r~ 

CO 

1 

1 

rH 

r- 

vO 

1 

CM 

t— t 

r-( 

r-l 

r-f 

0 

O 

O 

rH 

0 

rH 

T-\ 

r-t 

0 

0 

r-i 

O 

rH 

O 

rH 

O 
1 

0 
1 

0 
1 

O 
■ 

O 
I 

0 

O 

O 

O 

0 

O 
1 

O 
1 

O 
1 

0 

0 

O 

O 

O 

O 

O 

CO 

1 

LU 

1 

LU 

1 

LU 

1 

UJ 

1 

UJ 

LU 

LU 

UJ 

1 

LU 

LU 

1 

lU 

1 

LU 

1 

LU 

UJ 

UJ 

UJ 

LU 

LU 

UJ 

LU 

h- 

h- 

in 

S?" 

0 

r-l 

CO 

rH 

CO 

rH 

r- 

r- 

CO 

CO 

(O 

in 

CM 

CO 

r- 

0 

O 

z 

in 

CM 

0 

rH 

in 

0 

SJ- 

r-l 

CO 

0 

0 

CO 

O 

0 

co 

CO 

CT" 

CM 

CO 

O 

uu 

O- 

co 

r* 

co 

C\l 

rH 

in 

O 

CM 

0 

in 

in 

sf 

in 

CM 

CO 

CO 

CM 

r-l 

c- 

y- 1 

vO 

cm 

in 

Ln 

CO 

CO 

CO 

CM 

CM 

0 

CO 

sj- 

CO 

CM 

co 

CO 

r- 

0 

CM 

0 

O 

O 

O 

»— i 

O 

«fr 

CM 

r-l 

CO 

r- 

CM 

r-l 

in 

CO 

CO 

O 

rH 

rH 

sf 

si- 

r-t 

rH 

CM 

CO 

t-A 

<-\ 

II. 

in 

1 

1 

1 

1 

1 

O 

1 

1 

1 

1 

1 

1 

LL 

• 

ID 
O 
O 

f-H 

II 

11 

11 

11 

11 

II 

II 

11 

II 

II 

CM 

II 

II 

II 

ll 

II 

II 

II 

II 

II 

II 

II 

0 

f-( 

CM 

CO 

st- 

m 

>£> 

r- 

CO 

0 

II 

O 

r-t 

CM 

CO 

^ 

in 

vO 

h- 

CO 

0^ 

O 

r-t 

r-l 

r-< 

rH 

rH 

r-l 

rH 

1-1 

rH 

rH 

O 

rH 

rH 

rH 

rH 

r~< 

f-t 

r-{ 

rH 

rH 

r-i 

2T 

X 

»»«• 

«— ' 

^o 

%** 

«-* 

— 

w 

w* 

«-» 

«^» 

X 

**-- 

«-» 

s^ 

w* 

^»* 

%—* 

<w^ 

w^ 

+-» 

^_ 

O 

a: 

_l 

_J 

_J 

-J 

-_J 

•_J 

_J 

-J 

_J 

-J 

cc. 

_J 

_J 

_l 

_J 

_J 

_j 

_J 

_J 

_J 

_J 

t—t 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

UJ 

UJ 

LU 

LU 

UJ 

UJ 

LU 

UJ 

lU 

UJ 

UJ 

O 

D 

O 

O 

Q 

O 

D 

O 

O 

0 

0 

O 

O 

0 

0 

O 

a 

a 

0 

0 

2 
ID 

U- 

rH 

H 

CM 

r-l 

O 

O 

0 

r-l 

O 

0 

l-t 

rH 

O 

0 

0 

O 

rH 

rH 

0 

0 

UJ 

0 
2: 

O 
1 

O 

1 

O 
1 

O 
1 

O 

O 

0 

O 

O 

0 

0 

1 

O 
1 

O 

0 

0 

O 

O 

O 

0 

0 

1 

LU 

1 

LU 

1 

LU 

1 

LU 

UJ 

UJ 

LU 

LU 

LU 

lU 

1 

LU 

1 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

UJ 

0 

O 

ro 

CO 

CO 

r-i 

Sf- 

CO 

0 

t-- 

S»- 

0 

rH 

rH 

O 

CM 

O 

in 

CM 

co 

r-A 

CO 

a: 

0 

CM 

si- 

CM 

CO 

in 

co 

CM 

in 

-O 

O 

0 

rH 

O 

O 

CO 

CO 

r- 

CO 

CO 

CM 

co 

LU 

m 

in 

O 

O 

co 

co 

0 

O 

r-l 

CM 

CO 

in 

rH 

CO 

sr 

co 

O 

CO 

sf 

r~i 

in 

CM 

LL 

• 

vO 

CM 

sT 

r- 

r-l 

CM 

CM 

rH 

O 

CO 

• 

r-i 

& 

0 

CO 

01 

CO 

ro 

CO 

r- 

0> 

Ql 

0 

O 

LU 

m 

*t 

in 

st- 

r-l 

>t 

<f- 

rH 

sf- 

co 

CO 

vf- 

rH 

rH 

h- 

CO 

rH 

rH 

r- 

O 

h- 

11 

1 

1 

1 

1 

1 

ll 

1 

1 

1 

1 

z 

^ 

t— 1 

0 

cm 

r- I 

r-H 

rH 

r-l 

O 

O 

O 

rH 

0 

O 

r-i 

rH 

rH 

0 

0 

rH 

O 

rH 

0 

■— 1 

0 
1 

O 
1 

O 

| 

O 
1 

O 

O 

O 

O 
1 

0 

O 
1 

O 
1 

O 
1 

0 

0 

O 

O 

O 

0 

O 

1 

LU 

1 

LU 

LU 

1 

LU 

1 

LU 

UJ 

UJ 

LU 

1 

UJ 

LU 

1 

LU 

1 

UJ 

1 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

UJ 

LU 

H 

h- 

in 

St" 

O 

r-l 

CO 

rH 

CO 

T-i 

r- 

r- 

co 

co 

CO 

in 

CM 

CO 

r- 

0 

O 

H 

m 

CM 

O 

r-l 

in 

0 

SI- 

rH 

CO 

vO 

>o 

CO 

0 

0 

CO 

CO 

c> 

CM 

CO 

O 

> 

r- 

ro 

h- 

CO 

CM 

r-l 

LO 

vO 

CM 

0 

in 

in 

st- 

in 

CM 

CO 

CO 

CM 

r-K 

r- 

0 

CM 

in 

in 

CO 

CO 

CO 

CM 

CM 

0 

CO 

<* 

CO 

CM 

co 

CO 

h- 

O 

CM 

vO 

W 

< 

sf 

CM 

1 

rH 

CO 

I'- 

CM 

r-l 

in 

CO 

CO 

rH 

rH 

sf 

<f 

rH 

rH 

CM 

CO 

rH 

rH 

1 

1 

1 

1 

1 

1 

1 

1 

Eh 

0 

11 

II 

II 

11 

ll 

II 

II 

11 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

• 

0 

r-l 

CM 

CO 

-t 

in 

O 

r- 

CO 

O 

• 

O 

f-i 

CM 

CO 

sl- 

in 

0 

r- 

CO 

O 

0 

0 

O 

O 

0 

0 

0 

O 

0 

0 

0 

O 

O 

O 

0 

0 

0 

0 

0 

O 

O 

O 

11 

_j 

_J 

_J 

_j 

_i 

_i 

_J 

_i 

_J 

_J 

ll 

_J 

_l 

_J 

_J 

_J 

_j 

_l 

_i 

_) 

_l 

_j 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

UJ 

LU 

_i 

LU 

LU 

UJ 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

0 

Q 

O 

Q 

0 

O 

O 

O 

O 

a 

0 

O 

O 

O 

0 

a 

0 

a 

0 

O 

O 

a 

70 


O 

O 

rH 

i— < 

o 

o 

o 

r~l 

o 

r-t 

o 

r-l 

o 

CM 

ro 

i—l 

H 

i—t 

rH 

1—1 

o 

l-t 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

in 

1 

l 

O 

1 

1 

1 

1 

1 

1 

1 

1 

O 

UJ 

UJ 

LU 

UJ 

IU 

UJ 

LU 

UJ 

UJ 

LU 

i—t 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

• 

o 

<f- 

CO 

in 

r- 

in 

vt 

ro 

00 

r- 

• 

in 

co 

v* 

0s- 

CM 

CO 

i— 1 

(M 

o 

in 

o 

CM 

o 

h- 

<r 

h> 

0s- 

CO 

O 

oo 

00 

o 

rH 

CM 

0s- 

r- 

O 

rH 

CO 

m 

CO 

in 

00 

vO 

f- 

CM 

o 

Is- 

in 

O 

CO 

in 

0> 

vO 

Is- 

o 

<r 

0s- 

-t 

CM 

vO 

ro 

II 

in 

ro 

0s- 

CO 

vO 

O 

CO 

m 

o 

•4- 

II 

in 

rH 

ro 

in 

<i- 

o 

o 

O 

vO 

CO 

_J 

_l 

< 

i— i 
1 

in 
1 

r-t 

in 

ro 

I 

r-t 
1 

r-t 

CM 

f* 

r-t 
1 

<r 

r-t 
1 

vO 

rH 

1 

i 

CO 

1 

vO 

CO 

1 

co 
1 

<M 

rH 

o 

o 

o 

r-l 

o 

i—i 

o 

rH 

CM 

CM 

o 

o 

i— t 

T-i 

rH 

rH 

r-A 

rH 

O 
i 

O 
1 

o 

o 

o 

o 

o 

o 

o 

O 

O 
i 

o 

1 

o 

o 

o 

o 

O 

o 

o 

O 

oo 

I 
UJ 

1 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

U) 

I 

LU 

1 

LU 

LU 

LU 

UJ 

LU 

LU 

UJ 

UJ 

UJ 

h- 

0s- 

o 

vt 

0s 

Is- 

ro 

0s 

vO 

00 

Is- 

r~ 

CM 

O 

o 

CM 

rH 

CM 

CM 

r- 

o 

•z. 

CM 

>t 

in 

CM 

tji 

•* 

ro 

co 

r- 

rH 

0s 

CO 

in 

•4- 

CM 

CM 

oo 

CM 

cO 

0s" 

LU 

co 

CM 

CO 

»* 

r- 

0s- 

CO 

in 

o 

in 

o 

CO 

r-t 

o 

CM 

s0 

ro 

rH 

ro 

^t 

»— i 

co 

CO 

CM 

ro 

ro 

0s- 

vO 

0s- 

ro 

in 

o 

CM 

ro 

ro 

CO 

<*- 

rH 

o 

r- 

r- 

o 

o 

o 

• 

• 

• 

• 

e 

t 

• 

* 

« 

t 

t— « 

o 

r-t 

o 

CM 

o 

r- 

rH 

vO 

CM 

r-t 

rH 

o 

«fr 

& 

CM 

ro 

<~i 

CM 

<f 

in 

CM 

CO 

U_ 

o 

1 

l 

I 

1 

1 

o 

1 

1 

1 

1 

u. 

• 

UJ 

in 

II 

II 

II 

n 

ii 

II 

II 

li 

ii 

It 

o 

II 

II 

II 

ii 

II 

II 

II 

ii 

II 

ii 

o 

o 

n 

o 

rH 

CM 

CO 

■sJ- 

in 

O 

Is- 

co 

0s- 

II 

o 

rH 

CM 

ro 

<f 

in 

sO 

N- 

CO 

0s- 

o 

1— t 

i-< 

r-t 

rH 

l~l 

r-t 

r-t 

r-t 

r-t 

rH 

o 

r-t 

rH 

r-t 

rH 

r-i 

rH 

rH 

rH 

rH 

rH 

z 

X 

V** 

s.^ 

«-» 

w 

*»« 

<~» 

v,< 

>~» 

»_» 

»-» 

X 

«— ' 

w 

^» 

w 

w* 

^< 

«_« 

»-• 

w* 

..  * 

o 

a: 

_J 

_J 

_J 

-J 

_l 

_l 

_J 

_l 

_J 

_J 

a: 

_J 

_J 

_J 

_J 

_J 

_J 

_J 

-J 

_J 

-J 

r-t 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

UJ 

LU 

LU 

UJ 

o 

2 

Q 

o 

o 

D 

o 

O 

Q 

a 

o 

O 

o 

O 

Q 

Q 

Q 

a 

Q 

O 

a 

Q 

U- 

rH 

l-l 

o 

o 

o 

r-t 

O 

r-t 

o 

rH 

CM 

ro 

rH 

rH 

i—t 

rH 

•— l 

o 

rH 

O 

lu 

O 

o 

o 

o 

o 

o 

O 

O 

o 

o 

O 

o 

o 

o 

O 

o 

O 

o 

O 

o 

o 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

z 

UJ 

IU 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

UJ 

LU 

LU 

IU 

UJ 

o 

0> 

>d- 

oo 

in 

r- 

in 

st- 

ro 

CO 

Is- 

o 

in 

CO 

^ 

O 

CM 

CO 

rH 

t\J 

O 

in 

CC 

o 

CNJ 

o 

r- 

«t 

r- 

0s 

ro 

o 

CO 

00 

o 

rH 

CM 

o 

Is- 

O 

1—1 

CO 

in 

CO 

in 

UJ 

in 

00 

vO 

r- 

CM 

0s 

Is- 

in 

o 

CO 

in 

in 

0s- 

o 

Is- 

o 

«* 

o 

^J■ 

CM 

O 

ro 

LL 

t 

in 

ro 

o 

00 

o 

o 

ro 

in 

o 

>* 

• 

m 

rH 

ro 

in 

^ 

o 

vO 

0^ 

vO 

CO 

cl 

o 

o 

UJ 

f-l 

in 

r-t 

m 

ro 

r-l 

r-t 

(NJ 

r~t 

rH 

rH 

V0 

rH 

<r 

CM 

CO 

sO 

CO 

«t 

co 

h- 

H 

1 

1 

1 

1 

1 

n 

1 

1 

1 

1 

■z. 

^ 

t—t 

D 

CM 

r-t 

o 

o 

o 

r-t 

o 

r-t 

o 

r-t 

o 

CM 

CM 

o 

o 

r^ 

r~ 1 

rH 

r-i 

r-t 

i—t 

O 
t 

o 
■ 

o 

o 

o 

O 

o 

O 

o 

o 

i 

O 
1 

o 

o 

o 

o 

O 

o 

O 

o 

1 

UJ 

i 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

1 

LU 

1 

LU 

LU 

UJ 

LU 

LU 

UJ 

LU 

UJ 

UJ 

0s- 

o 

*J- 

0s 

r- 

CO 

o 

vO 

co 

Is- 

Is- 

CM 

O 

o 

CM 

rH 

<M 

C^J 

Is- 

o 

CM 

<t- 

in 

CM 

o> 

»J- 

ro 

CO 

Is- 

rH 

0s- 

CO 

in 

■4- 

CM 

(M 

CO 

(\J 

ro 

o> 

H 

ro 

CM 

00 

>t 

r» 

0s- 

00 

in 

o 

in 

vO 

00 

r-i 

O 

CM 

vO 

ro 

r-i 

ro 

«4- 

> 

co 

00 

CM 

ro 

ro 

0s- 

o 

0s 

ro 

in 

o 

CM 

ro 

ro 

CO 

>t 

r-t 

o 

Is- 

Is- 

w 

< 
Eh 


O 

o 

II 

_J 
o 


rH 

0s- 

1 

CM 
1 

O 

Is- 

rH 
1 

I 

CM 

rH 

rH 
1 

vf 

0s- 

1 

CM 

ro 

i— i 
1 

CM 
1 

>* 

in 

CM 
1 

ro 
1 

II 

II 

II 

II 

II 

II 

ii 

11 

II 

II 

II 

II 

II 

ii 

ii 

II 

ii 

ii 

ii 

ii 

o 

1—t 

CM 

ro 

-* 

in 

^0 

Is- 

00 

o 

• 

o 

r-t 

CM 

ro 

** 

in 

o 

r- 

00 

0s- 

o 

o 

O 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

o 

o 

LU   LU 
Q   O 


LU   LU   IU 

a  o  a 


UJ   UJ 


UJ 

Q 


LU 

a 


LU 

a 


ii 

_i 
o 


LU   LU 

a  a 


LU   LU   LU 

o  o  a 


UJ   UJ   LU 

Q   Q   O 


UJ 

O 


LU 
Q 


71 


O 

o 

o 

f— 1 

rH 

1— 1 

rH 

CM 

r-t 

CM 

r-l 

r-t 

O 

t-t 

r-4 

r-t 

r-l 

rH 

o 

o 

o 

o 

O 

LO 

o 

o 

o 

o 

O 

O 

o 

O 

o 

O 

O 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i— t 

1 

1 

1 

1 

1 

o 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

o 

UJ 

UJ 

LU 

LU 

LU 

LU 

UJ 

UJ 

LU 

LU 

• 

o 

o 

o 

CM 

CO 

o 

<£) 

CO 

■+ 

r- 

• 

o 

in 

<f 

in 

1*- 

co 

CM 

■* 

CM 

<f 

o 

o 

o 

o 

CM 

o 

CM 

CO 

m 

r- 

o 

o 

Q> 

h- 

CM 

r-t 

o 

o 

co 

CM 

h- 

r~i 

<«• 

o 

vO 

r-t 

o 

r- 

CO 

CM 

CM 

v0 

O 

-* 

CM 

r*- 

<t- 

O 

CO 

r-t 

o 

■4- 

II 

<f 

s0 

CD 

r-t 

in 

m 

CO 

o 

o 

in 

II 

in 

o 

in 

M" 

m 

r- 

CM 

n0 

o 

r-i 

_J 

_J 

< 

rH 

1 

in 
1 

<?• 

CM 

1 

F-l 

o 
1 

r-i 

1 

0^ 

CO 

CM 

1 

< 

1 

in 
1 

r-l 

CM 

o 

CM 

1 

CO 
1 

r- 

CO 

o 
1 

CM 

i— t 

•— t 

CM 

ro 

CM 

CM 

r-t 

CM 

CM 

l—l 

r-t 

l—l 

rH 

r-i 

r~< 

r-t 

o 

rH 

1— 1 

o 

1 

o 

1 

o 

1 

o 

1 

o 

1 

O 
i 

O 
i 

o 

1 

o 

1 

O 

i 

o 

1 

O 
1 

o 

O 
l 

a 

o 

1 

O 
1 

o 

o 

1 

o 

1 

CO 

1 

UJ 

1 

UJ 

1 

UJ 

1 
UJ 

1 
UJ 

1 

LU 

1 

UJ 

1 

UJ 

1 

UJ 

1 

LU 

1 
LU 

1 

LU 

1 

LU 

1 

UJ 

i 

LU 

1 
LU 

1 

UJ 

UJ 

1 

LU 

1 

LU 

»- 

<1- 

in 

i-H 

CM 

0> 

r-t 

>4- 

O 

l»- 

O 

l—l 

CM 

f- 

in 

CO 

in 

rH 

o 

o> 

r- 

z 

o 

i-i 

co 

CO 

o 

CM 

CO 

o 

r-l 

o 

vO 

CM 

0> 

r-t 

4- 

in 

CM 

o 

CO 

m 

UJ 

O 

-o 

CO 

O 

o 

O 

h- 

r> 

in 

in 

CM 

in 

o 

o 

CO 

CM 

h" 

o 

O 

o 

t— 1 

vO 

o 

CO 

r- 

in 

f— I 

00 

co 

o 

i— t 

«t 

CM 

CM 

o 

CO 

r-t 

O 

M" 

in 

a> 

o 

o 

o 

r-t 

o 

in 

o 

r\j 

CM 

rH 

cO 

in 

l—l 

r-t 

CO 

o 

CM 

CM 

CM 

l—l 

l—l 

h- 

rH 

rH 

s}- 

CM 

U- 

in 

1 

I 

1 

1 

o 

1 

1 

1 

1 

U- 

• 

UJ 

o 

ii 

II 

ll 

ii 

II 

II 

ii 

II 

II 

II 

l—l 

II 

ii 

II 

II 

II 

II 

II 

II 

II 

ii 

o 

**°* 

^, 

o 

II 

o 

rH 

CM 

co 

4- 

in 

NO 

f- 

CO 

o 

II 

O 

r-4 

CM 

CO 

0* 

in 

O 

r- 

CO 

o 

o 

i~» 

»H 

i-i 

r-l 

r-A 

r-t 

r-l 

r-l 

I— 1 

H 

D 

r-l 

r-t 

rH 

r-t 

r~\ 

r-t 

r-t 

H 

1— 1 

rH 

z 

X 

«.- 

*^* 

*■»* 

w 

«— 

w 

«— ' 

>— ' 

w 

^^ 

X 

v^  - 

«— 

w* 

w- 

«_> 

«— > 

*-^ 

— < 

w 

tarf 

o 

c£ 

_i 

_J 

-j 

_J 

_l 

_J 

_J 

_J 

_J 

_J 

a: 

-J 

_l 

_l 

_J 

_J 

_) 

_J 

_l 

—1 

^ 

t— 1 

UJ 

UJ 

LU 

LU 

LU 

UJ 

UJ 

lU 

LU 

LU 

UJ 

LU 

U) 

LU 

UJ 

LU 

UJ 

LU 

LU 

UJ 

1- 
o 

o 

O 

o 

o 

Q 

O 

Q 

Q 

Q 

O 

a 

O 

o 

o 

o 

o 

a 

a 

a 

o 

U_ 

o 

r-4 

r-t 

r-l 

rH 

CM 

r-t 

CM 

r-l 

r-t 

r-t 

r-l 

r-l 

l-t 

rH 

o 

o 

o 

o 

o 

UJ 

o 

O 

o 

O 

o 

O 

o 

O 

O 

o 

a 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

z 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

lU 

UJ 

LU 

LU 

LU 

UJ 

LU 

UJ 

LU 

LU 

LU 

LLI 

LU 

LLt 

o 

o 

>o 

o 

CM 

CO 

o 

o 

CO 

■4- 

r» 

o 

vO 

in 

^4- 

in 

r- 

CO 

CM 

vt 

CM 

>t 

OC 

o 

o 

o 

o 

CM 

o 

CM 

CO 

in 

r- 

vO 

o 

o 

r- 

CM 

i—i 

vO 

CJ> 

CO 

CM 

r- 

r-t 

UJ 

o 

<t- 

o 

o 

r-t 

o 

r- 

CO 

CM 

CM 

vO 

o 

o 

4- 

CM 

h- 

^4- 

o 

CO 

rH 

>o 

4- 

u. 

• 

•4" 

o 

co 

f-l 

in 

in 

CO 

O 

o 

in 

• 

in 

o 

m 

vj- 

in 

r- 

CM 

<0 

o 

r-t 

C* 

rH 

• 

• 

« 

* 

• 

• 

• 

• 

• 

• 

r-l 

UJ 

T-\ 

in 

^ 

CM 

r-t 

0^ 

r-t 

O 

co 

CM 

v0 

in 

i—i 

CM 

o> 

(M 

CO 

r- 

CO 

o 

H 

II 

1 

1 

1 

1 

II 

1 

1 

I 

l 

1 

Z 

x: 

^ 

•— i 

o 

CM 

r-t 

f— 1 

CM 

co 

CM 

CM 

r-t 

CM 

(M 

o 

•—I 

r-l 

r-t 

i—l 

r-t 

r-l 

•—I 

o 

r-l 

f— 1 

O 

UJ 

O 
1 

o 

1 

O 
i 

o 
i 

O 

o 

1 

o 

1 

O 
i 

O 

1 

o 

o 

1 

O 
1 

O 
I 

O 
1 

o 

1 

o 

1 

o 

o 

1 

o 

1 

1 

UJ 

1 

UJ 

1 

UJ 

i 

UJ 

1 

UJ 

1 

LU 

1 

UJ 

1 

UJ 

l 

LU 

i 

UJ 

1 

LU 

1 

UJ 

1 

lU 

1 

LU 

UJ 

1 

UJ 

UJ 

1 

UJ 

1 

UJ 

<J- 

in 

rH 

CM 

o> 

1-1 

«t 

o 

h- 

O 

r-l 

CM 

r- 

in 

CO 

in 

l—l 

0^ 

0> 

r- 

H 

& 

r-t 

CO 

ro 

0> 

CM 

CO 

o 

r-t 

o 

vO 

CM 

o 

r-l 

** 

in 

CM 

<J< 

00 

in 

H 

> 

O 

o 

CO 

■o 

o 

o 

h- 

r^ 

in 

in 

CM 

in 

o 

vO 

CO 

CM 

r- 

o 

vO 

o 

vO 

& 

CO 

t^- 

in 

r-l 

CO 

CO 

r- 

i— t 

sh 

CM 

CM 

^O 

CO 

r-l 

o 

«*■ 

m 

o 

in 

o 

CM 

CM 

r-l 

CO 

in 

r-l 

r-t 

co 

CM 

CM 

CM 

r-l 

r-t 

r- 

rH 

1— 1 

-t 

CM 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

< 

ii 

II 

II 

II 

II 

II 

ii 

II 

II 

ii 

II 

II 

II 

II 

II 

n 

II 

II 

ll 

II 

^ 

o 

f 

o 

rH 

CM 

co 

«4" 

in 

o 

h- 

CO 

o 

• 

o 

r-l 

CM 

CO 

>t 

in 

vO 

r- 

co 

a> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

II 

_i 

_l 

_J 

_i 

_J 

_j 

_J 

-J 

_l 

_J 

II 

_J 

-J 

_J 

_l 

_J 

_j 

_J 

_j 

_J 

_j 

_J 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

-J 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

o 

o 

o 

o 

Q 

o 

o 

o 

O 

Q 

Q 

o 

o 

a 

D 

D 

O 

O 

O 

o 

Q 

O 

72 


O 

o 

CM 

O 

o 

r-H 

rH 

r-< 

rH 

rH 

rH 

T-* 

o 

rH 

o 

o 

O 

o 

o 

O 

o 

o 

o 

LO 

o 

O 

o 

o 

O 

O 

o 

O 

O 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

rH 

1 

CM 

1 

o 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

O 

LU 

LU 

LU 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

• 

vO 

C\J 

O 

t— 1 

ro 

vO 

o 

o 

Sf 

r- 

• 

rM 

rM 

o 

CM 

sh- 

CO 

CM 

rH 

o 

o 

o 

o 

-H 

CM 

0^ 

CM 

CM 

CO 

o 

f» 

o 

o 

0s 

<t 

o 

rH 

oo 

co 

CM 

sf 

CM 

0> 

CM 

O 

i-t 

r-l 

Sf 

s0 

ro 

r~i 

C\J 

r«- 

o 

rH 

o 

CO 

r-l 

o 

O 

o 

CO 

o 

II 

-4- 

r- 

ro 

Sf 

Is- 

Is- 

O 

ro 

rH 

o 

II 

ro 

CO 

o 

vO 

CM 

CO 

0> 

vO 

Sf 

rH 

_J 

• 

• 

• 

• 

c 

• 

t 

• 

• 

» 

_J 

• 

* 

t 

• 

« 

• 

0 

0 

• 

• 

< 

CO 

1 

rH 
I 

i 

r-l 

CO 

Sf 
1 

vO 

1 

v£> 

ro 

CM 
1 

< 

sf 

rH 
1 

in 
I 

<|- 

o 

SD 

1 

Is- 
1 

CO 

ro 

in 
1 

rH 

o 

o 

r-4 

rH 

rH 

rH 

CM 

rH 

rH 

rH 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 
1 

o 

o 

O 

o 

O 

o 

O 

o 

o 

O 
1 

o 

o 

o 

O 

o 

o 

o 

o 

o 

00 

i 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

1 

UJ 

IU 

LU 

UJ 

UJ 

UJ 

LU 

UJ 

UI 

LU 

Y- 

O 

o 

CM 

CM 

v0 

O 

Sf 

r- 

sf 

rH 

o 

o 

in 

CM 

vO 

sf 

<f 

CM 

o 

CM 

•ZL 

O 

in 

O 

CM 

CM 

h- 

CM 

Is- 

in 

CM 

rH 

CO 

h- 

in 

o 

in 

rH 

o 

CO 

CM 

LU 

O 

to 

r»» 

CO 

CO 

ro 

O 

o 

O 

rM 

CM 

rH 

r- 

o 

ro 

rH 

CO 

in 

CM 

CO 

t-H 

in 

CO 

s0 

r-l 

rH 

ro 

CM 

o 

in 

o 

CM 

ro 

rH 

(M 

Is- 

ro 

o 

in 

O 

CM 

o 

o 

o 

t—4 

o 

Sf 

t— i 

<f 

CM 

CM 

Is- 

ro 

rH 

i— i 

in 

o 

r- 

CM 

ro 

r- 

^d- 

vO 

rH 

CM 

rH 

rH 

LL 

in 

1 

1 

1 

1 

1 

1 

o 

I 

1 

1 

1 

» 

U_ 

• 

• 

UJ 

O 
o 

i-i 

ii 

ii 

ii 

II 

II 

ll 

ll 

II 

ii 

ii 

CM 

ii 

n 

ii 

ii 

II 

ii 

II 

ii 

II 

II 

ii 

o 

rH 

CM 

ro 

sf 

in 

O 

r- 

co 

o 

II 

o 

rH 

CM 

ro 

sf 

in 

nO 

Is- 

CO 

0> 

o 

rH 

r-l 

r-; 

rH 

rH 

rH 

rH 

H 

«-i 

rH 

O 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

2 

X 

X 

•_• 

«~» 

<a* 

v^ 

ta« 

M 

*— 

w 

V-* 

^  * 

o 

aL 

_» 

_J 

_J 

_J 

_J 

_J 

_J 

_J 

_i 

_J 

a: 

_l 

_l 

_l 

_J 

_J 

_J 

_J 

_J 

_J 

_J 

»— < 

LU 

LU 

LU 

LU 

IU 

LU 

LU 

LU 

UJ 

UJ 

LU 

LU 

LU 

LU 

UJ 

LU 

UJ 

LU 

IU 

UJ 

I- 
o 

O 

O 

Q 

o 

Q 

O 

Q 

O 

Q 

Q 

o 

O 

O 

O 

Q 

Q 

O 

o 

o 

a 

Z> 

u_ 

cm 

O 

o 

r-l 

rH 

rH 

rH 

rH 

rH 

r-t 

rH 

o 

O 

o 

O 

o 

o 

o 

o 

o 

IXJ 

o 
■ 

o 

O 

o 

o 

O 

O 

o 

o 

o 

o 

1 

o 

o 

o 

O 

o 

o 

o 

o 

o 

i 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

iu 

LU 

IU 

1 
UJ 

LU 

LU 

UJ 

LU 

UJ 

LU 

UJ 

LU 

LU 

LJJ 

O 

O 

CM 

o» 

r-l 

CO 

vO 

o 

0s 

Sf 

Is- 

o 

CM 

r-* 

O 

CM 

<f 

ro 

CM 

rH 

o 

o 

cc 

O 

O 

r- ( 

CM 

O 

r\j 

CM 

CO 

o 

Is- 

o 

o 

O 

<t 

O 

rH 

co 

CO 

CM 

Sf 

CM 

o 

LU 

O 

(M 

O 

r-t 

rH 

sf 

O 

ro 

rH 

CM 

Is- 

o 

O 

rH 

o 

CO 

rH 

O 

O 

o 

ro 

v0 

LL 

• 

ST 

f- 

ro 

Sf 

Is- 

r- 

o 

(O 

rH 

o> 

c 

CO 

CO 

o 

o 

CM 

CO 

0> 

o 

sf 

rH 

c£ 

rH 

* 

t 

t 

• 

• 

• 

« 

* 

• 

• 

rH 

iu 

oo 

rH 

Is- 

r-H 

CO 

<t- 

vO 

vO 

CO 

(M 

>± 

r-i 

in 

>t 

O 

hO 

Is- 

CO 

ro 

in 

h- 

II 

1 

1 

1 

1 

II 

1 

1 

1 

1 

1 

2 

X. 

>-H 

Q 

rH 

O 

o 

rH 

r-l 

rH 

rH 

CM 

rH 

rH 

Q 

i— i 

O 

o 

o 

O 

o 

O 

o 

o 

o 

o 

1 

o 

o 

O 

o 

O 

o 

O 

O 

o 

o 

1 

O 

o 

o 

o 

o 

O 

o 

o 

o 

1 

LU 

LU 

LU 

LU 

LU 

LU 

Ui 

Uj 

LU 

UI 

1 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

UJ 

o 

O 

CM 

CM 

O 

O 

sf 

f- 

Sf 

rH 

o 

vO 

in 

CM 

sO 

sf 

Sf 

CM 

o 

CM 

o 

in 

vO 

CM 

CM 

h- 

CM 

Is- 

in 

CM 

rH 

CO 

Is- 

in 

o 

in 

rH 

O 

CO 

CM 

H 

o 

CO 

r- 

CO 

CO 

co 

O 

o 

o 

CM 

CM 

r-< 

r- 

o 

ro 

r~t 

ro 

in 

CM 

co 

> 

in 

CO 

nO 

rH 

r-l 

CO 

rM 

o 

in 

O 

CM 

ro 

rH 

CM 

h- 

CO 

-O 

in 

O 

CM 

Sf 

i— t 

vf 

CM 

CM 

I'- 

ro 

i— i 

rH 

in 

r» 

rM 

ro 

Is- 

vf 

vO 

rH 

CM 

rH 

r-i 

W 

1 

| 

1 

I 

1 

| 

1 

1 

Y 

J 

ff) 

ii 

ii 

II 

II 

II 

ll 

ll 

ii 

II 

n 

n 

ll 

ii 

II 

II 

II 

II 

ll 

II 

II 

< 

O 

&H 

t 

o 

i-i 

CM 

ro 

Sf 

in 

^O 

Is- 

CO 

a> 

• 

o 

r-t 

CM 

CO 

Sf 

in 

vO 

Is- 

co 

0> 

O 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

O 

II 

—I 

_j 

_J 

_i 

_J 

_J 

_l 

_J 

_J 

_j 

n 

_i 

-J 

_J 

_) 

_J 

_i 

_J 

_j 

_l 

_J 

_l 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

_j 

LU 

LU 

LU 

LU 

LU 

LU 

IU 

LU 

LU 

LU 

o 

D 

o 

o 

O 

o 

a 

o 

o 

o 

o 

o 

a 

o 

Q 

o 

o 

Q 

o 

a 

a 

O 

73 


o 

o 

o 

H 

r-l 

r-l 

r-l 

rH 

rH 

r-H 

r-t 

r-H 

O 

r-H 

rH 

.H 

o 

o 

rH 

o 

rH 

o 

rH 

in 

o 

o 

O 

O 

o 

O 

O 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

o 

1 

I 

1 

1 

1 

1 

1 

1 

1 

rH 

l 

1 

1 

o 

UJ 

UJ 

UJ 

111 

LU 

UJ 

IU 

UJ 

UJ 

Ul 

O 

LU 

LU 

LU 

LU 

LU 

UJ 

UJ 

LU 

LU 

IU 

• 

r- 

On 

CO 

CO 

Is- 

CO 

o 

in 

i-H 

r-H 

» 

o 

in 

o 

CM 

co 

•t 

<t 

r- 

rH 

ON 

o 

r-t 

o 

H 

r-l 

0N 

o 

o 

o 

r-H 

r-H 

o 

Is- 

(O 

CO 

in 

rH 

in 

CO 

Is- 

CO 

rH 

CO 

CO 

o 

rH 

o 

Is- 

in 

sO 

in 

CM 

«t 

CO 

sO 

CO 

•o 

CO 

o 

rH 

CO 

r- 

II 

rH 

rf) 

(\J 

in 

o 

Is- 

co 

CO 

r-H 

«* 

II 

o 

o 

r- 

Is- 

o 

in 

Is- 

Is- 

CM 

o 

_J 

~1 

• 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

< 

rH 
1 

1 

in 

C\l 

1 

CM 

r-l 
1 

CO 

1 

CO 

in 

in 
1 

< 

CM 
1 

CO 

» 

Is- 
1 

1 

V0 

CM 

CO 

1 

m 
1 

rH 

<f- 

rH 

rH 

r-l 

CM 

CM 

CM 

CO 

r-H 

rH 

CM 

rH 

H 

o 

o 

rH 

rH 

r-< 

CM 

H 

rH 

O 
| 

O 
I 

o 
■ 

O 
1 

O 
■ 

O 
1 

o 

1 

o 

1 

o 

1 

O 
l 

O 
■ 

o 

1 

o 

o 

O 

O 

o 

O 

o 

O 

l/> 

UJ 

1 

UJ 

i 

UJ 

i 
UJ 

1 

UJ 

i 

UJ 

1 

UJ 

1 

IU 

1 

UJ 

1 

LU 

1 

LU 

1 

IU 

IU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

1- 

rH 

in 

CM 

o 

CM 

r>- 

vO 

rH 

o 

CM 

in 

CO 

Is- 

in 

On 

CO 

St" 

co 

CM 

r-l 

2 

on 

Is- 

r-l 

r-l 

O 

co 

CO 

vO 

CO 

in 

rH 

vO 

rH 

rH 

CO 

rH 

CO 

Is- 

<(- 

sJ- 

UJ 

co 

CM 

c- 

o 

r- 

v0 

CO 

in 

CM 

•* 

co 

r-l 

CO 

eg 

CM 

CM 

in 

vj- 

rH 

* 

►—« 

in 

co 

r> 

CM 

vj" 

CM 

co 

o 

CO 

o 

& 

CO 

co 

rH 

<1- 

in 

o 

rH 

in 

CM 

o 

o 

o 

r— 1 

o 

<* 

in 

<C 

>0 

CO 

<j- 

CM 

rH 

rH 

r-l 

o 

in 

rH 

rH 

o 

rH 

in 

CO 

rH 

CM 

Is- 

u_ 

in 

1 

1 

1 

o 

1 

1 

1 

1 

1 

1 

u_ 

• 

UJ 

o 
o 

o 

n 

ii 

ii 

II 

II 

ii 

II 

II 

II 

II 

r-l 

ii 

II 

II 

II 

II 

ii 

ii 

II 

ii 

II 

II 

o 

rH 

CM 

CO 

<• 

in 

vO 

h- 

CO 

CT> 

II 

o 

rH 

CM 

ro 

>t 

in 

vO 

r- 

co 

ON 

o 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

f~\ 

rH 

o 

rH 

rH 

rH 

rH 

r-l 

rH 

rH 

rH 

rH 

rH 

2 

I 

*m*> 

«-- 

«-« 

X 

O 

OL 

_l 

-J 

_J 

_l 

-J 

_J 

_J 

_l 

_J 

_J 

CC 

_J 

_J 

-J 

-J 

_J 

_) 

_J 

_J 

_l 

_J 

►— i 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

UJ 

UJ 

LU 

UJ 

LU 

IU 

UJ 

UJ 

UJ 

LU 

LU 

UJ 

LU 

O 
2 

a 

O 

D 

o 

Q 

O 

o 

O 

D 

Q 

Q 

a 

o 

o 

a 

a 

a 

o 

Q 

O 

Z> 
UL 

o 

r-l 

rH 

rH 

r-l 

rH 

rH 

•— i 

r-l 

rH 

rH 

rH 

rH 

o 

o 

rH 

o 

rH 

O 

rH 

UJ 

O 
2 

o 

O 
1 

o 

1 

O 
I 

O 
1 

O 
1 

o 

o 

1 

o 

1 

O 
1 

o 

1 

O 

O 

1 

o 

o 

O 

o 

O 

O 

O 

LU 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

1 

IU 

1 

uj 

1 

UJ 

1 

IU 

1 

LU 

1 

UJ 

1 

LU 

IU 

UJ 

LU 

LU 

LU 

UJ 

LU 

UJ 

o 

Is- 

0> 

CO 

CO 

Is- 

CO 

o 

in 

r-i 

r-< 

o 

O 

in 

O 

CM 

CO 

«* 

•*■ 

Is- 

rH 

ON 

ct 

o 

rH 

o 

|H 

ft 

o 

s0 

o 

o 

rH 

rH 

o 

Is- 

CO 

CO 

in 

rH 

in 

CO 

Is- 

CO 

rH 

UJ 

in 

CO 

CO 

o 

1— 1 

o 

Is- 

in 

sO 

m 

CM 

in 

<*• 

CO 

vO 

CO 

o 

CO 

sO 

H 

CO 

Is- 

u_ 

• 

r-l 

co 

CM 

in 

o 

Is- 

CO 

CO 

r-l 

vt 

• 

o> 

o 

r- 

Is- 

o 

in 

Is- 

Is- 

CM 

O 

crl 

r- 1 

r-H 

UJ 

rH 

vfr 

in 

CM 

CM 

c— 1 

CO 

CO 

in 

in 

CM 

CO 

r- 

vt 

vO 

CM 

co 

in 

rH 

4- 

h- 

II 

1 

1 

1 

1 

II 

1 

1 

1 

1 

2 

^ 

i— < 

o 

r-l 

r-l 

r-l 

CM 

CM 

CM 

CO 

r-l 

t— i 

CM 

a 

rH 

rH 

o 

o 

rH 

rH 

rH 

CM 

rH 

H 

O 
1 

o 

1 

o 

1 

o 
■ 

O 
1 

O 
l 

o 

1 

o 

1 

o 

o 
■ 

O 
1 

O 
1 

o 

o 

o 

o 

O 

o 

o 

o 

1 

UJ 

1 

UJ 

1 

UJ 

i 

UJ 

l 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

i 

UJ 

i 

LU 

1 

LU 

1 

LU 

LU 

UJ 

UJ 

UJ 

LU 

LU 

UJ 

Ul 

r-t 

in 

cm 

o 

CM 

Is- 

vO 

rH 

o 

CM 

Ln 

CO 

f- 

in 

ON 

CO 

<!• 

co 

CM 

rH 

ON 

r~ 

rH 

r-l 

O 

co 

CO 

O 

CO 

in 

rH 

vO 

rH 

rH 

CO 

rH 

CO 

Is- 

»* 

«* 

co 

CM 

r- 

o 

r» 

v0 

CO 

in 

CM 

** 

CO 

rH 

CO 

CNJ 

CM 

CM 

in 

»t 

rH 

vt 

M 

in 

CO 

o 

CM 

«cf- 

CM 

CO 

O 

CO 

o 

CJ> 

co 

CO 

rH 

vj- 

in 

O 

rH 

in 

CM 

H 

>t 

in 

»}• 

n0 

co 

«*■ 

CM 

rH 

rH 

r-l 

in 

rH 

rH 

O 

rH 

in 

CO 

rH 

CM 

Is- 

> 

1 

1 

1 

1 

1 

1 

1 

1 

1 

r 

1 

W 

ii 

n 

n 

II 

ll 

11 

II 

II 

II 

II 

ii 

II 

II 

II 

II 

ii 

ii 

II 

ll 

II 

J 

o 

CQ 

• 

o 

r-l 

CM 

CO 

-* 

in 

vO 

r- 

CO 

o 

• 

o 

rH 

CM 

CO 

■4- 

in 

o 

Is- 

co 

ON 

< 
Eh 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

O 

o 

o 

o 

o 

o 

II 

>j 

_J 

_J 

_J 

_J 

_j 

_l 

_j 

_l 

_J 

n 

_j 

_J 

_l 

_l 

_l 

_i 

_l 

_l 

_J 

_l 

_J 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

_j 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

o 

O 

o 

O 

a 

Q 

o 

a 

D 

Q 

Q 

a 

o 

o 

O 

o 

a 

o 

a 

a 

Q 

Q 

74 


o 

o 

i—4 

O 

o 

rH 

O 

r-t 

o 

rH 

o 

O 

o 

r-l 

o 

rH 

o 

o 

rH 

O 

rH 

o 

rH 

in 

O 

O 

o 

o 

O 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

O 

O 

o 

o 

O 

rH 

1 

CNJ 

1 

1 

o 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

o 

UJ 

UJ 

IU 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ui 

• 

in 

h- 

in 

<f 

o 

Is- 

o> 

0> 

vO 

CJN 

• 

ro 

in 

vO 

o 

vO 

CNJ 

CO 

o 

ro 

CO 

o 

o 

C\J 

rH 

0- 

vO 

r-l 

r-t 

vO 

CO 

co 

o 

rH 

r~ 

ro 

o 

vO 

ON 

vj" 

rH 

o 

o 

co 

CO 

<t- 

ro 

rH 

vO 

in 

0> 

o 

O 

O 

in 

vO 

co 

h- 

^d• 

<r 

rH 

rH 

>d- 

II 

<s- 

ro 

CO 

r\j 

Q> 

vO 

Is- 

rH 

rH 

o 

II 

o 

o 

r- 

ON 

<J- 

in 

vj- 

rH 

ro 

ro 

_J 

_J 

< 

C\] 

co 
1 

ro 
1 

rH 

CO 

rH 

1 

co 
1 

rH 

4" 

<*- 
1 

< 

ro 

rH 

ro 
1 

ro 

rNj 

ro 

1 

1 

Is- 

CNJ 

1 

o 

o 

O 

O 

rH 

o 

r-i 

o 

O 

rH 

rH 

rH 

o 

o 

O 

o 

o 

O 

rH 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

1 

O 

r 

UJ 

O 
1 

o 

o 

O 

o 

o 

O 

o 

o 

CO 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ul 

Ul 

UJ 

1 

UJ 

1 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

1 

UJ 

UJ 

H 

cnj 

ro 

o 

r\i 

o> 

O 

in 

ro 

vO 

vO 

CO 

ro 

rjN 

Nl- 

vO 

-t 

in 

st 

rH 

o 

■z. 

Is- 

in 

ro 

vO 

in 

>t 

o 

•+ 

o 

CNJ 

CM 

H 

CNJ 

v0 

CNJ 

>*■ 

o 

O 

o 

H 

UJ 

co 

>o 

•4- 

Is- 

o 

in 

vO 

CNJ 

ON 

CNJ 

CO 

ro 

o 

o 

o 

o 

CJN 

o 

rH 

(NJ 

t— 1 

o 

vO 

m 

sf 

in 

CO 

<\i 

ro 

<r 

rH 

0^ 

>d- 

in 

CO 

<c 

o 

CO 

<*■ 

CO 

CNJ 

o 

o 

o 

t 

• 

• 

» 

« 

• 

• 

t 

• 

» 

r— < 

o 

rH 

rH 

r- 

O 

rH 

Is- 

r-l 

ro 

ro 

f- 

o 

rH 

<f 

rH 

rH 

rH 

r- 

CNJ 

Is- 

ro 

rH 

U- 

in 

1 

1 

I 

1 

1 

1 

o 

, 

1 

1 

1 

1 

1 

U_ 

• 

* 

<^ 

LU 
O 

rH 

II 

II 

ii 

II 

II 

II 

II 

ll 

ii 

II 

CNJ 

II 

ii 

II 

II 

II 

ii 

II 

n 

ii 

II 

II 

o 

rH 

CM 

CO 

«*t 

in 

vO 

l— 

CO 

fJN 

II 

o 

rH 

CNJ 

ro 

*t 

in 

vO 

Is- 

co 

0^ 

o 

rH 

■H 

>-i 

1— t 

rH 

r-l 

rH 

rH 

rH 

rH 

o 

rH 

rH 

r-l 

rH 

rH 

rH 

•H 

rH 

rH 

rH 

2 

X 

wo 

*-.•* 

«» 

wo 

wo 

wo 

wo 

wo 

wo 

wo 

I 

wo 

so 

K^ 

«~» 

>— ' 

«o> 

SkO 

wo 

wo 

wo 

o 

CC 

-J 

_J 

_i 

_J 

_J 

_J 

_J 

_J 

_l 

_J 

or 

_J 

_J 

_J 

_J 

_J 

-J 

_) 

_J 

_J 

-J 

»— 1 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

o 

O 

Q 

o 

o 

a 

Q 

O 

Q 

a 

o 

o 

D 

o 

o 

O 

Q 

Q 

o 

o 

a 

3 
u_ 

rH 

O 

o 

rH 

o 

rH 

o 

rH 

o 

o 

rH 

O 

rH 

o 

o 

rH 

O 

rH 

o 

rH 

UJ 

o 

o 

1 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 
1 

o 

O 
i 

o 

o 

O 

O 

O 

a 

O 

1 

IU 

UJ 

UJ 

UJ 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

1 

UJ 

UJ 

1 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

o 

in 

Is- 

in 

<f- 

o 

c- 

ON 

o 

o 

ON 

o 

ro 

in 

vO 

o 

o 

C\J 

CO 

o 

ro 

ro 

a: 

o 

o 

CNJ 

rH 

Sf 

vO 

r~t 

rH 

vO 

CO 

CO 

o 

rH 

Is- 

ro 

o 

o 

ON 

«t 

rH 

o 

o 

UJ 

in 

co 

CO 

<i- 

ro 

r-l 

vO 

in 

ON 

o 

vO 

in 

vO 

m 

x0 

00 

Is- 

■4- 

Nh 

rH 

rH 

<fr 

UL 

• 

«fr 

CO 

CO 

C\J 

0> 

o 

Is- 

rH 

rH 

o 

* 

o 

o 

r- 

CJN 

«st 

in 

vt 

rH 

ro 

ro 

cc 

rH 

rH 

Mi 

CNJ 

ro 

ro 

rH 

CO 

I— 1 

co 

rH 

<* 

Nt" 

ro 

rH 

ro 

CO 

CNJ 

ro 

«t 

r- 

CNJ 

«* 

\- 

II 

I 

1 

1 

1 

1 

II 

1 

1 

1 

1 

1 

1 

2 

x: 

s^ 

r— < 

o 

O 

O 

o 

O 

r-t 

o 

rH 

o 

o 

rH 

D 

rH 

rH 

o 

o 

o 

o 

o 

o 

rH 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

1 

o 

1 

o 

1 

o 

o 

o 

o 

o 

o 

O 
1 

o 

UJ 

UJ 

UJ 

IU 

UJ 

IU 

UJ 

UJ 

UJ 

1 

UJ 

1 

UJ 

1 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

1 

UJ 

UJ 

C\J 

ro 

ON 

c\J 

ON 

cr> 

in 

ro 

vO 

o 

CO 

ro 

O 

«fr 

vO 

•4- 

in 

<!- 

rH 

o 

Is- 

in 

ro 

•■o 

in 

S» 

0 

sf 

o 

CNJ 

CNJ 

rH 

(NJ 

o 

CNJ 

>t 

x0 

O 

o 

rH 

M 

co 

<J 

<T 

Is- 

vO 

in 

o 

r\j 

on 

<NJ 

co 

ro 

o 

o 

O 

o 

0s 

o 

rH 

CNJ 

> 

w 
r-3 

o 

s0 

in 

tf 

in 

CO 

cnj 

ro 

<f 

rH 

CJN 

4- 

in 

CO 

>t 

O 

co 

•fr 

co 

CNJ 

i-i 

r-l 

Is- 

vO 

r—l 

Is- 

rH 

ro 

ro 

Is- 

rH 

Nf 

rH 

rH 

rH 

Is- 

CNJ 

p~ 

ro 

rH 

1 

1 

1 

1 

1 

1 

1 

1 

CQ 

ii 

II 

ii 

II 

II 

ii 

II 

ll 

n 

ii 

II 

II 

II 

II 

II 

H 

ll 

II 

ll 

II 

< 

o 

£H 

• 

o 

rH 

CNJ 

ro 

v* 

in 

vO 

Is- 

CO 

0 

• 

o 

rH 

CNJ 

ro 

vt 

in 

o 

Is- 

CO 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

II 

—I 

_J 

_J 

_J 

_j 

_j 

_J 

_j 

_J 

_l 

II 

_J 

_J 

_J 

_i 

J 

_j 

_J 

_J 

_J 

_J 

-J 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

_l 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

D 

D 

o 

O 

Q 

o 

a 

o 

o 

D 

o 

Q 

o 

Q 

a 

a 

o 

Q 

a 

o 

D 

Q 

75 


O 

o 

t-i 

i—i 

rH 

rH 

rH 

rH 

rH 

O 

i—i 

o 

O 

r-l 

CM 

o 

o 

rH 

r-t 

rH 

rH 

rH 

rH 

in 

o 

O 

O 

O 

o 

O 

O 

O 

o 

o 

o 

o 

O 

o 

o 

O 

O 

o 

O 

o 

O 

o 

l 

1 

1 

1 

1 

1 

1 

1 

rH 

1 

1 

o 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

UJ 

UJ 

o 

LU 

LU 

Ul 

UJ 

LU 

UJ 

UJ 

LU 

LU 

LU 

• 

CO 

O 

f- 

o 

CM 

co 

v}- 

in 

CM 

O 

. 

<t- 

rH 

co 

CO 

CM 

o 

v0 

f— ! 

CO 

O 

o 

m 

in 

<t" 

CM 

vO 

•4- 

CO 

in 

CO 

CM 

o 

h- 

o 

rH 

CM 

r- 

o 

o 

CO 

P» 

<f- 

o 

00 

r- 

co 

O 

CM 

CO 

rH 

o 

O 

C0 

CO 

<*■ 

in 

cO 

r- 

o 

CO 

vO 

in 

II 

in 

CM 

r-l 

o 

CO 

h- 

CO 

o 

in 

O 

II 

o 

<f- 

rH 

in 

CO 

CO 

o 

vO 

O 

o 

_J 

t 

■ 

• 

0 

t 

» 

• 

• 

• 

0 

_J 

■ 

■ 

* 

• 

• 

* 

• 

« 

» 

• 

< 

l 

c<i 
1 

in 

co 

1 

CM 

co 
1 

1 

rH 

CO 

rH 
1 

< 

rH 
1 

r- 
1 

1 

1 

CM 

CM 

CO 

1 

CO 

1 

CM 

rH 

rH 

rH 

r-K 

CM 

CM 

CM 

rH 

rH 

T-i 

rH 

r-\ 

o 

o 

rH 

rH 

rH 

l-H 

rH 

rH 

rH 

o 

1 

O 
1 

O 

O 
1 

O 

o 

1 

o 

O 

o 

1 

O 

O 

o 

o 

O 

O 

O 

o 

O 

o 

O 

00 

LU 

LU 

1 

UJ 

LU 

1 

LU 

1 

IU 

1 

LU 

1 

UJ 

1 

LU 

1 
LU 

1 

LU 

LU 

LU 

LU 

UJ 

LU 

UJ 

UJ 

LU 

LU 

1- 

<* 

O 

CM 

O 

CT> 

o 

O 

in 

O 

O 

O* 

CO 

CM 

•tf 

rH 

CM 

CM 

o 

CM 

o 

z 

in 

m 

in 

•O 

in 

co 

o 

r- 

O 

in 

rH 

CO 

1- 

r-t 

>t 

■* 

CO 

CO 

CO 

o 

UJ 

o 

in 

co 

CO 

CM 

CM 

CM 

co 

CO 

r» 

O 

in 

00 

P- 

rH 

CM 

o 

o 

tf 

>o 

»— i 

co 

co 

C<1 

CM 

nO 

P- 

vO 

nO 

r-t 

r- 

CO 

r~t 

h- 

<T 

CO 

vO 

CO 

CO 

CO 

CM 

o 

o 

o 

►— < 

o 

P- 

Lf> 

f- 

0> 

r- 

>t 

rH 

rH 

<J' 

H 

o 

in 

rH 

rH 

r-i 

r-i 

in 

CM 

CO 

r-A 

vt 

u_ 

in 

1 

1 

1 

1 

1 

o 

1 

1 

1 

1 

1 

1 

LU- 

• 

LU 

o 

II 

ii 

ii 

11 

n 

II 

II 

II 

II 

II 

rH 

ii 

II 

II 

II 

II 

ii 

11 

II 

II 

II 

O 

— • 

»-» 

O 

ii 

o 

r-l 

CM 

CO 

<}- 

in 

vO 

r~ 

CO 

o 

II 

o 

rH 

CM 

m 

vt 

in 

vO 

p» 

CO 

o 

o 

rH 

rH 

r-4 

rH 

f-\ 

rH 

I— I 

i—t 

rH 

rH 

o 

rH 

tH 

r-A 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

2 

I 

\m* 

*-* 

>~» 

w« 

»■*# 

«■* 

W 

■w- 

«-* 

•«• 

X 

*-• 

*~* 

*-> 

•w* 

<— 

w 

N-« 

<_> 

*0> 

w 

o 

a: 

_J 

_J 

-J 

_J 

-J 

_J 

_J 

_j 

_l 

-J 

a: 

_J 

_J 

_J 

-J 

~J 

_l 

_J 

_J 

_J 

_J 

f— « 

LU 

LU 

UJ 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

uj 

UJ 

LU 

LU 

UJ 

UJ 

LU 

LU 

h- 

Q 

O 

Q 

O 

O 

O 

Q 

a 

a 

O 

O 

o 

Q 

a 

Q 

Q 

o 

o 

o 

Q 

O 

2 

3 

U_ 

LU 

O 

■z 

LU 
CC 
LU 
LL 
CC 
LU 


M 
H 
> 

W 

rH 

< 
Eh 


rHrHrHrHrHrHrHOrHO  rHCMOOrHrHrHrHrHrH 

OOOOOOOOOO  OOOOOOoOOO 

I      I      I      I      I      I      I           I  II 

LU   LU   UJ   LU   UJ   UJ   LU   LU   UJ   LU  LU   LU   LU   LU   IU   LU   LU   LU   LU   LU 

OcnOr-OcMcOvJ-incMO  o^-i-hcoodcmOnO'-'cOvO 

OinL^^tcMOsj-comcocM  or^-OrHCMr^OvOcor^-Nt 

oocoh-coOcMcOrHOO  ococo^tir^coh-ocOvOin 

•  LncMrHocof^-oooino  »Ovt--<incocoa.>ooc> 
CM        ...»        t        ....        «  (M        .........        . 

ocomcocMcoNO'-icorH  rHr*-<j-vOCMcMcococMrH 

II                 II                 I      I  II                            II                 It 

DrHrHrHCMCMCM1— (rHrHrH  QrHOOrHrHrH^HrHrHrH 

OOOOOOOOOO  OOOOOOOOOO 

I  I      I      I      I      I      I      I      I      I  I 

•       LULUUJUJLULULULULULU  IU       Ltl       LU       LU       UJ       LU        LU       UJ       UJ       UJ 

vtOCMOOOOinOO  OC0CM>4-rHCMCMOcMO 

incomOinrOOr^C^m  rHcor^rH^"<J"cocOfoO 

oincococMcMCMcocor-  omcor--rHCMOOvj-o 

cocococMor~-\OvOrHh-  corHh-^-cOvOcnc^icocM 

..........  .......... 

r-~mr-CT>r~^t'-<rHvd-rH  ir\rHrHpHrH(nCM00rH>t 

I      I           I      I                 I  II                 II 

II  II   II   II   II   II   II   II   II   II  II   II   II   II   II   II   It   II   II   II 
o~,-.  —  —  ,^_„^„„  o«  —  —  —  -»—  —  —  —  — 

•  OrHCMCO^-     m     vOf-cocj>  .OrHcMco^mvor—     coo 

OOOOOOOOOOO  OOOOOOOOOOO 

it     _j     t    — i     i     i     — i    i    — i_j     — i         ii     i     — i     i    i     i    i     — i    _i     _ i     _j 

-JLULULULUUJUJLULULULU  _ilululululululululuuj 

OOOOOOOOOOO  OOOOOOOOOOO 


76 


o 

o 

f-H 

o 

o 

o 

o 

O 

o 

r- 1 

o 

rH 

o 

r-t 

O 

o 

rH 

o 

CM 

o 

CM 

o 

rH 

in 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

o 

O 

o 

o 

O 

o 

O 

o 

o 

o 

o 

i—( 

1 

CM 

1 

o 

IU 

UJ 

IJU 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

o 

UJ 

UJ 

LU 

LU 

UJ 

LU 

UJ 

UJ 

LU 

LU 

■ 

vO 

o 

r-H 

o 

■-o 

r-H 

(NJ 

vO 

in 

o 

t 

On 

o 

ON 

m 

ro 

CM 

o 

CxJ 

vO 

ON 

o 

O 

in 

0^ 

Ni- 

CM 

CO 

vO 

nS- 

r— t 

o 

o 

<f 

c* 

CM 

o 

CO 

Is- 

o 

r- 

ro 

f- 

ro 

o 

r~ 

CO 

CO 

rH 

o 

-4- 

in 

ro 

co 

& 

rH 

ro 

o 

0N 

in 

ro 

O^ 

00 

II 

ON 

o 

o 

in 

in 

0> 

rH 

vi- 

O 

o 

II 

o 

ro 

CM 

CM 

CM 

o 

r~ 

<f- 

in 

<T 

_J 

_J 

* 

e 

• 

• 

■ 

• 

• 

• 

• 

• 

< 

0^ 

r-H 

i 

ro 
i 

in 

CO 
1 

1 

t— i 

ro 

rH 
| 

<x 

CM 
I 

ro 

ro 

ro 
■ 

o 

1 

r-l 

CM 

rH 
1 

CM 

O 

** 

*- 

■ 

»- 

1 

1 

•- 

w 

1 

l 

r- 

•« 

*• 

i 

1 

•* 

•» 

1 

r» 

•^ 

ro 

o 

o 

o 

o 

o 

o 

, j 

o 

o 

rH 

o 

a 

rH 

r-l 

>-l 

rH 

r~{ 

rH 

rH 

o 

1 

o 

o 

o 

o 

o 

o 

o 

1 

o 

o 

O 
1 

o 

o 

O 

o 

o 

O 

O 

O 

O 

(/> 

1 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

IU 

1 

UJ 

UJ 

UJ 

1 

UJ 

UJ 

LU 

LU 

LU 

UJ 

LU 

UJ 

LU 

UJ 

H- 

C\J 

ro 

co 

LOi 

r\l 

in 

o 

vt 

CM 

<f- 

rH 

ON 

rH 

o 

CM 

ON 

r-l 

rH 

CM 

ro 

Z 

rH 

f>- 

vO 

co 

in 

o 

r-- 

ro 

<J- 

rM 

rH 

in 

o 

rH 

O 

o 

CM 

<t" 

CM 

CO 

LU 

00 

f-i 

og 

rH 

m 

vO 

o 

o 

in 

CO 

rH 

o 

ON 

ro 

h- 

rM 

r- 

ON 

r> 

o 

►— i 

o> 

c- 

-4- 

o 

O 

o 

in 

<t 

0s- 

o 

O 

Is- 

Is- 

rH 

00 

o 

o 

v0 

r> 

CO 

O 

o 

o 

• 

* 

• 

* 

• 

• 

c 

e 

• 

• 

t— t 

o 

co 

rM 

ro 

vt 

<f- 

r-i 

ro 

ro 

CM 

rH 

o 

Is- 

rH 

ON 

rH 

<i- 

in 

CO 

CO 

-t 

Nf 

U_ 

in 

1 

1 

1 

1 

1 

1 

o 

1 

1 

1 

1 

LL 

* 

LU 

«-h 

II 

n 

II 

II 

li 

II 

n 

ii 

11 

II 

CM 

II 

II 

11 

II 

ii 

ii 

11 

II 

ii 

II 

O 

«-, 

»-» 

«-» 

O 

ii 

o 

rH 

rM 

ro 

<t- 

in 

vO 

Is- 

CO 

O 

II 

o 

rH 

rM 

CO 

-t 

in 

o 

r- 

CO 

ON 

o 

rH 

rH 

i— < 

r-A 

r-H 

.-h 

r-( 

T-i 

r-l 

rH 

o 

rH 

rH 

rH 

1— 1 

rH 

rH 

r-\ 

rH 

rH 

r-A 

2T 

X 

O 

or 

_l 

_l 

_i 

-J 

_J 

_j 

_J 

_J 

_J 

_J 

a: 

_1 

_J 

_J 

_J 

-J 

_J 

_J 

~J 

_J 

_J 

t—f 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

LU 

UJ 

LU 

LU 

IU 

LU 

LU 

y- 
o 
z 

o 

o 

o 

o 

O 

o 

a 

o 

o 

o 

o 

o 

O 

O 

Q 

a 

O 

O 

o 

O 

rH 

o 

o 

o 

O 

o 

o 

rH 

o 

1— 1 

rH 

o 

o 

i— I 

O 

CM 

o 

CM 

o 

rH 

UJ 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

1 

o 

o 

o 

O 

o 

o 

O 

o 

O 

1 

UJ 

UJ 

UJ 

UJ 

UJ 

IU 

UJ 

UJ 

UJ 

LU 

1 

UJ 

UJ 

LU 

UJ 

UJ 

LU 

UI 

LU 

UJ 

LU 

UJ 

o 

nO 

o 

rH 

o 

o 

rH 

CM 

o 

in 

o 

o 

ON 

vO 

ON 

in 

CO 

CM 

vO 

CM 

sO 

ON 

a: 

o 

0^ 

LH 

0^ 

sj" 

C\J 

CO 

O 

<t 

r-H 

o 

o 

<}- 

ON 

CM 

o 

co 

Is- 

vO 

r- 

ro 

h- 

UJ 

o 

ro 

O 

Is- 

00 

CO 

r-l 

o 

>t 

in 

ro 

o 

CO 

rjN 

i— 1 

CO 

sO 

0N 

m 

ro 

ON 

00 

U_ 

• 

o 

vO 

o 

in 

in 

O 

rH 

** 

v0 

o 

• 

o 

ro 

CM 

CM 

CM 

o 

Is- 

M- 

in 

<!- 

a: 

C\J 

CM 

UJ 

o 

r-( 

-J- 

ro 

in 

CO 

v£> 

T-A 

ro 

rH 

CM 

ro 

ro 

ro 

vO 

rH 

CM 

r-l 

rM 

vO 

»- 

ii 

1 

1 

1 

1 

1 

II 

1 

1 

1 

1 

1 

1 

z 

^ 

^ 

•—I 

o 

ro 

o 

o 

o 

o 

O 

O 

rH 

o 

o 

a 

rH 

o 

o 

r-A 

r-l 

r-l 

rH 

r-l 

r-l 

r-l 

o 

1 

o 

c~> 

o 

o 

o 

O 

o 

1 

o 

o 

o 

1 

o 

o 

o 

o 

o 

o 

o 

O 

o 

1 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

IU 

1 

Ui 

UJ 

U) 

1 

UJ 

UJ 

UJ 

LU 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

eg 

ro 

CO 

in 

(NJ 

in 

o 

<f- 

CM 

<t- 

rH 

ON 

rH 

vO 

CM 

0^ 

r-l 

r-l 

CM 

ro 

H 

rH 

r- 

nO 

CO 

in 

o 

Is- 

ro 

<j* 

CM 

rH 

in 

O 

rH 

O 

vO 

CM 

<t 

CM 

m 

M 

00 

rH 

(\l 

r- < 

CO 

-0 

c* 

o 

in 

CO 

rH 

o 

0> 

CO 

r- 

CM 

r- 

ON 

0s 

vO 

> 

On 

r- 

vt 

o 

o 

o 

in 

<t 

o 

o 

o 

r- 

r- 

rH 

co 

o 

o 

vO 

On 

co 

00 

(M 

ro 

<f 

<J- 

rH 

ro 

ro 

CM 

rH 

r«- 

i— i 

0N 

rH 

in 

CO 

CO 

>t 

st 

1 

1 

i 

t 

l 

1 

i 

I 

1 

l 

< 

II 

II 

II 

II 

ii 

II 

ii 

ii 

ii 

II 

n 

ii 

II 

II 

ii 

ii 

11 

ii 

II 

it 

e< 

o 

• 

o 

r— 1 

(NJ 

ro 

<!- 

in 

O 

Is- 

CO 

0^ 

• 

o 

rH 

CM 

CO 

Nj- 

in 

o 

r~ 

CO 

0> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

O 

O 

o 

o 

o 

o 

o 

II 

_J 

_j 

_J 

_j 

_l 

_j 

_i 

_i 

_l 

_J 

ii 

_J 

_J 

_J 

_J 

_l 

_J 

_J 

_j 

_i 

_J 

-J 

UJ 

UJ 

UJ 

LU 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

_i 

UJ 

UJ 

LU 

Uj 

LU 

LU 

UJ 

UJ 

UJ 

UJ 

o 

O 

O 

Q 

a 

Q 

o 

a 

o 

Q 

o 

o 

Q 

o 

o 

a 

O 

o 

o 

o 

Q 

Q 

77 


I/O 

H- 
■ZL 

UJ 

►-H 
<_) 
»— « 

u. 
u_ 

UJ 

o 
(_) 

o 
t— i 

O 


II  oooooooooo 

_j»cs«ttett* 
< OOOOOOOOOO 


OO^HH(\IMMhM 

oooooooooo 
I  I  I  I  I  I  I  I 

UJ  LU I U  U  J IJJ  UJ  UJ  UJ  LU  UJ 

rn<f'C^oo^orococovi- 
in  o  ,-i  o  ro  01s-  o>  o  o 
<j-ir\vo^-(vOcoLnh-oo 

0*t****«*«* 

O  r\|  r-4  m  r-<  i-i  >o  c\l  rO  i-i  (3^ 
O  ||  || 

rH   ||     ||     ||     II     ||     II     II     ||     ||     |! 

II  O'HCMeO^-mvOh-OOO 

CjrHrH.-4.-4i— IHHHHH 

CC"  _J  _J  _J  _  J  _J  — I  _J  _1  -J  _J 

UJUJUJUJUJUJUJUJUJUJ 
OOOQOOOOQQ 


Ncococo^a>o,o(^a> 
oooooooooo 

I  I  I  I  I  I  I  I  I  I 

O  UJ  UJ  UJ  LU  UJ  UJ  UJ  UJ  U-IUJ 

•  ^"•OvOOcO'-<fN~c\jf~o 

O  r~\  la  ro  co  1-1  c\|  .-f  vj-  rn  O 

■^■(Njr-iLrivd-vtr—  vD>-tco 

it  rOLTiLnor\;r-iLnv^i>.Ln 

„Jtc«t     (.(.«••« 

<— <lt\cm— ^-^•-(coCr^ 

II  I        III        I 

OOr-<— (r-<fNjc\Jrn'-ic\i 
OOOOOOOOOO 

I  I  I  I  I  I  I  I 

UJ  UJ  UJ  UJ  UJ  LU  U  J  UJ  U  J  UJ 

rO^tOOoOrOcocCK?- 
mOriOnCM^c>oo 
vtinvo--H<)coLr\r^-oo 
o  <\jo  lh  o  i*-  o  ro  ro  Is- 
o»«»«  **••*• 
0(NJ>-tro^-i^vOrNJm'-io 

O  II  || 

•   <•*•**  *-*  *■— *  *-*  »***  *tt*«  *■**  W*  *■*"  *— * 

--4  II     II     II     II     II     II     II     II     II     II 

ii  Or~i(Mm<j-ir>vors-coo 

D'-l<-<Hf-lHHHr-*'--IH 

a:  _j  _i  _J  _J  -  J  _j  -j  — i  _J  — i 

luujujujujujujujujuj 

doodoqqqqq 


UJ 

o 
z 

UJ 

cc 

UJ 

u_ 

UJ 


H 
H 
H 
> 

W 

< 

EH 


o 

•oooooooooo 

Ot****»«*c* 

oooooooooo 
II 

>^*f.    *-#•«-••    «*»*»-^ 

OOo>-4'-*'-<rM(\lr<)'-<r\J 
oooooooooo 
I  I  I  I  I  I  I  I 

UJ  UJ  UJ  IU  U  J  U  J  UJ  LU  UJ  UJ 

m^hOOoorncooOst 

moHOmCMsa<^o 

^■Lf\sO'-!sDCOLnix-OvO 

ocMOLnor^-ommr- 

o\j.~4crw>-i>ocMrQ!-<ai 
I   I       i       I    l  l 

0«-w«-^-w^^-w^— 

O  II  II  II  11  II  II  II  II  II  II 

CM  ^  -. ,-,  -.  ~*  „  ~,  ^,  .-,  ^ 

•  O^iCNjrOvj-Ln-oh-coo 
OOOOOOOOOOO 

II  _J— l_J_J_J_J_l_ l_l— I 

— I  UJ  UJ  UJ  LU  UJ  UJ  UJ  UJ  UJ  UJ 

OOOOOOOOOOO 


h~cocooDO^OoOa>o 
OOOOOOOOOO 

I  I  I  I  I  I  I  I  I  I 

LULU  UJ  Ul  LU  UJ  UJ  UJ  UJUJ 

ONt-sQvO>ocOf~tf--c\irs-o 
O  r-t  in  ro  co  i--i  cm  r-<  -J/  m  o 

O  v}-  CM  ft  LPi  -J/-  ^  r-  nO  r-l  CO 

*rOLnLr\or\i>-iLnvd-h-ir\ 

o«*<»***»<* 

r- iincM  — ifs->tHCO(?'>0 

II  I        I  I 

St**--*--*-* 

OOO^r- If-lCMCNJCOr-ICM 

OOOOOOOOOO 
I  I  I  I  I  I  I  I 

UJUJUJUJuJUJUJOjUJUJ 

m^-O^oO^orncocOvj- 
LnOrHoc^ic^r-ooo 

^J-LO^r-tvOCOLnh-OvO 

ocMOmor—  oroo*,is- 

CM— lrOr-tr-lxOCMC<*li— (0s 

I    I        I        I    I  I 

O w  — w  —  « —  - . —w w 

O  II    II    II    II    II    II    II    II    II    II 
(M^.-~— •—»-»— .~-.~» -— 

•  0'-trMro-4-LnvOrs-coC5>i 
OOOOOOOOOOO 

II  _J_J_J-J_J_J-J_J_J_J 

_JlUUJUJUJ|JJtUUJUJUJUJ 
OOOOOOOOOOO 


78 


V.   DISCUSSION  OF  RESULTS 

A  computer  program  was  designed  to  compute  interference 
coefficients  for  the  collocation  technique  developed  in 
Section  III-A.   For  the  steady  case,  the  resulting  coeffi- 
cients in  Table  II  were  applied  to  the  technique  developed 
in  Section  III-C  and  III-D  to  provide  the  values  of  the 
blade  surface  potentials  in  Table  III.   For  this  same  case, 
the  sonic  wind  tunnel  wall  interference  solution  of  Section 
III-E  provided  the  data  presented  in  Table  IV  which  verified 
the  collocation  technique  for  that  case. 

Although  no  "separate  source"  checks  were  made  for  the 
coefficients  presented  in  Table  VII  (The  case  of  simple 
oscillation  about  the  mid-chord  with  no  translation.),  every 
effort  was  made  to  verify  the  procedure.   Table  VIII  com- 
pares the  proven  steady  state  solution  with  a  general  small 
K  solution.   Unsteady  integral  approximation  techniques  were 
justified  in  Sections  III-B-1,  -2,  and  -3  as  summarized 
below. 

Difficulties  arising  from  the  solution  of  the  occurring 
integrals  in  the  unsteady  solution  were  not  insurmountable 
but  required  careful  handling.   Whenever  numerical  integra- 
tion of  functions  involving  ^  or  CQS  (y)  near  X  =  0  is 
attempted,  accuracy  is  suspect.   One  must  closely  inspect 
arguments,  limits,  and  results  to  ensure  validity.   For 
example,  for  the  unsteady  case  and  X  i-   0,  with  no  apparent 
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means  of  computing  the  integral  in  the  immediate  vicinity 
of  the  origin,  one  must  decide,  against  the  shortcomings 
in  computer  capacity  and  programming  methods,  just  what 
approximation  must  be  applied  and,  once  applied,  how 
accurate  it  is.   In  this  case,  the  question  arose  as  to 
how  closely  one  might  approach  the  origin  with  numerical 
integration  techniques  before  accuracy  was  degraded.   At 
the  same  time  the  required  computer  time  increased 
dramatically  as  the  integrand  became  more  ill-behaved. 

The  most  difficult  question  that  arose  in  this  investi- 
gation was  the  accuracy  of  the  correction  integral  in  the 
case  x  =  0 .   Once  this  approach  was  found  valid  for  J  =  1 
through  the  calculations  shown  in  Section  III-B-2,  there 
remained  the  need  to  improve  the  accuracy  of  the  technique 
and  ensure  its  validity  for  J  >  1.   In  question  were  the 
number  of  terms  needed  in  the  truncated  series  shown  in 
Equations  (111-40)  and  (III-41),  the  validity  of  the  inte- 
gral reduction  formula  used  in  subroutine  RED,  and  the 
selection  of  the  upper  limit  of  the  correction  integral. 

Although  four  terms  were  used  in  the  truncated  series 
for  the  case  J  =  1  (Refer  to  Appendix  A-IV) ,  only  three 
terms  affected  the  significant  figures  of  the  results  for 
the  parameters  used  in  this  investigation.   For  J  =  2,  3, 
and  4,  the  number  of  terms  used  was  decreased  by  one  for 
each  computation.   The  primary  motivation  for  this  step  was 
the  inaccuracy  of  the  integral  reduction  formula  for  large 
J  and  large  R. 
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Referring  to  Tables, VA  and  VB ,  it  was  determined  that 
numerical  inaccuracies  began  to  occur  for  large  J  in  the 
integral  reduction  routine.   The  effect  of  this  inaccuracy 
was  not  important  for  small  R  but  for  the  case  of  R  =  5, 
for  example,  it  overwhelmed  the  results.   For  integrals 
of  the  form 

/    UJ  z   ^(E/U)dU, 

as  J  becomes  large  and  for  AL  <<  1,  the  value  of  the  integral 
must  become  small.   However,  in  the  case  of  R  =  0.5,  the 
magnitude  of  the  integral  values  diminished  with  increasing 
J  up  to  J  =  9  but  increased  for  J  >  9.   In  the  case  of 
R  =  5,  the  anomaly  in  the  computation  affected  integrals 
for  J's  as  low  as  five.   (In  this  latter  case,  it  is  obvious 
that  the  erroneous  answer  for  J  -  13  must  affect  the  results.) 
The  cause  of  this  error  was  determined  to  be  differences 
taken  between  large  numbers  in  the  reduction  formula.   The 
value  obtained  for  the  Fresnel  integrals  was  only  as  accurate 
as  the  curve  fit  used  in  Subroutine  CS.   This  error  and  the 
error  generated  in  computing  the  SIN  and  COS  functions 
occurring  in  the  formula  resulted  in  erroneous  answers  of 
a  magnitude  that  was  directly  related  to  J  and  R.   To 
correct  the  problem,  the  use  of  the  integral  reduction 
routine  for  J  _>  5  was  eliminated. 

For  J  _>  5,  numerical  techniques  were  resorted  to.   Here, 
for  AL  <<  1,  the  amplitude  of  the  integrand  was  small 
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permitting  one  to  neglect  that  portion  of  the  integral  near 
zero.   This  technique  seemed  to  produce  reasonable  results, 
at  least  the  results  did  not  have  large  magnitudes  that 
would  have  affected  the  final  answers.   The  resulting 
integral  values  were  in  some  cases  less  than  the  convergence 
limits  placed  on  the  numerical  integration  subroutine,  ZSUM, 
and  were,  therefore,  assumed  to  contain  error  not  greater 
than  the  said  limits.   Further,  as  these  integral  values 
became  small,  it  should  have  been  possible  to  neglect  them 
altogether  since  the  values  of  those  integrals  with  small 
(J)  predominated  in  the  series. 

Comparative  results  for  the  case  K  =  0.5  and  R  =  0.5 
and  5.0  are  presented  before  and  after  the  error  was 
detected  in  Tables  VIA  and  VIB. 

The  upper  limit  of  the  correction  integral  was  selected 

such  that  the  approximation  cos^1^  =  C0S^DU  +  E/U^  is 
valid  in  the  range  0  <_  U  <_  AL.   The  point  AL  also  needed  to 
be  within  the  capabilities  of  the  numerical  integration 
scheme  on  the  limits  AL  <_   U  <_  X  +  1.   Such  a  point  was 
picked  at  AL  =  R/100  which  ensures  that  (For  U  <_  AL)  E/U 
was  at  least  1.0EA  larger  than  DU.   Such  an  AL  was  also 
within  the  capabilities  of  the  numerical  technique  over  the 
range  of  parameters  used  in  this  investigation.   Table  IB 
demonstrates  that  the  use  of  the  truncated  series  in  this 
approximation  decreases  its  sensitivity  to  AL. 

Computer  time  could  have  been  reduced  with  a  larger  AL 
if  one  were  willing  to  accept  reasonable  inaccuracies  as 
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the  term  DU  became  more  .significant.   A  smaller  AL  increased 
computer  time  and  the  accuracy  of  the  numerical  integration 
technique  became  more  and  more  questionable  as  the  integral 
became  more  ill-behaved.   Computer  time  for  a  single  run 
of  the  program  (for  X  =  0)  on  the  IBM  360/67  computer 
system  was  on  the  order  of  15  minutes. 

It  should  be  noted  that  in  the  absence  of  a  good  error 
measurement  yardstick,  statements  concerning  more  or  less 
error  are  indeed  relative.   It  is  difficult  to  assess 
where  computer  numerical  error  or  approximation  error'  be- 
comes more  important.   It  is  for  this  reason  that  an  inde- 
pendent check  of  these  results  against  another  source  is 
necessary.   The  Laplace  transform  solution  for  the  sonic 
unsteady  wind  tunnel  wall  interference  problem  would,  in 
theory,  be  such  a  check;  however,  this  technique  will 
result  in  integrals  of  the  same  form  as  those  used  in  the 
collocation  solution  and  will  require  similar  (or  the  same) 
solution  techniques.   It  will  be  difficult,  then,  to  obtain 
a  completely  independent  check  of  these  results. 

Once  verification  or  modification  of  results  presented 
in  this  investigation  is  in  hand,  application  of  Spreiter's 
local  linearization  technique  alluded  to  in  Section  II 
should  be  applicable  and  lead  to  reasonable  analytic 
results  for  transonic  pressure  and  moment  distribution 
prediction. 
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VI.   SUMMARY 

Following  the  basic  approach  offered  by  Gorelov  for 
subsonic  and  low  supersonic  flows,  a  collocation  technique 
was  applied  to  the  purely  transonic  case  superimposing  the 
isolated  foil  solution  provided  by  Teipel  with  unknown 
interference  potentials.   A  computer  routine  was  developed 
which  would  generate  the  coefficients  for  the  unknown 
interference  source  distributions  in  the  steady  case.   The 
interference  coefficients  in  hand,  the  potential  along  the 
blade  was  determined  through  a  second  routine.   These 
results  were  favorably  compared  with  the  potential  solution 
derived  from  a  routine  based  on  the  Laplace  transform 
treatment  of  the  sonic  wind  tunnel  wall  interference  problem 
for  the  steady  case.   All  three  of  these  computer  routines 
initially  utilized  series  expansions  which  reduced  the 
occurring  integrals  to  functions  of  the  well  known  error 
function.   Confirmation  of  these  results  led  to  an  extension 
of  the  collocation  approach  into  the  unsteady  regime. 

First  attempts  to  utilise  series  approximations  for 
the  unsteady  case  resulted  in  failure  due  to  nonconvergence 
of  the  series.   Faced  with  a  lack  of  an  apparently  simple 
solution  to  the  resulting  integrals,  numerical  integration 
was  resorted  to.   Here  difficulties  arose  centering  around 
the  pathological  behavior  of  the  integrands  near  the  origin 
as  pointed  out  in  the  Discussion  of  Results. 
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As  a  basic  numerical  integration  technique,  the  "pyramided" 
Simpson's  rule  method  was  incorporated.   For  the  steady  case 
where  the  integrals  were  easily  computed,  the  results  of 
this  method  were  checked  against  the  steady  state  case 
already  in  hand  and  found  to  agree  perfectly. 

The  greatest  difficulty  occurred  in  the  unsteady  case 
for  X  =  0  where  one  could  not  reasonably  neglect  the  inte- 
grals in  the  interval  near  zero.   To  overcome  this  problem 
an  analytical  approximation  (for  the  integral  evaluated 
near  zero)  was  introduced  and  checked  against  a  known 
analytical  solution  shown  in  Section  III.   The  technique 
involved  the  reduction  of  the  integrals  to  functions  of 
the  Presnel  integrals.   However,  this  procedure  resulted  in 
additional  problems  which  were  not  immediately  detected. 
Although  the  technique  was  theoretically  sound,  the  computer 
numerical  error  pointed  out  in  the  Discussion  of  Results 
became  apparent  after  study  of  the  output.   Once  the  cause 
was  detected,  the  anomaly  was  remedied  in  a  simple  manner. 

The  coefficients  for  the  unsteady  interference  functions 
as  presented  here  shall  require  a  "separate  source"  check. 
The  Laplace  transform  solution  for  the  oscillatory  transonic 
wind  tunnel  interference  problem  should  provide  such  a 
check.   If  these  techniques  verify  one  another,  a  third 
independent  check  could  be  the  Hamamoto  linearized  high 
frequency  solution. 

Development  of  required  computer  routines  for  the 
potential  and  pressure  coefficients  from  the  collocation 
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and  Laplace  techniques  should  be  facilitated  by  the  use  of 
subroutines  derived  for  the  unsteady  coefficient  routine 
presented  here. 

It  is  felt  that  the  collocation  technique  will  prove  a 
versatile  tool  in  transonic  flow  studies  in  that  the  free- 
dom which  one  has  in  specifying  boundary  conditions  will 
permit  treatment  of  problems  which  defy  solution  by  other 
techniques.   For  example,  the  sonic  wind  tunnel  wall  inter- 
ference technique  is  only  good  in  the  it  phase  shift  case 
while  the  collocation  technique  will  provide  results  for 
any  phase  relationship.   Further,  it  should  be  valid  for 
more  complex  boundary  conditions  such  as  those  posed  by 
wind  tunnels  with  porous  walls. 
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APPENDIX  A 
FORMULATION  OF  THE  COLLOCATION  COEFFICIENT  ROUTINE 

A-I.   MAIN  PROGRAM 

The  purpose  of  the  main  program  is  the  solution  of 
Equations  (III-16)  and  (111-17)  for  the  unknown  inter- 
ference function  coefficients. 

2^ONA#6IN+4N^/X7^EXP^(X-S)-f'Ss]dS}=0 

(III-16) 


£j?ix*$&e0K*fm\fj£tfk  exp[^(x-s)-^  sL]as}=0 


L-l^    „  „  f77      ^      X   0N 
N= 

(III-17) 

Equation  (II.I-I6)  must  be  satisfied  for  each  station  along 
the  lower  blade  while  simultaneously  equation  (111-17)  must 
be  satisfied  at  each  station  along  the  upper  blade.   For 
the  purposes  of  this  routine,  ten  stations  have  been  selected 
(X  =  -0.9+0.2r,  r  =  0,1,2, ... ,9) .   Computations  near  the 
leading  edge  are  avoided  due  to  the  singularity  there. 
For  a  given  K,  A,  and  R,  the  term 


|J?  '    S  gEXP[4-(X-S)-S  ^-oldS  =  SR(X,N)+iSM(X,N) 

2U-i     (x-s)4    Q     ^   x"s 


(A-I-l) 
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is  a  complex  function  of  the  parameters  X  and  N.   The 
complex  value  of  this  function  at  each  (X,N)  is  provided 
by  Subroutine  VISUM  and  will  be  dealt  with  under  Appendix  A-II 
WSUM  returns  the  value  of  (A-I-l)  as  real  and  imaginary 
parts  in  arrays  SR  and  SM  respectively. 

The  problem  then  is  to  evaluate  the  various  known  or 
computed  values  and  insert  them  in  proper  form  in  a  suitable 
matrix  for  simultaneous  solution. 

Observe  (111-16)  and  (111-17)  become 

9  (  . 

Z  [^0NX  +  (0in+4n)(SR(X>N)  +  isM<x*N))J  =  ° 
N=0  (A-I-2) 


and 


9  r 

I  [^INXN  +  (60N+P0N)(SR(X,N)  +  iSM(X,N))}  =  0 

(A-I-3) 

With  6TKf  and  8nTJ  determined  by  boundary  conditions  at  each 
of  the  ten  stations,  (A-I-2)  and  (A-I-3)  form  a  series  of 
twenty  equations  in  unknowns  Pq^   anci^TN>  ^  =  0,1,2, ...  ,9-. 
Note,  for  the  unsteady  case,  6M„  and /)„  are  complex. 

Taking  the  set  of  coefficients,  Q^vr,  from  (111-21)  , 
(A-I-2)  and  (A-I-3)  are  rewritten  for  the  in  phase  case  as 


N=0 


9    c 

I    [/>0NXN+/>IN(SR(X,N)  +  iSM(X,N))]    =    (SR(X,0)    -   DK   SM(X,D) 

+    i(SM(X,0)    +   DK   SR(X,1))  (A-I-4) 
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and 

I   {^INxN+^ON(SR(X»N)+iSM(X'N))j  =  (SR(X,0)-DK  SM(X,1)) 
+  i(SM(X,0)  +  DK  SR(X,1))  (A-I-5) 

where  DK  replaces  K  in  the  computer  routine. 

For  the  tt  phase  shift  case,  the  coefficients  (111-20) 
apply  and  the  lower  blade  boundary  equation  (A-I-2) 
becomes 

1    {/>oNxN+^IN(SR(X»N)  +  iSM(X»N))}  =  ~(SR(X>°)-DK  SM(X,1)) 
-  i(SM(X,0)  +  DK  SR(X,1)).  (A-I-6) 

The  system  of  twenty  equations  may  be  written  in  matrix 
form 

X  D  =  B  CA-I-7) 

y 

where  X,  D,  and  B  are  the  matrices  of  coefficients,  unknown 
/)'s,  and  right  hand  sides.  X,  D,  and  B  are  complex  and  may 
be  written 


X  =  X,  +  X2i 

5  -  *!  +  V  (A-I-8) 

and  B  =  B x  +  B2±. 
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Then  in  (A-I-7) 


(X1  +  X2i)(D1  +  D2i)  =  Bx  +  B2i 


or  (D1X1~D2^2}  +  i(D2Xl+DlX2)  =  Bl  +  B2i#     (A~I-9) 


It  may  be  concluded  that 


or 


D1X1-D2X2  =  B]_ 


D1X2+D2X1  -  B2 


{}_  -})  ( })  ■  (}) 


x2       viV      ^B2 


The  twenty  complex  equations  in  twenty  unknowns 
represented  by  (A-I-7)  have  been  transformed  into  40  real 
equations  in  ^0  real  unknowns.   One  has  only  to  solve 
for  D..  and  D^  and  set 

D  =  D  +  iD2. 

PROGRAM  STEP  ANALYSIS 

(Refer  to  Routine  MAIN  at  the  end  of  this  Appendix. ) 
Double  precision  is  used  throughout  the  routine  for 
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numerical  accuracy.   The  use  of  the  various  arrays  and 
common  storage  areas  that  are  declared  will  be  detailed 
below. 

Various  arrays  and  fixed  parameters  are  initialized 
in  MAIN  to  minimize  computational  redundancy  in  later  pro- 
gramming.  The  arrays  COC(IO),  COS(IO),  AS(10),  and  AC(10), 
zeroized  in  DO  LOOP  100,  are  utilized  in  Subroutine  WSUM. 
DL  =  A,  DK  =  K,  and  RHO  =  R,  the  controlling  parameters, 
are  initialized  and  PI  =  tt  is  set.   The  value  for 


A  =  (Mi! (A.LH) 

(DLT  +  MDK)d 


is  computed.  For  the  case  DK  =  0,  A  is  set  to  zero  to 
avoid  possible  computer  problems  involving  division  by 
zero.   Additional  parameters  are  computed  as  follov;s: 


B  -  PL (RHO)2  (A-I_12) 


C  =     D^(DK)2 .  (A-I-13) 

(DLT  +  4(DK)^ 


D  =  21DKJ (A-I-14) 


(DL)2  +  4(DK)2 


E  -   D™)2  (A-I-15) 


G  =  |  TANC2^)  (A-I-16) 
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Here    G   is    set    to 


G   =   tt/4  (A-I-l6a) 


for   the    case   DL  =    0   to   avoid   division   by    zero. 


AK   =    SIN(G)  (A-I-17) 


BK   =    COS(G)  (A-I-18) 

z     =  Rho  [(dl)2  +  Mdk)2]7'  (A-I-19) 

2    >/pT 


For  future  ease  in  changing  the  number  of  stations  on 
the  blades,  the  parameter  L  is  used.   Here  L  = 10  is  set. 
The  initial  station,  X,  =  -0.9,  and  the  interval  between 
stations,  XINT  =0.2,  are  set. 

The  first  write  statement  is  included  to  indicate  the 
controlling  parameters  DL,  DK,  and  RHO. 

Finally,  the  matrix  manipulation  section  begins. 
Referencing  (A-I-10),  array  AA(40,40)  is  dimensioned  to 
hold  the  coefficient  matrix.   BB(^O)  holds  the  matrix  of 
right  hand  sides  before  solution  and  the  matrix  of  unknowns 
after  solution.   (A-I-10)  reflects  the  diagonal  duplication 
of  the  real  and  imaginary  submatrices,  X,  and  X?  .   Note 
that  equations  (A-I-4)  and  (A-I-5)  form  adjacent  rows  in 
the  array  and  that  there  is  duplication  between  them  for  a 
given  X.   The  first  two  rows  of  X,  are  as  follows: 
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J-  A-,         •  •  •    '  A 


9 


^  SR(X1,0)   SR(X1,1)  ...  SR(X1,9) 


SR(X1,0)  SR(X1,1)  ...  SR(X199)      1         X1  X^ 

(A-I-20) 


corresponding  to  the  unknowns 


^ROO      HR01    •*•  ^R09       ^RIO       ^Rll   *•'  " 


R19 


where  7>m   =  PRMN  +  ^imn*   The  first  two  rows  of  x2  are : 

0      o     ...    o     si(x1,o)  si(x1,i)  ...  si(x1,9) 

SI(X]50)  si(x1,i)  ...  SI(X  ,9)      0         0    ...    0 

(A-I-21) 

corresponding  to 
^100      ^101    '••  ^109       "ill       ^112   **•  ^119' 

Reference  to  (A-I-10)  shows  the  placement  of  these  rows. 
The  first  row  of  X   fills  AA(1,1)  to  AA(1,20)  as  well  as 
AA(21,21)  to  AA(21,*»0).   The  first  row  of  X"2  fills  AA(1,21) 
to  AA(1,*J0),  and,  with  signs  reversed,  AA(21,1)  to  AA(21,20). 
For  the  it  phase  shift  case,  the  right  hand  side  matrix  is 
filled  four  rows  at  a  time  as  follows: 
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BB(1)  =  -(SR(X1>0)  -  DK  SI(X  ,1)) 

BB(2)  =  (SR(Xl90)  -  DK  SI(X  ,1)) 

BB(21)  =  -(S1(X1,0)  +  DK  SR(X  ,1)) 

BB(22)  =  (SI(X1,0)  +  DK  SR(X  ,1)) 

and  (A-I-22) 

BB(3)  =  -(SR(X2,0)  -  DK  SI(X2,1)) 

BB(4)  =  (SR(X2,0)  -  DK  SI(X2,1)) 

BB(23)  =  -(SI(X2,0)  +  DK  SR(X2,1)) 

BB(24)  =  (SI(X„,0)  +  DK  SR(X0,1)) 


and  so  forth. 

DO  LOOP  1000  on  J  =  1,1/ varies  the  station  from  X.  to 
X-  Q.   It  sets  two  rows  in  submatrices  X-,  and  X  ? ,  or, 
equivalent ly,  four  rows  in  AA  and  four  in  BB.   A  is  used 
to  set  the  power  of  X  in  submatrix  X,  and  is  initially  set 
to  A  =  X   =  1.   JJ  is  a  counter  used  to  select  every  other 
row,  i.e.,  JJ  =  1,3. ...,19  as  J  =  1,2,. ..,10. 

With  X   set  at  X  =  -0.9,  V/SUM  is  called  to  provide  the 
integral  term  in  the  boundary  equations.   WSUM  computes 
the  integral  for  X  and  for  N  =  0,1,2,. ..,9.   These  ten 
values  are  returned  in  common  arrays  SR  and  SM.   SR  holds 
the  real  parts  and  SM  the  imaginary  as  shown  in  functional 
notation  in  (A-I-l).   The  arrays  are  filled  as  follows: 

9^ 


SR(1)   =   SR(X,0)       SM(1)   =   SM(X,0) 
SR(2)   =   SR(X,1)       SM(2)   =   SM(X,1) 


(A-l-23) 


SR(10)  =   SR(X,9)       SM(IO)  =   SM(X,9) 

DO  LOOP  900  on  I  =  1,10  controls  the  columns.   The 
columns  of  submatrices  X,  and  X~  are  filled  in  sequence  for 
each  station  counted  by  DO  LOOP  1000.   In  each  step  in 
LOOP  900,  four  positons  in  submatrix  X,  are  filled  with  the 
power  of  X.   Four  positions  in  X?  are  zeroized.   The  real 
part  of  the  applicable  integral  is  placed  in  four  positions 
in  X,  via  the  dummy  variable  SUMR.   Similarly,  the  imaginary 
part  is  placed  into  Xp  via  SUMI.   Care  is  taken  to  obtain 
the  proper  signs.   Once  the  LHS  is  set  for  each  station, 
four  rows  in  the  RHS  matrix  BB  are  set  according  to  (A-I-22). 
The  station  is  then  augmented  by  the  interval,  XINT,  and 
the  procedure  is  repeated  until  all  positions  in  AA  and  BB 
are  filled. 

For  verification,  the  matrix  prior  to  solution  is  printed 
out.   Matrix  BB  is  destroyed  in  the  solution  process. 

The  matrix  is  solved  for  the  unknowns  using  library 
Subroutine  DGELG.   The  inputs  are  the  coefficient  matrix, 
RHS  matrix,  coefficient  matrix  size  (here  set  to  40  by  M=^L), 
RHS  matrix  size,  and  a  significance  measure  (here  set  to  1.0E-6) 
The  solution  is  returned  on  BB  and  any  error  indication  on  IER. 
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BB(1)  to  BB(10)  contain  the  real  parts  of  the  coefficients 

and  BB(ll)  to  BB(20)  contain  the  imaginary  parts.   The  parts 

are  combined  in  complex  form  and  placed  in  complex  array  ZB. 
The  form  of  ZB  is 

ZB(1)   -   DELT(0,0)   = 
ZB(2)   =   DELT(0,1)   = 


ZB(ll)  =   DELT(1,0)   =  (A-I-24) 


ZB(20)  =   DELT(1,9)   = 

Finally  the  desired  coefficients  are  printed  out  on 
computer  paper  output  and  on  cards. 
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A-II.   SUBROUTINE  WSUM 

The  purpose  of  WSUM  is  to  provide  the  control  for  the 
computation  of  the  integral 

r-    X  SN      K2      OR2 

e   2  l/7r  -1  Cx-s)%    Q     *  x"s    ik_llmml) 

As  is  demonstrated  in  Section  III-B,  (A-II-1)  may  be 
reduced  to 

Xe  "   I  /I"  {hn  -   "Nil  (A-"-2> 


where 


I         -    ^(-Q^W^1       rX+1fuJ-/EXP[AAU-B/U] 
-"NR         Tin(J-l)!(N-J+l)!~   _j  I 

J_1  °  (H-30) 


•   COS  [2KAU+E/U]]  dU 


and 


41/.UJ_1KIIVK-J+1  IHf    ,     f 

fJj-lJKN^Dl       „'         (UJ"?EXP[XAU-B/U] 


SIN  [2KAU+E/U]  ]  dU 


where  the  fixed  parameters  A,B,C,D,  and  E  are  set  in  MAIN 
[Refer  to  Eqns.  (A-I-ll,  -12,  -13,  -14,  and  -15).]  and 
carried  to  WSUM  in  common  space  HOLD. 
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Given  the  integral  values,  the  summation  is  a  simple 
computer  manipulation.   The  solution  of  these  integrals 
involves  three  variations  of  the  same  basic  numerical  inte- 
gration technique.   The  numerical  integration  of  (111-30) 
and  (III-31)  is  accomplished  by  a  "pyramided"  Simpson's 
rule  using  end  correction.   This  technique  is  presented  in 
Appendix  A-III  where  Subroutine  ZSUM  is  described.   ZSUM 
prepares  the  numerical  solution  of  the  basic  integral  for 
whatever  values  of  J,K,A,  and  upper  and  lower  limits,  XO 
and  XA,  are  specified. 

The  difficulty  encountered  with  the  numerical  integration 
of  this  basic  integral  occurs  as  the  result  of  the  lower 
limit,  XA  =  0.   For  J  <  2.5,  the  amplitude  of  the  integrand 
may  be  unbounded  near  XA  =  0 .   Further,  except  in  the  special 
case  where  K  =  0,  the  SIN  and  COS  functions  in  the  integrand 
become  unmanageable  near  zero  where  the  controlling  term 
in  their  arguments,  E/U,  becomes  large.   The  resulting 
infinitesimal  period  cannot  be  treated  by  numerical  tech- 
niques.  For  a  reasonable   X,  one  may  neglect  the  troublesome 
portion  of  the  integral  near  zero.   This  is  permissable 
since  the  exponential  function  of  (-B/U)  forces  the  ampli- 
tude  to  zero  in  the  interval  where  the  integrand  is  ill- 
behaved.   When  X   =  0,  for  J  <  2.5,  one  may  not  neglect  the 
integral  near  zero.   Fortunately,  in  this  case,  one  may 
approximate  the  ill-behaved  portion  of  the  integral  by 
analytical  methods. 
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For  A  =  0,  the  basic  integrals  in  Eqns.  (111-30)  and 
(111-31)  become 


X+1  J-f. 


SUMR  =   /    U   *C0S[DU  +  E/U]dU      (111-33) 
0 


and 

SUMI  =   /    U   £SIN[DU  +  E/U]dU      (III-31) 
0 

Introducing  the  correction  term  of  Section  III-B,  one  may 
write  (111-33)  and  (111-3*0  as 


X+1  ,  sr 
SUMR  =  COC  +    .         U  *C0S[DU  +  E/U]dU    (A-II-3) 
AL; 


and 


X+1  T  r 
SUMI  =  COS  +    /    U  *SIN[DU  +  E/U]dU    (A-II-4) 
AL 


where  the  correction  terms  are 


AL  ,  r 


COC  =   /   UJ  *C0S[DU  +  E/U]dU        (A-II-5) 


and 


AL  T  r 
COS  =   /   UJ~2SIN[DU  +  E/U]dU        (A-II-6) 
0 

Integrals  COC  and  COS  are  prepared  by  Subroutine  COR. 
The  technique  as  presented  in  Appendix  A-V  provides  an 
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approximation  which  reduces  the  correction  integral  to  the 
well  known  Fresnel  integrals. 

The  limit  AL  is  selected  such  that  (1)  the  numerical 
integration  technique  may  be  applied  to  the  second  term 
in  (A-II-3)  and  (A-II-4),  and  (2)  the  correction  approximation 
is  valid  over  its  interval. 

STEP  ANALYSIS 

(Refer  to  the  Subroutine  at  the  end  of  this  Appendix.) 

Subroutine  WSUM  is  called  to  provide  the  block  of  L 
integrals  at  each  station  on  the  blades  where  the  boundary 
conditions  are  enforced.   (In  this  investigation  L  is  set 
at  ten.)   The  subroutine  inputs  are  DK  =  K,  DL  =  X}    RHO  =  R, 
and  X. 

In  the  first  steps  several  fixed  parameters  are  estab- 
lished for  convenience.   EPS  =  1.0E-8  is  the  convergence 
parameter  used  in  the  numerical  integration  process  in 
Subroutine  ZSUM.   LIM  =  50  is  a  step  limiter  also  used  by 
ZSUM  [Refer  to  A-III].   PI  =  it  is  used  in  later  calculations 
in  WSUM. 

A  test  Is  made  to  determine  which  of  the  three  cases 
is  being  dealt  with.   If  K  =  0,  branching  is  executed  to 
Step  1.   If  not,  the  next  step  is  encountered.   If  K  ?*  0 
and  X   =  0,  the  routine  branches  to  step  2.   The  case  K  ^  0 
and  A  j-   0  follows  at  step  3. 

Step  1,  K  =  0 

The  intention  here  is  to  use  the  numerical  integration 
process  over  the  entire  interval.   Since  no  SIN  or  COS 

10H 


functions  enter,  one  may,  approach  closely  to  zero  at  the 
lower  limit  (XA  is  set  to  1.0E-8).   The  upper  limit  is  set 
to  XO  =  X  +  1.   Subroutine  ZSUM  is  called  to  compute  the 
integrals  J  =  1,2,..., 10.   The  Integra]  in  (11-30)  is 
returned  on  SUMR  and  the  negative  of  the  integral  in  (11-31) 
is  returned  on  SUMI.   (Note  that  in  this  case  SUMI  =  0.) 
The  values  are  stored  in  the  vector  arrays  AC  and  AS.   The 
routine  then  jumps  to  step  1000  where  the  summation 
indicated  in  Eqns.  (11-30)  and  (11-31)  takes  place. 

Step  2,  X  =  0,  K  j-    0 

In  this  step  the  special  case  of  A  =  0  is  computed. 
The  Presnel  integral  is  obtained,  and  where  applicable,  the 
correction  is  computed  and  summed  with  that  portion  of  the 
integral  obtained  by  numerical  methods. 

Observing  that  the  value  of  X  does  not  affect  the 
correction  integrals,  they  need  be  computed  only  once. 
Hence  a  test  is  made: 

IF(I.GT.l)  GO  TO  22. 
Step  22  bypasses  the  correction  computation. 

The  first  steps  in  computing  the  correction  integrals 
includes  the  preparation  of  inputs  for  Subroutine  COR  and 
Subroutine  Red.   The  rationale  for  these  functions  is  given 
in  A-IV  and  A-V.   They  include  the  following: 

AL   =   RHO/100  (A-II-7) 

AR   =   E/AL  (A-II-8) 
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SS   =   SIN(E/AL)  (A-II-9) 

CC   =   COS(E/AL)  (A-II-10) 

The  Fresnel  integrals  for  upper  limit  AR  are  obtained 
from  library  Subroutine  CS.   The  explanation  for  CS  is 
given  on  call  from  the  library  and  a  copy  is  found  at  the 
end  of  this  Appendix.   The  only  argument  required  is  AR. 
The  following  is  returned. 


AR 

CCC      =fCOSW    dT  (A-II-lla) 

J  /2ttT 


o 

AR 


=  C  SINI 
J  /2ir5 


SSS      =1    "ni?)    dT.  (A-II-llb) 


In  preparation  for  their  use  in  Subroutines  COR  and  RED, 
(A-II-lla, b)  are  modified  as  follows: 


CCC  =  j/5  [i  _  CCC]  (A-II-12a) 


and 


SSS  =  J^r   [~  -  SSS].      -  (A-II-12b) 


Refer  to  A-V  for  an  explanation  of  these  functions. 

The  various  inputs  in  hand,  Subroutine  COR  is  called  to 
prepare  the  block  of  L  =  10  correction  integrals  for  Eqns . 
(A-II-5,  -6).   The  arguments  of  COR  are  AL,  D,  E,  CC,  SS , 
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CCC,  and  SSS.   The  value  of  (A-II-5,-6)  for  each  J  is 
returned  on  CCOR  and  SCOR  respectively  and  stored  in  arrays 
COC  and  COS.  COC  and  COS  are  of  vector  size  ten  and  are  held 
in  common  with  MAIN  for  permanency.   Also  note  that  (-SCOR) 
is  stored  in  COR  for  compatibility  with  the  numerical  solu- 
tion of  the  remaining  integral.   The  correction  terms  are 
printed  for  reference. 

Finally  the  remainder  of  the  integrals  are  computed 
numerically  in  the  same  manner  as  was  used  in  step  1  but 
with  lower  limit  AL.   The  result  is  summed  with  its  correction 
and  stored  in  AC  and  AS. 

Step  3.   K  y*  0,  A  ?    0 

In  this  step,  one  relies  on  the  exponential  function  to 
reduce  the  amplitude  of  the  integrand  near  zero.   That  portion 
of  the  integral  from  zero  to  1.0E-3  is  neglected.   The 
numerical  solution  is  called  as  before  and  the  results 
stored  in  AC  and  AS. 

Once  the  values  of  the  integrals  for  J  =  1,2,..., 10  are 
in  hand,  the  indicated  summations,  Eqns.  (111-30,-31) s 
must  be  carried  out.   Step  1000  begins  this  operation.   The 
summation  is  carried  out  in  a  straight  forward  manner  for 
each  of  the  L  =  10  integrals.   I„R  and  IWT  are  computed  on 
SUMC  and  SUMS  respectively.   (Note  that  IFAC(M)  is  a  function 
subroutine  generating  the  factorial  of  M.   Refer  to  A-VI.) 
The  completed  integrals  are  stored  in  arrays  SR  and  SM  and 
are  printed  for  reference. 
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Recalling  that  equation  (A-II-2)  contained  Q  =  A+2iK, 
a  final  complex  manipulation  must  be  made.   Observe  that 


where 


Then 


where 


and 


for 


Q 


VA2  +  4K2  [COS(to)  +  iSIN(w)]     (A-II-13) 


co  =  ARCTAN(2K/X). 


|^=  Z[BK  +  iAK]  (A-II-14) 


R  [A2  +  4K2]"* 


(A-I-19) 


2/P 
AK  =  SIN(G)  (A-I-17) 


BK  =  COS(G)  (A-I-18) 


G  =  |  ARCTAN(2K/A).  (A-I-16) 


Finally,  (A-II-2)  may  be  written 


IB  =  Z(INRBK  -  INIAK)  +  iZClNRAK  .+  INIBK).     (A-II-15) 
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Z,  AK,  and  BK  are  computed  in  MAIN  to  eliminate 
repetition  and  are  carried  to  WSUM  in  common.   DO  LOOP  1030 
computes  the  real  and  imaginary  parts  of  (A-II-15)  for  the 
ten  integrals  and  returns  these  values  to  MAIN  in  arrays 
SR  and  SM. 
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A-III.   SUBROUTINE  ZSUM 

Subroutine  ZSUM  computes  the  value  of  the  integral  in 
the  general  Eqns.  (111-30)  and  (III-3D  by  numerical 
methods . 


XO  T_r 
SUMR  =    /   U  2"EXP[CU-B/U]C0S[DU+E/U]dU   (A-III-1) 
XA 


XO    5- 
SUMI      /   U  *EXP[CU-B/U]SIN[DU+E/U]dU   (A-III-2) 
XA 


The  technique  utilized  is  Simpson's  rule  with  end  correction. 

Noting  the  behavior  of  the  integrand,  especially  near 
XA  =  0,  it  is  not  satisfactory  to  compute  the  integral  using 
a  fixed  interval.   Instead  a  "pyramided"  scheme  is  used 
which  automatically  reduces  the  size  of  the  subintervals 
until  the  scheme  can  compute  the  integral  to  within  the 
limits  of  a  satisfactory  convergence  criterion. 

The  general  form  of  Simpson's  rule  with  end  correction 
is 


SUM  =   ~(7Y1+16Y2  +  14Y3+16Y1{+.  ..+16YN-1+7YN) 

w  •  i 

+   —  (Y      -Y      ) 

or  N^l  N-3 

SUM   =   ^(TCY^Y^  +  ieJY^m^^^^CY^-Y^) 

(A-III-3) 
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where  Y±   and  Y±      are  the  value  of  the  integrand  and  the 
derivative  of  the  integrand  evaluated  at  station  (i),  and 
H  is  the  length  of  the  subinterval  [Refer  to  Fig.  8]. 


?*-x 


Figure  8.   INTEGRATION  DIAGRAM 


Figure  9  demonstrates  the  "pyramid"  algorithm.   X,  and 
Xj.  are  the  first  and  last  station  of  the  interval  over  which 
Y(X)  is  to  be  integrated,  i.e.,  X-,  =  XA  and  XN  =  XO.   N  is 
the  number  of  stations  at  which  the  integrand  must  be 
evaluated.   There  are  (N-l)  equal  sub-intervals  of  length  H 
in  INT  =  XN  -  X1. 


Step  N   H 


coefficients  of  X 


"■X 


N 


1 

3 

INT/2 

7                                               16                                             7 

2 

5 

INT/4 

7                      16                      14                      16                    7 

3 

9 

INT/8 

7          16          14          16          14          16          14          16        7 

4 

17 

INT/16 

7    16    14    16   14    16    14    16    14    16   14    16    14    16   14  16  7 

PYRAMID  SCHEME 
Figure  9 
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The  accuracy  of  Simpson's  approximation  increases  with 
decreased  H.   The  pyramided  technique  progressively  computes 
the  integral,  doubling  the  number  of  sub-intervals  each  time, 
until  the  value  of  the  computation  ceases  to  change  to 
within  some  desired  criterion,  EPS. 

Figure  9  shows  the  change  in  N  and  H  with  each  step. 
The  coefficients  for  each  step  are  shown  in  their  relative 
positions  according  to  station.   Observe  that  once  a  Y.  is 
evaluated  for  a  given  step,  it  is  used  in  each  succeeding 
step  and  need  not  be  reevaluated.   Further  note  that  the 

first  time  a  Y.  is  used,  its  coefficient  is  16;  thereafter 

1        '  ' 

it  is  14. 

The  integrand  in  (A--III-1)  is 


YD(U)   =   UJ  *EXP(CU-B/U)COS(DU+E/U) ,     (A-III-4) 
n 


and  it  may  be  shown  that 


Y  *  (U)  =UJ"2EXP(CU-B/U)[cOS(DU+E/U)(J-5/2+CU+B/U) 

(A-III-5) 
-(DU~E/U)SIN(DU+E/U)?  . 


The  integrand  in  (A-III-2)  is 


t  I 

Y  (U)   =   U  *EXP(CU-B/U)SIN(DU+E/U)      (A-III-6) 


and 
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Yj'CU)   =   UJ  JEXP(CU-B/U)[siN(DU+E/U)(J-5/2+CU+B/U) 

(A-III-7) 

+  (DU~E/U)C0S(DU+E/U)7  . 

STEP  ANALYSIS 

(Refer  to  the  subroutine  at  the  end  of  this  section.) 

Subroutine  ZSUM  arguments  are  as  follows: 

XA  =  Lower  limit 

XO  =  Upper  limit 

J   =  Summation  index  set  in  WSUM. [Refer  to  Eqns . 
(11-30,-31).] 

The  following  parameters  are  brought  to  ZSUM  via  common 
space  HOLD:  B,C,D,E,EPS,  and  LIM.   B,C,D,  and  E  are  computed 
in  MAIN  [Refer  to  Eqns.  (A-I-12 ,-13,-iH ,  and  -15).].   EPS 
and  LIM  are  set  in  WSUM.   XA,  XO,  and  J  are  also  set  in  WSUM. 
The  values  of  (A-III-1,  -2)  are  returned  on  SUMR  and  SUMI . 

Equations  (A-III-*!  ,~5  ,-6,  and  -7)  are  first  evaluated 
at  the  left  endpoint.   AA  =  CU-B/U  is  the  argument  of  the 
exponential  function.   If  AA  is  less  than  (-172),  it  is 
assumed  that  the  exponential  term  drives  the  integrand  and 
its  derivative  to  zero  by  branching  to  statement  6  (Computer 
underflow  errors  are  generated  for  smaller  values  of  the 
exponent.).   BB  =  DU+E/U,  the  argument  of  the  SIN  and  COS 
functions,  and  CC  =  DU-E/U,  a  term  used  in  Y.  ,  are  computed. 
The  evaluations  of  YR(XA)  and  Y-^XA)  are  obtained  on  YAR 
and  YAI  from  the  arithmetic  statement  functions  YR  and  YI 
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where  YR  is  (A-III-*0  and  YI  is  (A-III-6).   EE  and  FP  are 
additional  terms  used  in  the  derivative.   Finally  the  values 
of  YR  (XA)  and  Yj  (XA)  are  computed  on  YARP  and  YAIP. 

Statement  6  sets  all  values  YR(XA),  Y-^XA),  YR'(XA), 

t 

and  Y_  (XA)  to  zero  for  small  values  of  the  exponential 

function  (those  cases  where  AA  <_  -172). 

The  same  process  is  executed  for  the  right  end  point. 
Here  YR(XO)  =  YOR,  Y].(XO)  =  YOI,  YR' (XO)  =  YORP,  and 
Yj' (XO)  =  YOIP. 

The  first  step  in  computing  the  internal  points  is  to 
set  the  interval  of  integration  in  DINT  =  XO  -  XA .  Next, 
various  summation  parameters  are  zeroed  including  SOR  and 

SOI.  the  sums  of  the  Y„ . ' s  and  YT.'s  which  will  be  multiplied 

5  Ri        Ii  * 

by  the  coefficient  Ik    [Refer  to  Eqn .  (A-III-3)]  and  SER 
and  SEI,  those  sums  to  be  multiplied  by  16.   SP  sets  the 
number  of  subintervals  in  each  step  (initially  set  at  2) 
and  K  sets  the  number  of  terms  in  SER  and  SEI  to  be  computed 
in  each  step  (initially  set  at  1). 

DO  LOOP  100  controls  the  "pyramid"  level  or  step  from 
one  to  LIM.   H  is  the  length  of  each  subinterval  in  each 
step,  H  =  DINT/SP. 

On  each  step,  the  Y.'s  computed  the  previous  step  for 
coefficient  16  (summed  in  SER  and  SEI)  are  added  to  those 
with  coefficient  lh    (SOR  and  SOI).   SER  and  SEI  are  then 
zeroed  in  preparation  for  the  next  computation.   XE  is  the 
dummy  variable  in  this  portion  of  the  calculation.   The 
initial  Y.  is  computed  at  XE  =  XA  +  H. 
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DO  LOOP  10  computes,  the  K  new  Y. ' s  for  SER  and  SEI. 
AA  and  BB  are  computed  as  before.   Y„(XE)  and  YT(XE)  are 
computed  on  YER  and  YEI  and  summed  on  SER  and  SEI  respectively 
Finally  the  dummy  variable  is  augmented  by  2H  and  the  loop 
continues . 

SUMR  and  SUMI  hold  the  values  of  (A-III-1)  and  (A-III-2) 
for  any  given  step.   TEMPR  and  TEMPI  are  used  as  temporary 
storage  for  the  values  in  SUMR  and  SUMI  computed  in  the 
previous  step.   Once  a  new  set  of  values  are  computed  in 
SUMR  and  SUMI,  they  are  compared  with  the  previous  values. 
If  both  do  not  differ  from  the  previous  values  by  more  than 
EPS,  convergence  criteria  is  met  and  the  loop  terminates. 
Should  the  process  reach  LIM  Iterations  without  convergence, 
a  warning  is  printed  and  the  existing  values  of  the  integrals 
are  accepted. 
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A-IV.   SUBROUTINE  COR 

Subroutine  COR  controls  the  computation  of  the  correc- 
tion integrals  (A-II-5)  and  (A-II-6)  for  the  case  X   =  0. 
The  integrals  are 


AL  T_£ 
CCOR      /  U  *COS(DU+E/U)dU        (A-IV-1) 
0 


and 


AL  J-*" 

SCOR       /   U  «SIN(DU+E/U)dU        (A-IV-2) 

0 


where  D  and  E  are  specified  in  (A-I-14 ,-15) .   An  analytic 
approximation  for  these  integrals  exists  in  theory;  however, 
computer  numerical  error  prohibits  reliable  use  of  this 
technique  for  J  >_  5  and  R  >  2  [Refer  to  the  Discussion  of 
Results  for  an  explanation].   For  this  reason,  the  approxi- 
mation is  used  only  for  J  <_  4 .   This  range  includes  those 
troublesome  integrals  having  integrands  with  infinite  ampli- 
tude near  the  origin.   Standard  numerical  integration  is 
used  for  J  >  4  where  the  amplitude  of  the  integrands  becomes 
small  near  zero  thus  enabling  one  to  approach  zero  sufficiently 
close  to  ensure  that  the  omitted  interval  is  negligable. 

CASE  FOR  J  <  4 

Under  the  assumption  that,  for  a  properly  selected  AL, 
as  U  approaches  zero  the  term  E/U  dominates  DU  in  the 
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arguments  of  SIM  and  COS  then  equations  (111-40)  and 
(III-41)  are  reasonable  approximations  of  (A-IV-1)  and 
(A-IV-2)  even  v;hen  a  reduced  number  of  terms  is  included. 
Such  a  selection  of  AL  is  AL  =  R/100  which  ensures 
(E/U)/DU  >  1.0E4  for  U  <  AL. 

Subroutine  COR  carries  out  the  operation  of  (111-40) 
and  (TII-41)  where 

AL  T_£  AL  t_3 

CCOR  =   /  U°   *COS(E/U)dU  -  D  /   U  *SIN(E/U)dU 

0  0 

D2   AL   j   i  n3   AL  T+1 

~-     f      U  *COS(E/U)dU  +  Z-     f      IT  *SIN(E/U)dU 

*     0  °  0 


(A-IV-3) 


and 


AL  T  5  AL  T_3 

SCOR  =   /   IT  *SIN(E/U)dU  +  D   /   IT  *  COS  (E/U  )dU 
0  0 

D2    AL  T-±  D3    AL  J+-i 

-~     f      U     *SIN(E/U)dU  -  t-  /   IT  *COS(E/U)dU. 

^  0  °  0 

(A-IV-4) 

Due  to  the  inaccuracy  in  the  computation  of  the  integrals 
by  Subroutine  RED,  a  power  of  U  greater  than  1.5  is  avoided 
by  reducing  the  number  of  terms  in  the  truncated  series  by 
one  each  time  the  series  is  used.   That  is,  for  J  =  1,  four 
terms  are  used;  for  J  =  2,  three  terms  are  used;  and  so 
forth  until  only  one  term  is  used  for  J  =  4. 
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CASE  FOR  J    >    H 

For  this  case  Subroutine  ZSUM  is  used  to  numerically 
compute  the  integral  over  the  interval  XA  to  AL  where  XA 
is  picked  as  close  to  the  origin  as  computer  time  and 
accuracy  warrant  (XA  =  1.0E-4  was  seclected  for  these 
results).   It  should  be  noted,  however,  that  in  some  cases 
the  value  of  the  integral  over  this  interval  may  be  less 
than  the  convergence  limit,  EPS,  used  in  ZSUM  in  which  case 
the  results  may  suffer  from  inaccuracy.   Provided  the  error 
is  on  the  order  of  EPS  (EFS  should  be  the  upper  limit)  such 
inaccuracies  are  not  significant  and  should  not  affect  the 
end  result . 

STEP  ANALYSIS 

(Refer  to  the  subroutine  at  the  end  of  this  Appendix.) 
The  inputs  to  Subroutine  COR  include: 
AL  =   Upper  limit  of  the  correction  integral 
J   =   Summation  index  in  equation  (111-30)  and  (111-31) 
3 


D 


2K" 


A   +  4K 


2 

p   _  kr" 

E    -  — 

C    =   COS(E/AL) 
S    =   SIN(E/AL) 

AL, 

SS   =    S      ~   SIN(E/U)dU 
0    U* 

AL  , 
CC   =    /   -~  COS(E/U)dU 

0   u3/* 
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D  and  E  are  fixed  parameters  computed  in  MAIN;  C,  S, 
SS,  and  CC  are  formed  in  Subroutine  WSUM.  CCOR  and  SCOR 
return  the  values  of  (A-IV-3)  and  (A-IV-4)  to  WSUM. 

The  first  step  is  to  determine  if  J  =  1,  2  <  J  <  i|, 
or  J  >  4.   For  the  case  J  =  1,  four  integrals  are  computed 
in  (A-IV-3, -4 )  of  the  form 


AL  M-* 

CINT  -   /   IT1  *COS(E/U)dU  (A-IV-5) 

0 


and 


AL  M-f 

SINT  -   /   UH  *SIN(E/U)dU  (A-IV-6) 

0 


where  M  =  J,  J+l,  J+2,  and  J+3.   These  integrals  are  easily 
computed  in  four  step  DO  LOOP  10  utilizing  Subroutine  RED. 
RED  reduces  (A-IV-5, -6)  to  the  general  form 


Integral  =  f(AL,N,E)  +   /   -~  ^(E/U)dU .     (A-IV-7) 

0     u%^o 


The  values  of  these  integrals  are  stored  in  two  arrays, 
SI  and  CI,  of  dimension  four.   Note  that  should  a  full  four 
term  truncated  series  be  desired  for  J  >  1,  only  one  addi- 
tional integral  is  needed  in  (A-IV-3, -4)  each  time.   That 
is,  in  terms  of  the  J  of  the  preceding  calculation,  values 
of  CINT  and  SINT  are  required  with  M  =  J+l,  J+2,  J+3,  and 
J+4.   Since  the  first  three  are  already  on  hand,  DO  LOOP 
13  sets  SI(I+1)  into  SI(I)  and  CI(I+1)  into  CI(I).   The 
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fourth  position  is  filled  with  the  newly  calculated  values 
of  CINT  and  SINT  for  M  =  J+4.   However,  assuming  error  in 
RED  for  J  >  ^1 ,  it  was  elected  to  truncate  the  series  by 
one  additional  term  for  each  case,  J  =  2,3,  and  4 .   Hence, 
zero  is  set  into  SI (4)  and  01(4). 

One  assumes  that  this  reduction  in  accuracy  is  still 
acceptable  compared  to  the  error  introduced  by  an  attempt 
to  use  numerical  integration. 

In  the  case  J  >  H ,  a  branch  is  provided  to  statement  21 
where  subroutine  ZSUM  is  used  to  calculate  the  integral 
between  AL  and  XA.   Although  this  procedure  simplifies 
the  handling  of  COR  results  by  WSUM,  inaccuracies  of  the 
order  of  EPS  are  expected  as  previously  mentioned.   One 
might  expect  greater  accuracy  using  another  technique 
which  did  not  compute  the  integral  piecemeal. 
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A-V.   SUBROUTINE  RED 


Subroutine  RED  provides  the  value  of  the  general  integrals 


A  t-£ 
SINT  =   /  \J°    *SIN(E/U)dU  (A-V-l) 


and 


A  J-* 

CINT  =   /  \J°    *COS(E/U')dU  (A-V-2) 


by  reducing  them  to  functions  of  the  form 


A  n 


SINT  =  P.  (A,J,E)  +   /  \   ^™(E/U)dU  (A-V-3) 

1  Q      ^  UUb 


and 

CINT  =  P.  (A,J,E)  +   /  \   ^(E/U)dU.  (A-V-4) 

0  u  * 

The  integrals  in  (A-V-39-*0  are  functions  of  the  Fresnel 
integrals  as  follows: 


CC  =   /  -i.  COS(E/U)dU  =  M-(k  ~    C(E/A)]        (A-V-5) 
n  \V'i  V  *    <L 


and 


S3   -      /   Ar/SIN(E/U)dU   =   Mi\  -  S(E/A)]  (A-V-6) 

o    u  /z  Y  E   2 
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where 


C(E/A)      -         / 


E/A 


'2ttT 


COS(T)dT  (A-V-7) 


and 


5(E/A)      =        / 


E/A 


0  fz¥F 


SIN(T)dT.  (A-V-8) 


That  is,  SINT  and  CINT  for  any  J  may  be  reduced  by 
integration  by  parts  to  functions  of  the  basic  Fresnel 
integrals. 

Integration  by  parts  provides  the  following  relations 
for  (A-V-l). 
If  J  =  1, 


CINT  =   /  -~  COS(E/U)dU. 
•0   UV* 


.(A-V-9) 


For  J  =  2, 


A 


CINT  =  2[A^C0S(E/A)  -  E  /  -%SIN(E/b )dU] .  (A-V-10) 


0   U 


% 


And  for  J  >  3, 
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CINT    =    (J"2>!      f«* 


(2J-3) 


rP 


L  N=l  ^J=2N7T 

[2~~'l3rc),T    „M      .n    .    1vN-l04N-i.J-2N-l_2N-l: 

cTKT/'T?/a>1    v  L2(J-2N-1)  J  !  (-1)         2  A h 

-    bllUE/A;^      (J-2N-1)! 


+  INTf  (A-V-ll) 


J 


where 


INT  =  22J  3EJ  1(-1)[2]  /  Aj  COS(E/U)dU.  (-1)J  <  0 

0   U^ 


(A-V-12a) 


or 


INT  ■  22J"3EJ  1(-1)(2)  /  -~  SIN(E/U)dU,  (-1)J  >  0. 

0  it2 


(A-V-12b) 


Similarly,  it  may  be  shown  that  for  (A-V-2)  one  obtains 
For  J  =  1, 


A  1 
SINT  =   /  A?,  SIN(E/U)dU,  (A-V-13) 

0  u3/2 


if  J  «  2, 


1  A 

SINT   =    2[A2SIN(E/A)    +   E   /   \  C0S(E/U)dU] ,  (A-V-14) 

0      U/2 
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and  for  J  >  3, 


SINT=  (J"2)!   ^ 


"T2J 


r        [24"'1]r  ,   ohm,.  1vN-l04N-3..J-2N_2(N-l) 
•|SIN(E/A)  E     IiLJ-Z?N)I]i^J_^^TrA E 

+COS(E/A)'Z     r2(J-2N-l)]!(-l)    2A E 

N=1  (J-2N-1)! 


+  INT 


J 


(A-V-15) 


where 


INT  =  22J"3EJ~1(-1)  2   /  i  SIN(E/U)dU,   (~1)J  <  0 


0   U 


(A-V~l6a) 


or 


C--1)   A 
INT  =  2^J  JE    (-1)        /  -~  COS(E/U)dU,   (-1)   >  0. 

0   U  z 

(A-V-l6b) 

The  notation  [ARG]  indicates  the  sign  of  ARG  with  the 
largest  integer  less  than  the  absolute  value  of  ARG.   The 
use  of  [ARG  ±  0.1]  ensures  the  proper  value  is  selected  and 
excludes  the  possibility  of  computer  round  off  error. 

STEP  ANALYSIS 

(Refer  to  the  subroutine  at  the  end  of  this  Appendix.) 

The  arguments  of  Subroutine  RED  include: 

A  =   Upper  limit  of  the  correction  integral. 
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E  =  KR2/2 

J  =  Related  to  the  summation  index  in  (111-30,-31) 

S  =  SIN(E/A) 

C  =  C0S(E/A) 

Al 
SS      =        /   -i-  SIN(E/U)dU 

0      U% 

Al 
CC      =         f   ~  COS(E/U)dU 

0  u3/2 

E  is  set  in  MAIN  and  carried  in  common  to  WSUM  and  as 
an  argument  through  COR  to  RED.   A,S,C,SS,  and  CC  are 
formed  in  Subroutine  WSUM  and  are  carried  as  arguments  to 
RED  via  COR.   J  is  set  by  COR.   The  values  of  (A-V-1,-2) 
are  returned  to  Subroutine  COR  on  SIMT  and  CINT. 

A  test  is  made  to  determine  if  J  <  3-   If  it  is,  the 
routine  branches  to  statement  1000. 

If  J  >  3,  equations  (A-V-ll)  and  (A-V-15)  apply.   II 
and  12  are  the  summation  limits  where  INT(ARG)  =  [ARG] 
described  previously.   The  primary  series  DO  LOOP  computes 
the  first  summation  in  (A-V-ll,  -15).   The  secondary  series 
DO  LOOP  computes  the  second  summation.   A  test  is  made  to 
check  the  criteria  of  (A-V-12,-l6)  and  proper  branching  is 
initiated.   Statement  31  computes  (A-V-12a,-l6a)  for 
(-1)J  <  0,  and  statement  32  computes  (A-V-12b ,-l6b)  for 
(-1)   >  0.   Finally,  a  statement  *J0  combines  the  terms  to 
form  (A-V-ll, -15)  on  SINT  and  CINT. 
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For  the  case  J  <  3, -a  test  is  made  to  further  restrict 
the  case  to  J  =  2  or  J  =  1.   For  J  =  2,  (A-V-10,-l4)  are 
used.   For  J  =  1,  (A-V-9,-13)  are  computed. 
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A -VI .   SUBROUTINE  IFAC 

Function  subroutine  IFAC  returns  the  value  of  the 
argument  factorialed,  e.g.,  IFAC(5)  =  5  •  ^  •  3  *  2  •  1 .   It  has 
provisions  for  IFAC(O)  =  1.   A  copy  of  the  subroutine 
follows . 
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APPENDIX  B 
STEADY  STATE  POTENTIAL  SOLUTION 

B-1'       COLLOCATION  TECHNIQUE 

In  order  to  investigate  the  perturbation  velocity 
potential  along  a  cascade  blade  for  the  steady  state  and 
7T  phase  shift  in  adjacent  blades,  one  may  use  the  following: 


*v.n(X,0)  -  -±r\    /"JL_dS  +    /  -z^EXPC-^^-^dS 


K:=0      /7T  l-i  /x=T 


^  N=o\IN-l   /x^         MX-^ 
X    N     7 


The  second  and  third  integrals  are  of  the  form 

X  N  2 

I  =    /  ■— — EXP[-B/(X-S)]  dS,   where  B  =  XR  /4 . 

-1   V^S"  (B-I-l) 

Letting 

U  =  B/(X-S),  (B-I-2) 

one  may  write  (B-I-l)  as 


I  .  /B     /  -t>  (X-B/U)N  e~U  dU  (B-I-3) 

B    UV2 
X+l 
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or 


I  =  /i1  XN    /  -\(1-B/XU)N  e"U 
B    U 
X+l 


«.J 


The  binomial  expansion  provides 


(B-I-4) 


n-B/xu)N  -  ?  ("1)K  N! 


(B/XU) 


K 


(B-I-5) 


Hence  (B-I-4)  becomes 


T  _  fa   y\l\      V   ("1)KBK  r°°TT-(K+l)    ~U   ,„    ,-  -  ,s 

I   -  </B   X  N!    E    (N_KyTXfXK  '  e        dU*    (B-i-6) 


K=0 


B 
X+l 


It  may  be  shown  that  integration  by  parts  provides 


M 


,      °°  -x 

1     f-$—   dX,   M  >  0 
1   0  Ym-1    ' 


(B-I-7) 


where  M  may  be  any  half  odd  integer  greater  than  zero, 

e.g.,  M  =  1/2,  3/2,  5/2,  etc.   Continuation  of  the  integration 

process  yields  a  general  relationship 


00  -x 

/  ^r— dX 


r 


,2   X 


m 


P2  ™~2  (M-I-DK-l)1"1 


(M-l)!(F2)m_I 


] 


1 
(-l)m"2"^P  ERFC(F) 
+  (M-l) ! 


(B-I-8) 
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where 


,     ».  -U       ., 
ERPC(F)  =  —   /  £_-  dU  «  ~-[l  -  erf(F)].       (B-I-9) 


For  the  case  N  =  0,  (B-I-6)  is 


1 "  ^i  S  dU-  (B-I-10) 

X+1U 
The  last  integral  in  (111-47)  may  be  written 

/  -A-  dS  =  2/X+T  £  ^,;   ,!    1}    ,       (B-I-ll) 
-1  /x=S*  K=0(2K-l-l)K!(N-K)! 

and  observe  that  for  the  case  N  =  0,  (B-I-ll)  becomes 


x   1 
/   _i__  dS  =  2/X+T.  (B-I-12) 

-1   /T-^S1 


Combining  Eqns .  (B-I-6, -10,-11,  and  -12)  in  (111-47) 
and  replacing  N  with  J-l  and  K  with  K-l,  one  arrives  at 
the  final  computer  compatible  equation 


T|»„n(X,0)    =   -^-[2/x+7  +    2/B1  SSUM[-3/2,    b/(X+l)] 


t        r    2(-l)K"1(J-l)!XJ~Kr^  (X+l)K"1/xTT 

X0,J-1 


1(K-1)!(J-K)!  1  0.J-1  2K-1 


"   Vl   J-!8*"""^  SSUM[-(K+|),B/(X+1)]}]  (B-I-13) 
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where 


SSUM(A,FSQ)  -  i     /  XAEXP(-X)dX.         (B-I-14) 

FSQ 


SSUM  is  computed  in  Subroutine  SSUM  and  will  be  dealt  with 
in  Appendix  B-I-B. 

B-I-A.   THE  MAIN  PROGRAM 

The  purpose  of  MAIN  is  the  computation  of  (B-I-13). 
Inputs  to  the  program  are  the  controlling  parameters,  blade 
spacing,  R,  and  thickness  parameter,   X,  as  well  as  the 
values  of  the  interference  coefficients  obtained  from  the 
collocation  coefficient  solution  of  Appendix  A. 

Step  Analysis 

(Refer  to  the  routine  at  the  end  of  Appendix  B-I.) 
Double  precision  is  used  in  Subroutine  SSUM  and  single 

precision  is  used  in  MAIN.   Double  precision  arguments  to 

and  from  SSUM  are  carried  on  variables  beginning  with  W. 
The  inputs  DLAM  =  A,  DK  =  K,  and  RKO  =  R,  and  the 

coefficients  of  the  upper  and  lower  interference  functions 

(P,   and  PA  ,  m  =  1,2,..., 10)  are  read  into  the  routine  from 
lm      0m*       '  i         3 

data  cards  and  printed  out  immediately  for  verification. 
(Note  that  K,  the  reduced  frequency,  is  used  for  accounting 
purposes  on  the  output  only.)   Array  DEL(20)  carries  the 
coefficients  with  the  lower  blade  coefficients  in  the  first 
ten  positions. 
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The  inputs  positioned,  various  fixed  parameters  are 
formed  and  the  station,  X,  is  initialized. 
B   =   AR2/4 
BH  =   /F 

c  =  V/iFX 
D   =   -3/2 
X   =   -0.9 


DO  LOOP  100  is  the  station  control  loop  which  generates 
a  potential  value  at  each  of  20  positions  on  the  blade 
(X  =  -0.9,  -0.8,. ..,  1.0). 

In  preparation  for  the  use  of  Subroutine  SSUM,  the  lower 
limit  is  formed, 

WPSQ  =  B/(X+1). 

SSUM  is  called  to  form 

SSUM[D,B/(X+1)]  =  i    /  ?—   dU. 

2   B    \J3/z 
X+l 

With  AA  =  /X+l" j  the  first  two  terms  in  (B-I-13)  are  formed: 

BB  =  2/X+T  +  2/B1  SSUM[-3/2,B/(X+l)]. 

DO  LOOP  200  accomplishes  the  J  summation  in  (B-I-13) 
while  LOOP  300  carries  out  the  K  summation.  (Recall  that 


1H5 


L  =  10.)   Each  term  in  the  summation  is  built  up  in  a 
straightforward  manner.   Subroutine  SSUM  is  required  to 
compute 


where 


SSUM(A,WFSQ)  =  \   WFSQ;  U   e   dU> 


A  =  -(K+|). 


Function  subroutine  IFAC  is  also  used.   This  subroutine 
produces  the  factorials  of  its  argument  and  is  explained 
in  Appendix  A-VI . 

The  summation  loop  completed,  the  summation  term  is 
combined  with  BB  and  C  to  form  POT  =  ^Ks_q(X,0)  of  equation 
(B-I-13).   Finally  POT  is  printed  together  with  its  station. 

B-I-B.   SUBROUTINE  SSUM 

The  purpose  of  SSUM  is  to  produce  the  value  of  the 
general  integral 

00 
\        !    XAEXP(-X)dX,   A  =  ±1/2,  ±3/2,  ±5/2, (B-I-15) 

2 

F^ 

Only  the  derivation  of  the  procedure  for  A  <  0  will  be 
presented  in  detail  here.   The  case  for  A  >  0  is  derived  in 
a  similar  manner.   This  latter  case  is  not  of  interest  in 
the  potential  solution. 
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The  on  line  function  subroutine,  INT(A),  generates  the 
sign  of  A  with  the  largest  integer  less  than  A,  e.g., 

INT (-3/2)  =  -1.   The  subroutine,  ERFC(F),  produces  the 

p 
compliment  of  the  error  function  using  P   as  the  lower  limit. 

It  should  be  understood  that  (N/2)!  =  |  •  (|  -  1) •  (|  -  2)  • . .  | 

and  that  (-1/2) !  =  1. 

Now  for  demonstration  purposes,  (B-I-8)  is  expanded 
for  the  cases  A  =  -3/2  and  A  =  -5/2. 
A  =  -  3/2 

4   /  ~T   dX  =  4[-~-  ~  -2v/7ERPC(F)]  (B-I-16) 

2  ^2  %Vz  2  ^p2   * 

A  =  -  5/2 


dY   -    X     T     1     f    2  22  }   +    22^ERFC(F)1 

dX   -    2     L  "72/^3"  "   "STlFJ    +  T7L         J 


-,       co  -X 

-     /  — 

2  P2   XV2  PF^(3F" 

*  e  (B-I-17) 


A  routine  may  be  generated  in  the  following  manner. 
Observe  that  if  one  defines  I  =  INT(A),  there  are  I  terms 
in  the  series  (Those  terms  not  including  the  error  function 
term.).   The  numerator  is  increased  by  two  in  each  succeeding 
term.   The  denominator  of  the  first  term  includes  (21-1); 
each  succeeding  denominator  in  the  series  is  multiplied  in 
turn  by  (2I-D-2,  (2I-D-4,  ....   The  power  of  F  in  the 
denominator  begins  with  (21-1)  and  is  decreased  by  two  in 
each  succeeding  term.   The  sign  changes  with  each  term.   Also 
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note  that  the  numerical^ coefficient  of  the  error  function 
term  is  the  same  as  that  for  the  last  series  term  with  the 
sign  reversed. 

Step  Analysis 

(Refer  to  the  subroutine  at  the  end  of  Appendix  B-I.) 

p 

The  inputs  to  the  subroutine  are  A  and  FSQ  =  F  .   The 

value  of  (B-I-15)  is  returned  on  SUM.   In  the  first  steps 
various  fixed  parameters  are  formed,  the  summation  parameter 
is  zeroed,  and  INT(A)  is  determined.   A  test  is  made  to 
determine  if  A  >  0  or  A  <_  0 .   Here  only  the  branch  A  <_  0 
is  of  interest  although  the  structure  of  the  unused  portion 
of  the  subroutine  for  A  >  0  is  similar. 

For  A  <_  0,  statement  10  tests  for  the  special  case 
A  = ~   (or  1=0).   This  case  cannot  be  generated  by  the 
algorithm  used  for  other  cases.   If  A  =  -  «  ,,  statements  12 
and  1*1  set 


SUM  =  ^~   ERFC(F) 


and  the  routine  returns  to  MAIN. 

If  A  <  — p,  statement  13  begins  the  algorithm.   After 
several  preliminary  steps,  DO  LOOP  11  forms  the  series 
pointed  out  previously.   The  method  takes  advantage  of  the 
observations  made  in  the  introductory  material.   The  series 
formed  on  SUM,  it  is  divided  by  the  exponential  and  the 
error  function  term  is  added.. 

The  result  is  returned  to  MAIN. 
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B-II.   SONIC  WIND  TUNNEL  WALL  INTERFERENCE  TECHNIQUE 

Development  of  the  Sandeman  transonic  wind  tunnel  wall 
interference  solution  for  the  steady  case  resulted  in  the 
following  perturbation  potential  function  from  Section  III-E 


!f»(X,0)  =  - 


] 


v'ttX 


/ 

-1 


X 


v(S) 

/x-s' 


dS 


v(S) 


2  2 
N  d 


+  2  I        [      -^f  EXP[-X  ^3  dS   . 
N=l  -1  VX-S  A~h 


(111-60) 


To  compare  this  solution  with  the  cascade  problem,  (d)  is 
related  to  blade  spacing  R  by  d  =  R/2  [Refer  to  Fig.  7]. 
Further,  as  shown  in  Eqn .  (III-18), 


v(S)  =  \\>    (X,Y  =  0)  =  -1,   for  K  =  0. 


It  follows  that  the  potential  may  be  written: 


<KX,0)    = 


/ttT 


X 


/ 

■1     v/X-S 


dS 


N=l      -1        /X-S 


dS 


(B-II-1) 


Integrating  the  first  term  and  making  the  proper  substitution 
in  the  second  term  one  may  write 


i|/(X,0)  = 


2/x+T  +  2  i 
N=l 


XN2R2 


-U 


2J2   ^dI 
4(X+1) 


(B-II-2) 
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The  basic  integral  within  the  summation  may  be  expressed  in 
terms  of  the  error  function,  i.e., 


oo    _u 
e 


/   tj  ™ 


,2   U 


2  e 


-F' 


P 


-  2VtTeRPC(P) 


(B-II-3) 


where 


„2  _  XN2R2 


(B-II-4) 


and 


ERPC(P)  =  1  -  erf(P) 


(B-II-5) 


Hence,  the  potential  function 


*(X,0)  = 


?  2 

N=l 


—  X/tt 


r2^2 


XN2R2 


AN^   ERFCC^^-y])  • 


(B-II-6) 


Computation  of  this  function  is  easily  carried  out  by 
computer.   The  series  is  truncated  when  the  inclusion  of 
additional  terms  only  changes  the  sum  by  a  value  less  than 
the  convergence  criterion,  EPS. 

Step  Analysis 

(Refer  to  the  routine  at  the  end  of  this  appendix.) 
The  inputs  to  the  routine  are  DLAM  =  X  and  RHO  =  R. 
DK  =  K  is  inputed  for  accouting  purposes  in  the  output  only 
The  output  includes  the  above  parameters  and  the  value  of 
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the  perturbation  velocity  potential  at  each  of  twenty 
stations  along  the  blade. 

The  control  parameters  and  various  fixed  parameters  are 
calculated  prior  to  entering  the  control  loops.   DO  LOOP  100 
provides  control  to  calculate  the  potential  at  each  of  the 
stations  beginning  at  X  =  -0.9  and  increasing  by  increments 
of  AX  =0.2.   DO  LOOP  50  computes  the  series  summing  terms 
until  convergence  is  reached.   Again,  the  function  subrou- 
tine ERFC  is  a  standard  on  line  library  routine  performing 
the  operation  in  (B-II-5). 

Once  the  series  is  computed,  statement  51  forms  the 
potential  on  POT  which  is  subsequently  printed  along  with 
the  applicable  station  location. 
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